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Preface 


The area of science to which this book is devoted 
is very young. Its basic ideas were formulated 
as recently as 1957 in a paper by two British 
scientists, S. R. Broadbent and J. M. Hammers- 
ley. During the mid-1950’s Broadbent was work- 
ing at the British Coal Utilization Research 
Association on the design of gas masks for use 
in coal mines. He came across an interesting 
problem and presented it to the mathematician 
Hammersley. 

The principal element of a mask is filled with 
carbon granules through which the gas must 
flow. Carbon contains pores that are connected 
together in an intricate manner to form a sort of 
complicated maze. A gas can enter the pores by 
being adsorbed on their inner surface. It was 
found that if the pores are wide and well con- 
nected, the gas penetrates deep into the carbon 
filter. Otherwise the gas cannot get beyond the 
outer surface of the carbon. The motion of a gas 
through the maze is a new type of process, which 
differs from diffusion. 

Broadbent and Hammersley called the phe- 
nomenon percolation, and the theory underlying 
processes of this type is referred to as percolation 
theory, 
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Since Broadbent and Hammersley published 
their pioneer paper 28 years ago, it was discovered 
that percolation theory can be used lo interpret 
an exceptionally wide variety of physical and 
chemical phenomena. The electrical properties 
of disordered systems, such as amorphous semi- 
conductors, crystalline semiconductors with im- 
purities, and materials formed as mixtures of 
a diclectric and a metal, are probably the best 
understood application of percolation theory. 

The phenomena best described by percolation 
theory are critical phenomena. They are charac- 
terized by a critical point at which some of the prop- 
erties of the system undergo abrupt changes. 
Critical phenomena include second-order phase 
transitions (e.g. the transition of a metal from 
its normal to its superconducting phase when its 
temperature is lowered). The physics of all 
critical phenomena is very unusual, but there 
are some common features, the most important 
of which is that in the neighborhood of the critical 
point, the system appears to break into blocks 
which differ in their properties, with the size 
of the individual blocks growing until the system 
approaches the critical point. The blocks are 
quite randomly shaped. In some phenomena the 
whole configuration changes chaotically because 
of thermal motion, while in other phenomena the 
configuration may be frozen in time but changes 
from specimen to specimen. The blocks are in 
complete disorder, so that no regularity is dis- 
cernable from an instantaneous photograph. How- 
ever, this geometry, which can be called the 


geometry of disorder, has quite definite properties 
“on the average”. 
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Actually, geometry is inseparable from physical 
properties. For instance, the physical properties 
of a crystal are determined by the geometry of 
its lattice. Likewise, the “geometry of disorder” 
determines a number of properties of a system 
in the vicinity of a critical point. The most 
interesting feature is that owing to the large size 
of the blocks the geometry is virtually independ- 
ent of the atomic structure of the material and 
thus possesses properties common to a number 
of quite dissimilar systems; hence, the univer- 
sality of the physical properties that we find in 
the neighborhood of critical points. 

This type of relation between physics and geo- 
metry can be traced in percolation theory, and 
in fact, making this relationship clear is the 
main objective of this book. Percolation theory 
is formulated in terms ofsimple geometric images, 
such as wire nets, spheres or crystal lattices. 
The theory does not operate wilh the concept of 
temperature, and this makes it possible to clarify 
the idea of critical phenomena for readers not 
familiar with statistical physics. 

Percolation theory, as a theory of critical 
phenomena, is not yet a mathematically rigorous 
science. A large number of important proposi- 
tions have not yet been proved, and certain 
questions have not been answered. However, if 
a rigorous proof exists but appears to be com- 
plicated, T decided to replace it here by argu- 
ments that will explain the result rather than 
prove it. However, an effort was made to separate 
those propositions that have not yet been proved 
from those that have. 

The book offers a detailed presentation of the 
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theory of percolation and ils various applications. 
However, a definition of what is understood by 
percolation theory and what processes it de- 
scribes is deferred until the last page of the book. 
The definition must encompass so many complex 
concepts that to give it at the beginning would 
be meaningless. Nearly every chapter concen- 
trates on one specific problem whose analysis 
leads to a problem in percolation theory. It is 
assumed that, having read several chapters, the 
readers will become aware of the common ground 
the various percolation problems have, and also 
become aware of the relevance of the title of the 
book. 

As a rule, the problems 1 chose represent im- 
portant applications of percolation theory. How- 
ever, some (such as laying out an orchard in 
Chapter 5 and the propagation of rumors in 
Chapter 44) are illustrative or even slightly 
humorous. ot 

The fundamentals of elementary probability 
theory required for understanding the material 
are also covered, Chapter 1 gives a general notion 
of probability and random variables, Chapter 2 


introduces the rules for the addition and multi- 
plication of probabilities and defines the distri- 
bution function, An “easy reading” version of the 
book is obtained by skipping Chapter 2 and the 
other chapters and sections marked with two 
asterisks, although the reader will not then be 
able to follow the derivation of certain of the 
quantitative resulls or a few of the exercises. 
Nevertheless, this should not hamper the under- 
standing (even if slightly Jess comprehensive) 
of the other chaplers. 
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| consider the exercises in the text important. 
As a rule, the exercises are quite simple, and 
I recommend that the reader do them without 
looking at the answers section beforehand unless 
this is specifically adv 

B. 1. Shklovskii played an important role in 
the creation of this book because together we dis- 
cussed ils structure and title, and then he read 
the manuscript. I am extremely grateful to him 
for his contribution. 

I am also grateful to my colleagues L. G. Asla- 
mazoy, N. B. Vasilyev, Yu. F. Berkovskaya, and 
M. E. Raikh because they read the manuscript 
and suggested a number of useful changes. 

I am especially grateful to my wife N. I. Efros 
for carrying the heavy burden of preparing the 
manuscript for publication, 


A. Efros 


Part I 


Site Percolation Problem 


Chapter 4 


Percolation Threshold 


Two Pundits Shred a Wire Mesh 


It is not frequent nowadays for a scientific journal 
to publish a report on experiments done with, 
for example, a piece of ordinary wire mesh bought 
in the nearest hardware store. A paper written 
by two American physicists Watson and Leath, 
published in the Physical Review in 1974, was 
definitely not the first one in the realm of per- 
colation theory; nevertheless, it is with this study 
that we start our story. 

The piece of wire mesh used by Watson and 
Leath was a square containing 137 x 137 = 
= 18 769 nodes, or sites, with the neighbor sites 
separated by a spacing of 1/4 inch = 6.35 mm. 
The scientists soldered copper electrodes to two 
opposite sides of the square and connected the 
mesh to an electric resistance-measuring circuit 
(Fig. 1a). Then they started to block some sites 
in the mesh and monitored the electric resistance 
as a function of the fraction of blocked sites. As 
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b) (c) 


Fig. 1. Schematic of the experiment of Watson and Leath: 
(a) an initial wire mesh (the number of sites in the figure 
is greatly reduced); (b) a piece of wire mesh with blocked 
sites (blocked sites are shown by black circles, and non- 
blocked sites by open (white) circles); (c) a black site sig- 
nifies that the four wires coming to the site are not in 
contact, and a white site signifies intact contacts (no elec- 
tric current flows through black sites in any direction, 
while it flows through white sites in any direction). 


16 Part 1. Site Percolation Problem 


shown in Fig. 1b, c, a site was blocked by simply 
cutting all four wires joining at this site. 

The choice of a new site to be blocked was made 
randomly among those still in the mesh. In 
principle, this could be done by writing the 
coordinates of cach sile on separate slips of 
paper, pulling all these slips into a hat, stirring 
well the contents, and then extracting them one 
by one. However, such a procedure, as well as 
any other mechanical method of sequence allot- 
ment, is highly impractical if the number of sites 
is large; consequently, the scientists used a ran- 
dom sequence of site coordinates generated in 
a computer, Later we shall describe how a com- 
puter can be “made to generate” random numbers, 
but for the time being we can think that a com- 
puter was replaced, with no loss to clarity, by 
a hat-and-paper-slips combination. 

Obviously, the electric conductance of the 
mesh decreased with increasing number of blocked 
sites. (Electric conductance is defined as 
a quantity reciprocal of resistance, Resistance 
is measured in ohms, and electric conductance 
in inverse ohms (ohm-'),) Furthermore, if < 
denotes the ratio of nonblocked sites to the total 
number of sites (137), the electric conductance 
vanished at a certain value of x that hereafter 
we call the threshold (critical) value, or percolation 
threshold, and denote by to. This vanishing 
occurred when the last path connecting the left- 
and right-hand electrodes had been cut. One of 
the objectives of the experiment was to deter- 
mine z. It was found that ze = 0.59. 

Probably, the first question that must be 
answered is whether the variable zg is random, 
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irreproducible from one experiment to another 
or ze is quite definite. Let us assume that the 
experiment is repeated with another square piece 
of wire mesh and with a different random sequence 
of blocked sites. It is common sense to expect 
that the configuration of blocked and intact 
sites being quite different al each stage of the 
second experiment in comparison with the first, 
the cutting of the last path connecting the elec- 
trodes must also occur at a different value of z, 
so that the second experiment must give a new 2. 
This is certainly correct. 

The threshold value xe in the experiment we 
are now discussing is a random quantity, or 
variable. Quantities of this sort will be encoun- 
tered throughout the book, and so it is useful 
to answer from the very beginning the following 
question: 


What Is a Random Variable? 


In mathematics a variable is said to be random 
if the values it takes and the frequency at which 
it takes them are known, but what particular 
value the variable will assume in each particular 
case is unknown (and cannot be known in the 
framework of the given mathematical prob- 
lem). 

Here is a classical example of a random vari- 
able: a cube (a die) with numbered faces is being 
thrown onto a table. The numeral read off the 
topmost face is a random variable. Such a vari- 
able is said to be discrete because it assumes only 
certain discrete values (in the chosen example 
these are six numerals: 1, 2, 3, 4, 5, 6). It would 


2-0240 
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be impossible to predict the specific numeral 
obtained in each given trial (i.e. in each throw), 
but the probability of getting a specific number 
(e.g. 4) is predictable. Assume that the number 
of trials was Q, and face 4 appeared in Q, cases. 
The ratio Q,/Q is called the relative frequency of 
recording a given value of the random variable 
(here, the numeral 4). If the total number of 
trials is not very large, this ratio varies: in 
another series of Q trials the ratio Q,/Q may be 
substantially different. However, as the number 
of trials Q increases, these fluctuations progres- 
sively diminish. The relative frequency of record- 
ing a given value of the random variable ap- 
proaches the limit which is called the probability 
of this value. 

Let us denote the probability of obtaining 
face 4 by P (4). If the cube is “honest” (nol 
“loaded”), that is, if all its faces are equivalent, 
the quantity P (4) is easily predicted. On the 
average, any one of the six faces of the cube will 
turn up topmost the same number of times, so 
that Q,/Q = Q,/Q =... = 1/6 provided Q is 
large. Therefore, the probabilities of recording 
the six numerals are identical and equal to 4/6. 

We thus find that after a large number of 
throws, chance melts into the background and 
clears the stage for a regularity, namely, the 
symmetry of cube’s faces. 


Mean Value and Variance 


Let us return to the experiment with the wire 
mesh. We concluded that since the experiment 
operated with a random sequence of blocked 
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siles, the critical concentration ze at which the 
current between the left- and right-hand elec- 
trodes is interrupted is also a random variable, 
and it would be impossible to predict the value 
of ze in each specific experiment. 

The theoretical approach to the situation 
could be a study of the “mean” properties of the 
quantity Te, that is, of the properties apparent 
after a sufficiently large number of experiments 
run under identical conditions. These conditions 
are, first, the total number of sites in the wire 
mesh, Æ (A = 137? in the experiment described 
above), and second, the properties of the generator 
of random numbers that prescribes a random se- 
quence of blocked sites. The requirement that the 
properties of the generator should not change 
from one experiment to another does not mean 
at all that the sequences of blocked sites must be 
identical. (In this case all values of z_ would be 
identical!) All we need is to use the same method 
of generating the random sequence of sites to be 
blocked in all experiments (e.g. a hat with slips 
of paper). 

Having conducted Q experiments with a wire 
mesh containing f° sites, we obtain Q values zi, 
where the subscript i is the number of the exper- 
imental run. For example, zı, denotes the ze 
obtained in the fifteenth series of measurements. 
The most important of the mean values is the 
arithmetic mean, or the average, zg, obtained 
by adding up all z; and dividing the sum by 
the number Q of experiments: 


To= ateit (1) 


Qe 
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The bar over aand the subscript Q denote averag- 
ing over the results of Q experiments. The quanti- 
ty XQ is still a random variable. If we conduct 
another series of Q experiments under the same 
conditions and calculate the resultant Xe, this 
Tå will be somewhat different. Ilowever, the 
larger the number of experiments, Q, within the 
series, the smaller the differences between the 
averages representing different series. The point 
is that random fluctuations of x; cancel each other 
in the long run, so that as Q increases, the arithmetic 
mean ig tends to a quite definite value independent 
of Q but dependent on the conditions under which 
the experiments were run. This limiting value is 
called the mean value of a random variable. (In 
probability theory this limiting value is also 
called the expectation of a random variable, but 
this term will not be used in this book.) 

The mean value of the percolation threshold 
of a wire mesh made up of . sites will be denoted 
by ze (S°). The quantity ze (S°) is not a random 
variable but a certain quantity. Its dependence 
on fis a regular property worth thinking 
about. 

An important characteristic of a random vari- 
able ze is also found in the deviations ô; of the 
values z; from the mean value: 


6; = t; — Te (S`) (2) 


The deviations ô; vary from one experiment to 
another, and we want to select a quantity char- 
acterizing the properties of ô; “on the average”. We 
cannot select the arithmetic mean for this quantity 
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because it tends to zero when Q —> œ. Indeed, 


Sitôt... tõe _ 21+22+...-} 29 > 
a = 0 Te (M) 
But the greater Q is, the less the difference is 
between the first term on the right-hand side 
of this equality and the second term; this proves 
the proposilion made above. This result is ob- 
tained because the individual values of x; neces- 
sarily fall on both sides of the mean value, so 
that on the average the deviations cancel out. 
We could also take the arithmetic mean of the 
nonnegative quantity |6;|; however, the stand- 
ard procedure is to calculate the variance 8° (S°) 
which is the arithmetic mean of the squares 
of deviations (for Q — co) which, of course, are 
also nonnegative quantities: 


EN) 834-52 opi 8 (3) | 


The quantity 6 (4°) = [6 (4C) is called the 
root-mean-square deviation of a random variable, 
or rms deviation. It is the rms deviation ô (f°) 
that characterizes the typical deviation of the 
values z; from their mean value ze (°). Obvi- 
ously, the quantity 6(.f) is also a function of 
the total number of sites, ., in the wire mesh. 

Strictly speaking, x, is a discrete quantity 
because it is obtained by dividing the number 
of nonblocked sites by the total number of sites 
f°, and therefore, it assumes only those values 
that convert to integers when multiplied by f% 
Let us denote all possible distinct values of/ (die 
random variable, Zo hy win rm ARY 


Dai fe, oF 
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The mean value z, (#*) can be expressed in 
terms of the probability P (tx) of the random 
variable £e assuming the value zp. Remember 
that expression (1) contains the sum of z; obtained 
in Q experiments. Each value of z; can appear 
many times. Formula (1) can be written in the 
form 


Ig= oo (4) 


where the summation is carried out over all 
distinct values x, that Ze Can assume (no one value 
of zp is encountered twice in this sum!), The 
factor Q, is equal to the number of times the 
value z, turned up in a series of Q experiments. 

The quantity Q,/Q is a relative frequency of 
obtaining Try as a result. When Q is very high, 
this ratio becomes equal to the probability 
P (x). By definition, for large values of Q the 
left-hand side of expression (4) turns into re (SP) 
Consequently, 


Zo (W) = 2P (2) + 2P (e) + .., (5) 


i.e. the mean value equals the sum of all the 
values that a random variable can assume, 
multiplied by their probabilities, Likewise, 


8 (SM) = (t1 — xe (9)? P (x) 

F (£a — te (W) P (a) +... (6) 
The summation in formulas (5) and (6) is carried 
over all possible values that the random vari- 


able zp can assume, with each value encountered 
only once in the sum, 
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By definition, Q, + Qe +... =Q, so that 
the definition of P (xp) yields 
P (e) +P (y+... =1 (7) 


The sum of the probabilities of all the values 
that a random variable can assume equals unity. 


Why a Large Wire Mesh? 


It was very simple to calculate, in the problem 
of “honest” cube, the probability for the random 
variable to assume a specific value. The properties 
of the random variable x, are incomparably more 
complex. 

It will be shown at the end of the next chapter 
how to solve the problem for a square network of 
four nodes (sites) (2 X 2 sites, f° = 4). The 
result of the solution: the random variable 
ze can assume only two values, namely, 41/2 
and 4/4. The first of these is assumed with prob- 
ability P (1/2) = 2/3, and the second with 
probability P (1/4) = 1/3. According to for- 
mulas (5) and (6) (where the sums consist of only 
two terms), ze (4) = 5/12 and ô (4) = V 2/12. 

With sufficient patience, we can also solve 
the problem with 3 x 3 sites (4° = 9). Experi- 
ence shows that the required efforts increase 
enormously if the side of the square is augmented 
by only one site. At the same time, the networks 
whose properties are of especially high interest 
contain a very large number of sites (e.g. 10"). 
Such networks can serve as models of films con- 
sisting of atoms. Indeed, the spacings between 
the atoms in condensed materials (liquids and 
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crystals) are as a rule about 3-10-8 em. Conse- 
quently, a film covering 1 cm? witha thickness of a 
single layer of atoms consists of roughly 
10 atoms. 

The problem of finding the probability for the 
percolation threshold of a network made up 
of a very large number of sites, f°, lo assume 
a specific value is the central problem of per- 
colation theory. It will be discussed, in one form 
or another, throughout the book. Here we want 
to single out, giving practically no proof, the 
most important property of this problem that 
offers the key to understanding the problem as 
a whole: 


_ The root-mean-square deviation ô (4) dimin- 
ishes with increasing number’ 4^ of sites in power- 
law fashion, tending to zero as WS —> œ. 


This property is expressed by the formula 
a Cc 
ôN LETE (8) 


where C x 0.54 and v x 1.3. (The quantity v 
is called the exponent of correlation radius. It is 
discussed in detail in Part III of the book.) 

Watson and Leath could not arrive at for- 
mula (8) as a result of their experiment. To 
obtain (8), il was necessary lo use wire meshes 
with different values of f° and to run numerous 
experiments at the same f°. Besides, expression 
(8) is a result of the theoretical studics to be 
discussed in Part III. 

As follows from formula (8), the results of 
experiments with different random sequences of 
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blocked sites differ from one another the less, 
the greater the number of sites in the wire mesh. 

Why should this be so? The point is that a 
sufficiently large wire mesh manifests all or 
nearly all the possible configurations of intact 
and blocked sites. These configurations appear 
to change places in different experimental runs. 
Consequently, the role of chance is the smaller, 
the greater ° is. An infinite net contains an in- 
finilely large number of large subnets, so that 
randomness ceases to play any role at all, and 
the value of xo is not a random variable but a certain 
quantity equal to 


t= lin pels) 
M +00 

It is this limiting value that is in fact defined 
as percolation threshold, and it was for finding 
this threshold that Watson and Leath have con- 
ducted their experiment. Otherwise why should 
they have taken a wire mesh with nearly 19 000 
sites? They could have taken a 2 X 2 wire mesh! 

Now let us formulate the most important 
result of this chapter: 


The concept of a sharply defined percolation 
threshold independent of the choice of a random 
sequence of blocked sites used in the experiment 
holds in an infinitely Jarge system. No sharply 
defined threshold exists in a finite system, but 
there is the so-called critical region with width of 
the order of 6 (./*) into which the values of te, 
which are recorded in the majority of experimental 
runs with different random sequences, fall. As the 
size of the system increases, this region contracts 
to a point. 
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It must be borne in mind, however, that the 
dependence on the size of the system is important 
only as long as we attempt to model a phenome- 
non artificially (e.g. by using a wire mesh). As 
a rule, percolation theory is applied to systems 
in which individual elements are extremely small 
(e.g. atoms, for as we have already mentioned 
(see Chapter 3) 1 cm? of a monatomic layer con- 
tains Æ = 10% elements, and 1 cm? contains 
MS = 10™!), and so to a very good accuracy, the 
systems can be considered infinite, and the un- 
certainty in percolation threshold due to a sys- 
tem’s size can be ignored. 

The problem that Watson and Leath were 
solving is referred to as the site percolation prob- 
lem, or simply site problem (because the random 
elements are represented by sites). Quite a few 
problems in science reduce to site percolation, 
and one of them (a doped ferromagnetic) is dis- 
cussed in Chapter 3. 

The exact value of percolation threshold for 
this problem has not been found yet. The quantity 
Te (I) is found approximately for large values 
of . in computer simulations or in so-called 
analogue experiments similar to that of Watson 
and Leath. (The techniques involved may vary 
substantially.) 

The degree to which the result obtained devi- 
ates from the sought limiting value can be evalu- 
ated from the change in ze (f°) with changing Se. 
A comparison of results obtained by different 
methods makes it possible to conclude that the 
number 0.59 is correct to within two decimal 
places (although it is far from obvious beforehand 
that the number of sites f° = 137? is sufficient). 
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Obviously, there can be no end to improving 
the accuracy of ze in the subsequent decimal 


places (see Exercise 4). 


Exercises 


1. Define a discrete random variable a as the 
numeral on the topmost face of a cube after 
a throw. Find the mean value of a. 

2. Define percolation threshold as a value of x 
at which percolation sets in not from left to 
right, but from top to bottom. Will this change 
the results of individual experimental runs, 
Lo (S), £o? Assume the network to be square. 

3. Answer the same question but take into 
account that percolation threshold is defined as 
the minimum value of z at which percolation 
exists both from left to right and from top to 
bottom. 

4. Answer the same question but in the case 
when percolation threshold is defined as the 
maximum value of x at which there is no per- 
colation both from left to right and from top 
to bottom. 

5. Make use of formula (8) and calculate the 
root-mean-square deviation in the conditions 
of the Watson-Leath experiment (4° = 137°). 
What accuracy can be expected if only one ex- 
perimental run was carried out? 

Note. In principle, the result of a single experi- 
ment may differ very much from the mean value 
Te (f°). However, using the distribution function 
given below for percolation thresholds (see for- 
mula (6) in Chapter 2), it can be proved that 
the probability for the result of a randomly se- 
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lected experiment to fall in the interval from 
Lo (NM) — ô bo ze (H) + ô is roughly 0.7. The 
greater J is, the smaller the “typical deviation” 
from the mean value is. 


Chapter 2** 


Basic Rules 
for Calculating Probabilities. 
Continuous Random Variables 


This book is devoted to the laws governing disorder, 
and widely uses the concepts of probability and 
random quantity, or variable. These concepts 
have been partly introduced in the preceding 
chapter, so that the reader who is not inclined 
to delve into the mathematical side of the prob- 
lem is welcome to skip Chapter 2, as well as 
all the subsequent chapters and sections marked 
with double asterisks. The reader, who wants 
to follow the solution of a number of beautiful 
mathematical problems given in the book and to 
form a more profound picture of percolation 
theory, must know the rules for the addition 
and multiplication of probabilities presented 
in this chapter. 


Events and Their Probabilities 


The concept of probability is used not only when 
we deal with numerical values taken up by a 
random variable. Any experiment with a random 
outcome can be discussed in these terms, The 
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distinct resulls of experiments are referred to 
as events. The relative frequency of an event is 
defined as the ratio of the number of trials that 
led to this event, to the total number of trials. The 
probability of an event is defined as the limit 
which the relative frequency of the event approaches 
when the number of trials tends to infinity. 

Example. Red, green, and blue balls in iden- 
tical numbers are put in a box. The balls are then 
stirred, and one of them is extracted in a grab-bag 
fashion. What is the probability of the event con- 
sisling in extracting a red ball? In contrast to 
the experiment with dies, here events differ not 
in quantity bul in quality (the color of balls). 
llowever, the argument runs along the same line. 
The number of balls being identical in each color, 
a red ball will be chosen in 1/3 of all trials. Gon- 
sequently, the sought probability equals 1/3. 
The probability of drawing a blue ball or a green 
ball is also 1/3. 

By definition, probability is a quantity that 
varies from zero to unity. An event with zero 
probability is, for instance, the extraction of 
a blue ball from a box with only red balls. The 
probability of drawing a red ball from the same 
box is unity. An event with probability equal to 
unity is said to be a certain event, not a random 
event. 

The concept of probability plays an enormous 
role in clarifying the regularities in the world 
of random processes. A law is quite often virtu- 
ally buried under chance. Imagine that you try 
to find a regularity in birthrates of boys and girls 
from the record of one maternity. You see a ran- 
dom sequence of the type BGBBBBGBGG.... 
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Sometimes it seems that boys appear much more 
often, and sometimes your conclusion has to be 
reversed. Your friend assures you that “only girls 
are born nowadays”. His impression may slem 
from the fact that girls were born in three families 
that he visits. 

However, a regular pattern does exist. The 
ratio of the probability for a boy to be born 
to that for a girl to be born is 51.5 to 48.5.The 
regularity is observed to hold quite well in 
such large countries as the USSR or the USA, 
even if the data cover only one year. 

In contrast to the problem with balls of differ- 
ent colors, the problem of male/female ratio in 
childbirth is very difficult to solve theoretically. 
However, the statistical data involved do reveal 
certain thoroughly analyzed properties of human 
physiology. 


Addition of Probabilities 


Events are said to be incompatible if they cannot 
be observed in the same trial. For instance, the 
event consisting in drawing a red ball is incom- 
patible with the event consisting in drawing 
a blue ball because, by the conditions of the 
problem, only one ball can be drawn in any trial: 
either red, or blue, or green. The events consist- 
ing in getting the numerals 5 and 2 in one throw 
of our cube are incompatible. 

We shall prove two important properties of 
probability. 

1. Addition rule. The probability for one (no 
matter which) event out of several incompatible 
events equals the sum of the probabilities of these 
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events. Suppose that we want to find the proba- 
bility for a cube to show either 3 or 4. The number 
of trials in which these two outcomes were 
obtained is equal to the number of trials that 
gave 3 and the number of trials that gave 4. 
By definition, the sought probability will be 
found if this sum is divided by the total number 
Q of trials, and the limit is found for Q> œ. 
The limit of each term of the sum divided by Q 
equals the probability of obtaining one of the 
numerals of interest, the sought probability 
indeed equals the sum of the probabilities of 
cach of them. Therefore, the probability of 
obtaining either 3 or 4 is 1/6 + 1/6 = 1/3. The 
probability of obtaining either 1 or 2 or 3 or4 
is 1/6 -+ 4/6 + 1/6 + 1/6 = 2/3. And the prob- 
ability of getting either 1 or 2 or 3 or 4 or 5 or 6 
is 1/6 + 1/6 + 4/6 + 1/6 + 1/6 + 1/6 = 1. This 
result is a particular case of the second 
property. 

2. Let us refer to the set of incompatible events 
that covers all possible outcomes of a given trial 
as to the complete system of events. For instance, 
in the experiment with a cube the complete 
system of events is made up of the events con- 
sisting in drawing the numerals 1, 2, 3, 4, 5, 6. 
The second property states: 

The sum of the probabilities of the events forming 
the complete system equals unity. By virtue of the 
first property, this sum equals the probability 
for one of the events forming the complete sys- 
tem to occur. But by the definition of the com- 
plete system, one of the events necessarily does 
take place. (In the example with a cube this 
means that one numeral out of the six possible 
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numerals has to turn up.) An event that cannot 
but occur is a certain event, and its probability 
equals unity. This proves the second property. 
(In the case of a cube it states that the sum of the 
probabilities for the six possible numerals equals 
unity.) 

As far as the probabilities of different values 
that can be assumed by a random variable are 
concerned, this property was formulated as for- 
mula (7) of Chapter 4. 

In some cases the reservation citing the incom- 
patibility of events may prove quite substantial 
when the rule for the addition of probabilities 
is applied. Consider the following example. 

Example. Five marksmen fire simultaneously 
at one target. The marksmen have identical 
skills: each hits the target with the probability 
of 1/3. What is the probability of at least one 
marksman hitting the target? 

We must find the probability of at least one 
of the following five events (no matter which) 
to take place: the target is hit by the first marks- 
man, it is hit by the second one, and so forth. 
A luring thought is to resort to the rule for the 
addition of probabilities. This rule states that 
the probability of one of the marksmen to hit 
the target equals the sum of the probabilities: 


P = 4/3 + 1/3 + 1/3 + 1/3 + 1/3 = 5/3 


The result is blatantly absurd. The sum is greater 
than unity, and this is nonsense. Where did we 
make a mistake? Recall that the rule for addition 
is formulated only for incompatible events. Is 
il excluded that several marksmen hit the target 
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simultaneously? Of course, they can. This is 
a typical example of compatible events. Hence, 
the rule for addition is inadmissible here. 

To solve the marksmen problem, we have 
bo make use of the rule for the multiplication 
of probabilities formulated in the next section. 


Multiplication of Probabilities 


We again throw a cube. The question that has 
to be answered is the following. We throw the 
cube twice and obtain two numbers. What is 
the probability for these numbers to be 6 and 4, 
and precisely in this order: first 6 and then 4? 

The approach to the solution is standard. We 
make Q trials (each consisting of two throws) 
and determine the number of trials giving the 
desired result. First we select all those trials 
in which the first throw gave 6 regardless of the 
result of the second trial. This problem is already 
familiar. All the faces of the cube being identical, 
the number 6 (or any other number from 14 to 6) 
appears on the first throw in 1/6 of the trials, that 
is, at the first stage we select Q, = Q/6 trials. 
(The number Q is assumed to be quite large, so 
that random deviations from Q, are small.) 
Now we want to select those trials in which the 
second throw gave 4. Each number appears on 
the second throw again with equal probability. 
Consequently, the number 4 was found in 4/6 
of the trials. Hence, the number of trials in 
which 4 followed 6 is Q, = 1/6-1/6-Q, and the 
probability of this event is 


Q/Q = 1/6-1/6 = 1/36 
8—0240 
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Now we introduce complications. Let a trial 
consist of three throws, and let us find the prob- 
ability of the trial to yield three numbers in 
a predetermined order, for instance, 4, 5, 1 
or 6, 6, 6. Arguing in a similar manner, we shall 
find that the number of trials giving the sought 
result is Q = 1/6-Q, = 1/6-4/6-1/6-Q, and the 
probability of this result is 


Q/Q = 1/6-4/6-1/6 = 1/216 


Let us analyze another example. Assume that 
one bicycle per every ten thousand bicycles 
manufactured by a bicycle plant has an unsound 
front-wheel axle, and two have unsound rear- 
wheel axles. Therefore, the probability for a 
bicycle taken at random to have an unsound 
front-wheel axle is 1/40 000, and that for the 
rear-wheel axle is 2/10 000. Assume now that 
the front- and rear-wheel axles are manufactured 
in different workshops, so that a defect in one 
of them does not increase or diminish the prob- 
ability. of a defect in the other. We want to 
find the probability for a bicycle picked at ran- 
dom to have two unsound axles. We have to 
argue as in the preceding cases. Among Q bicycles, 
select those with unsound front-wheel axle. Their 
number is Q/10 000. Among these select the 
bicycle with unsound rear-wheel axle. This gives 
(Q/10 000). (2/10 000). The sought probability is 
(1/10 000). (2/10 000) = 2-140-8, 

In both examples we fixed the probabilities 
of several events and wanted to find the proba- 
bility for these events to occur jointly, that is, 
in the same trial. The results obtained can be 
formulated in a general statement: 
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The probability for several events to occur jointly 
is equal to the product of the probabilities of these 
events. 

This rule must be supplemented with an impor- 
tantclarification. All the above examples operated 
with independent events. Two events are said 
to be independent if the probability for one of 
the events to occur is not affected by whether another 
event has occurred or not. For instance, the fact 
that the first throw of the cube gave 6 does not 
change in any way the probability to get 4 on 
the second throw, just as a defect in the front- 
wheel axle does not change the probability of 
a defect in the rear-wheel axle. 

It is not difficult to understand that the inde- 
pendence of events is essential for deriving the 
multiplication rule for probabilities. 

Let us consider again the example with bicycles 
and assume that the independence of events is 
distorted in the following manner. The front- 
and rear-wheel axles of each bicycle are assembled 
simultaneously, with the probability of defective 
axles being higher on certain days than on others. 
Then the presence of a defect in one of the axles 
increases the probability of the second axle 
being unsound, because it increases the probabil- 
ity of the bicycle as a whole to be assembled on 
an unlucky date. Hence, the probability of 
both axles being unsound increases. 

In order to better understand this point, con- 
sider an extreme situation: assume that all 
defective axles are assembled on certain days. 
All bicycles manufactured on these days have 
defective rear-wheel axles, and half of them 
additionally have defective front-wheel axles. 


3e 
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Then the probability for a bicycle selected 
randomly from the total annual output of the 
plant to have two defective axles equals the 
probability of the front-wheel axle to be defec- 
live, that is, equals 1/10 000 and not 2-1078, The 
multiplication rule for probabilities is thus 
valid only for independent events. 

The multiplication rule yields a straightforward 
solution of the problem of five marksmen for- 
mulated in the preceding section. Let us recall 
it: five marksmen shoot simultaneously at a 
target, and the probability for each of them to 
hit the target is 1/3. What is probability of at 
least one marksman hilting the target? 

The solution is most easily obtained if we 
calculate the probability for all five marksmen 
missing (let us denote this probability by Po) 
Since the missing or hitting the targel by each 
marksman must be regarded as independent 
events, the probability Py equals the product of 
the probabilities for each marksman to miss. 
The event consisting in a marksman hitting 
the target and the event consisting in this marks- 
man missing it form a complete system of events. 
The sum of the probabilities of these two events 
is unily. If the probability for a marksman to 
hit the target equals 1/3, the probability for 
him to miss it is 1 — 1/3 = 2/3. The probability 
for all five marksmen to miss is 


Py = 2/3-2/3-2/3-2/3-2/3 = (2/3) 


The event consisting in all five marksmen mis- 
sing and the event consisting in at least one of 
them hitting the target form a complete system 
of events. The sum of the probabilities of these 
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events equals unity. Therefore, the sought 
probability P satisfies the equation 


P+P,=1 
whence 
P=1—P, = 1 — (2/3) ~ 0.87 


Exercises 


1. g% sites of a network consisting of f° sites 
are blocked. What is the probability of a ran- 
domly selected site to be blocked? nonblocked? 

2. Find the probability for three consecutive 
throws of a cube to give (a) 1, 2, 3 in any order; 
(b) 4, 2, 2 in an arbitrary order. 

3. Workshop A manufactures high-quality 
parts at a probability 0.8, and workshop 2 
manufactures them at a probability 0.9. Three 
parts made by workshop A and four parts made 
by workshop B were taken at random. Find tho 


probability for all seven parts to be good, 


Percolation Threshold in a 22 Network 


The excerpts from probability theory given above 
are quile adequate for an analysis of a pereola 
tion problem in a square network consisting of 
four sites (J = 4). i 
Figure 2 gives the schematic of the experiment 
with a 2 X 2 network. The figures numbering 
the four sites are wrillen on separate slips of 
paper, the slips are pul ina hat, and the contents 
of the hat are well shaken, Assume now that the 
first extracted piece of paper has a numeral /, 


38 Part I. Site Percolation Problom 


and the site Z is blocked (Fig. 26). (The line of 
argument and the final results do not change 
in the least if the first to be blocked is a site 
with a different number. The point is that all 
sites in a four-sited network occupy equivalent 


Fig. 2. Calculations for a 2 X 2 network: 

(a) an initial network; (b) one network site is blocked; 
(c), (d) and (e) two sites are blocked. In case (c), current is 
interrupted only after a third site has been blocked, so 
that x, = 4/4. In cases (d) and (e), current is interrupted 
after a second site has been blocked, so that Z = 1/2. 
The three cases (c), (d) and (e) are equally probable. 


positions.) If the next site to be blocked is site 2, 
the current will not be cut (Fig. 2c): it will 
flow through the lower wire. After the third 
silte is blocked (2 or 4), the current is certainly 
interrupted, and one should state that the crit- 
ical fraction of nonblocked sites is 1/4. If, how- 
ever, the second to be blocked is site 3 or 4, the 
current is interrupted, and the critical fraction 
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is found to be 41/2 (Fig. 2d, e). The percolation 
threshold x, is therefore a discrete random vari- 
able assuming the values 1/4 and 1/2. Let us cal- 
culate the probability P for this variable to as- 
sume each of these values: P (1/4) and P (4/2). 

The decisive factor is which site will be blocked 
after the first. If it is site 2, then ze = 1/4, and 
if it is site 3 or 4, then x, = 1/2. Therefore, the 
probability P (1/4) equals the probability of 
site 2 being the second, and P (1/2) equals the 
probability of site 3 or 4 being the second. After 
site Z has been blocked, all three remaining sites 
have equal probabilities of being blocked at the 
next step. The sum of the three probabilities 
equals unity because these three events form a 
complete system. Hence, each of these proba- 
bilities equals 1/3. 

The probability for site 2 to be blocked next 
is thus 1/3. But if site 2 is the next, then ze = 1/4. 
Hence, the probability of £e = 4/4 is 1/3, that 
is, P (1/4) = 1/3. Now we want to find the prob- 
ability for either site 3 or 4 to be blocked at 
the second step. By virtue of the addition rule 
for probabilities, it equals the sum of probabil- 
ities: 41/3 + 1/3 = 2/3. This is precisely the 
probability for a, to take on the value 4/2; 
therefore, P (1/2) = 2/3. We must have 


P (1/2) + P (1/4) =4 
because only two values of £e are possible. Indeed, 
P (1/2) + P (1/4) = 28 +18 = 14 


The mean value of the percolation threshold 
Ze (4) is readily predictable. According to for- 
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mula (5) of Chapter 4, 


ae (4) = 1/2-P (1/2) + 1/4-P (1/4) 
= 1/2-2/3 + 1/4-4/3 = 5/12 


This figure differs quite significantly from the 
threshold value xp = lim ze (4^) thal was found 
—>co 
to be, as we have mentioned already, roughly 

0.59. 


There are no difficulties in calculating the 
variance of percolation threshold. According 
to formula (6) of Chapter 4, 


ô (4) = (4/2 — 5/12)?-2/3 + (1/4 — 5/12)? - 4/3 
= 1/72 


The root-mean-square deviation is 


6(4)=V 2/12 


Exercise 


4. Retrace the argument, assuming that the 
frst sile to be blocked was site 3. 


Continuous Random Variables 


So far we have been discussing only discrete ran- 
dom variables. However, there are also continuous 
random variables Wat can assume any value with- 
in a certain interval of the numerical avis. 
Assume that a random variable a is allowed 
to take on any value y in the range from A to B 
(A <y < B), but some of these values occur 
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frequently, while others only rarely. In order 
to describe this in terms of mathematics, the 
distribution function f (y) of a random variable a 
is introduced. 

The main property of distribution function is 
as follows: if the points A, and B, fall within the 
interval (A, B) and A; < Bı, the probability 
for the value of the random variable to be within 
the interval Ay S y < B, is equal to the area 


Fig. 3. 


encompassed by the graph of the function f (y), the 
abscissa axis, and the perpendiculars erected at 
points A, and B, (this area is shaded in Fig. 3). 
The reader familiar with integral calculus will 
recognize that this probability (let us denote 
it by P (A,, B,)) is given by the formula 
Bı 
P(A, B) = | fay 


Ay 


Since all the values of the random variable 
are within the interval (A, B) and one of these 
values is necessarily assumed by the variable, 
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the sought area equals unity. In other words, 


B 


P(A, B)= | fyay=1 (1) 
A 


Sometimes this equality is called the condition 
of normalization of distribution function. 

A figure whose area is given by the integral (1) 
is called a curvilinear trapezoid (see Fig. 3). 
If the interval (A,, B,) is so small that the dis- 
tribution function within this interval remains 
practically unchanged, the curvilinear trapezoid 
can be successfully replaced by a rectangle with 
height f (y,), where y, is an arbitrary point within 
the interval (4,, B,). Then 


P (Ay, By) = f(y) A (2) 


where A = B, — A, is the width of the interval 
(Ai, B,). 

In the mathematical literature the function 
Í (y) is referred to as the probability density. As 
we see from formula (2), if the interval is suffi- 
ciently narrow (otherwise this formula is not 
valid at all!), the probability for the random 
variable to fall within the interval is directly 
proportional to the interval width. The function 
f(y) is the probability divided by the interval 
width, or it is the probability per unit length 
of the interval, or in other words, it is the prob- 
ability density. Nevertheless, physicists often 
prefer the term “distribution function”. 

Formulas (5) and (6) of Chapter 4 for the mean 
value and variance are rewritten for a continuous 


Ch. 2. Basic Rules for Calculating Probabilities 43 


random variable in the form 


B 
a= j yf (y) dy (3) 
A 
R — 
& = | w— a (o) dy (4) 
A 


where a is the mean value of the continuous ran- 
dom variable a. 

Let us consider an example of distribution 
function. 

Uniform distribution. A continuous random 
Variable assumes all values from zero to unity 


f 

2 

1 

0 l 7 
Fig. 4. 


with equal probabilities and cannot assume any 
other value. Obviously, within the interval (0, 1) 
the function f (y) is independent of y and vanishes 
outside this interval (Fig. 4). Its value within 
the interval is easily obtained from the condition 
of normalization (1). In this case A 0,2 iN 
and the curvilinear trapezoid is transformed into 
a rectangle with area fy-1, where fọ is the value 
of the function within the interval, and the 
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i i rmalization c ition 
interval width is 1. The normalization condi 
dictates that fo-1 = 1, that is, fg = 1. Therefore 
1 at O<y<1 6) 
0 at y>1 and y<0 


’ 


iw] 


Exercise 


5. A continuous random variable a assumes 
with equal probability all values from —1 to +-4. 


Find the probability for it to fall within the 
interval from —3/4 to —1/4. 


Percolation Threshold 
as a Continuous Random Variable 


Strictly speaking, percolation threshold is a dis- 
crete random variable because all the values that 
it can assume convert to integers after being 
multiplied by the total number f° of sites. 
However, if f is very large, the difference be- 
tween the nearest allowed values of this 
random variable is very small (it equals 
ot) Therefore, in this most important 
case of a very large number of Sites, the percola- 
tion threshold z, can be considered with high 
accuracy as a continuous random variable that 
lakes on all possible values within a certain 
interval of the numerical axis. Then the quantity 
% must be characterized by the distribution 
function f (y). In this section we describe the 
shape of f (y) at large values of Ml 

The distribution function for threshold values 
Ye must depend on the number of network sites, 
«$°, with which we experiment. Therefore, it is 
more correct to denote the distribution function 
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by f p(y). It will be more convenient to define y 


nol as the threshold value as such, but as its 
deviation from the mean value ze (.f’). Then 
Tir) A is the probability for the threshold 
value found in a specific experiment to deviate 


Sh 


Fig. 5. Functions f ,.(y). The number of sites, 4^, increases 
with increasing curve number. The quantity A wv is half- 


width of curve 3 at half-height marked by the dashed 
line, 


from the mean value ze (S) by a quantity that 
falls within a small interval A in the neigh- 
borhood of y. By definition, the mean value 
calculated by using the function fil) and for- 
mula (3) equals zero. 

Figure 5 represents the function f gp(y) for three 
different values of f°. As we see from the graph, 
the distribution function becomes sharper as the 
number of sites, °, increases. This means that as 
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MWS grows, the deviations from the mean value 
(we remind the reader that it was assumed equal 
to zero!) become gradually less probable. Accord- 
ing to formula (1) of the preceding section, the 
areas under all bell-shaped curves must be iden- 
tical. With increasing i the maximum height 
of the curves increases, but the width diminishes. 
The width of a bell-shaped curve can be defined 
as the distance between the points at which this 
curve intersects the horizontal line drawn at a 
distance from the abscissa axis equal to half the 
maximum height of the curve (see Fig. 5). We 
denote this width, usually called half-width, 
by Ay 

The probability for the values of percolation 
threshold to fall beyond the half-width of the 
curve is lower than that of the most probable 
value of the threshold by a factor not less than 2. 
Hence, the half-width characterizes a typical 
spread of percolation thresholds, including the 
deviations whose probability is half that at the 
maximum of the curve ly U) 

Remember that the root-mean-square deviation 
carries essentially the same information (Chap- 
ter 1). It does not define a deflection whose prob- 
ability is exactly half as large as the maximum 
probability, but at the same time characterizes 
the typical spread of values of percolation 
thresholds. 

The quantities Aw and ôy are proportional to 
each other for any bell-shaped curve, but the 
proportionality factor depends on the shape of 
the curve. Computer calculations demonstrate 
that the percolation threshold distribution func- 
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tion is Gaussian (named after the great mathe- 
matician Karl Friedrich Gauss). This function 
has the following form: 


pig (i) 0 


where exp a =e", and e æ 2.72 is the base of 
the natural logarithm. It is plotted in Fig. 5 for 
different values of ôr- As the function is sym- 
metric with respect to the point y = 0 at which 
it reaches the maximum, the half-width Ay can 


be found from the relation (see Fig. 5) 


A 
x $ 
ty (£)= zir 0 
Making use of formula (6), we obtain 
Apy =2(21n 2) y 


By virtue of formula (8) of Chapter 1, the quanti- 
ty bp vanishes as a power function when JS > 
> oo. This means that as the number of sites 
grows infinitely, the half-width of percolation 
threshold distribution function tends to zero, 
that is, the function itself degencrates to a sharp 
peak. All values of the percolation threshold, 
except one, have zero probability. In this con- 
nection, we repeat the most important state- 
ment of the preceding chapter: as 4’ — œ, per- 
colation threshold converts from a random vari- 
able into a certain quantity. 
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Exercise 

6. (For readers familiar with integral calculus!) 
Substitute the function f(u) defined by for- 
mula (6) into formulas (3) and (4) and prove 
that the mean value @ calculated by means of 
this function is zero, and the variance ô? equals 


Bt, 


Chapter 3 
Infinite Cluster 


This chapter also deals with the site problem of 
percolation theory, but this time we formulate 
it in a different langu the language of clus- 
ters. Furthermore, we discuss a different object: 
instead of a network with blocked sites we dis- 
cuss a doped ferromagnetic, i.e. that with im- 
purity atoms. This is a much more complicated 
object, and so it must be described at least briefly. 


Permanent Magnet 


The reason why iron, nickel, cobalt, and some 
other materials can form permanent magnets 
is probably known to almost everyone. The 
explanation of this phenomenon is that the 
atoms of which such materials are composed are 
themselves elementary magnets. They possess 
magnetic moments. 

The magnetic needle of the compass is a well- 
known system possessing a magnetic moment. 
A magnetic moment is a vector. The needle of 
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the compass has southern and northern poles, 
and its magnetic moment is directed from the 
southern to the northern pole. The external mag- 
netic field makes the compass needle to turn 
so that it is oriented along the magnetic lines 
of force. Of course, any magnetic moment rotates 
in an external field in the same way. The com- 
pass needle produces an external magnetic field. 


w=(1/c)JS 


Fig. 6. Current loop and its magnetic moment. 


Any magnetic moment produces a completely 
similar magnetic field. 

It was found as early as the beginning of the 
49th century that the source of magnetism is the 
motion of electric charges, that is, the electric 
current. The magnetic moment is produced by 
the current. The magnetic moment u of a planar 
current loop shown in Fig. 6 is given by the for- 
mula 


u = (1/c) IS 


where J is the electric current, S is the loop area, 
and c is the speed of light (in the CGS system of 
units). The orientation of the vector is perpen- 
dicular to the plane of the loop and such that 
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the current flows counterclockwise if we look 
in the direction pointed by the vector’s arrow. 

If the system consists of more than one current 
loop, we can use the heat-to-tail method for the 
addition of the magnetic moments of the loops 
and find the total magnetic moment of the system. 

What is the origin of the alomic magnetic 
moment? Any atom is known to consist of a heavy 
nucleus and an electron shell. The magnetism 
of solids originates precisely from the magnetic 
moment of the shell (the atomic nucleus also 
can have a magnetic moment, but it is approx- 
imately a thousandth of that of the shell). 

First, the magnetic moment of the shell grows 
from the motion of electrons around a heavy 
nucleus. This rotation can be pul in correspond- 
ence with a certain current J and an effective 
area S. Besides, quantum mechanics ascribes 
to each electron an additional magnetic moment 
called the spin moment. This last moment is in 
no way related to the characteristics of motion 
of the electron, but represents its inherent prop- 
erly. However, the spin moment produces 
a magnetic field just as ordinary moment does. 
Most often the net magnetic moment of the 
electron shells of the atoms of which a solid 
consists is zero. However, in some materials 
such as iron, nickel, cobalt, and some others 
the electron shells possess a nonzero magnetic 
moment. 

The magnetic moments of neighboring atoms 
in a solid interact with one another. In principle, 
this interaction is similar to the interaction 
between two compass needles placed close togeth- 
er. Each needle produces a magnetic field acting 
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on the other needle. However, the situation 
becomes significantly more complicated because 
the interaction takes place not in the vacuum. 
The outer electron shells of atoms decisively 
affect the character of the interaction, up to 
reversing the directions of the applied forces. 
h Isperiments show that in some materials the 
interaction between magnetic moments is such 


Fig. 7. Fragment of the crystal lattice of a ferromagnetic. 
The arrows show the orientations of magnetic moments. 


that the forces between them make the moments 
align in the same direction. These materials 
are called ferromagnetics (Fig. 7). 

If the magnetic moments of all atoms are 
oriented in the same direction, the total magnetic 
momente# equals the arithmetic sum of individual 
moments: of = wf”, where M is the number of 
atoms in the solid, and p is the magnetic moment 
per atom. 

As the size of a body increases, its magnetic 
moment grows proportionally to the volume of 


4e 
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the body (the number of atoms, .f, is propor- 
tional to volume). A specific characteristic of 
magnetic properties, that is, the quantity inde- 
pendent of the size and dependent only on the 
properties of the atoms composing the body, is 
the spontaneous magnetization M. It is defined 
as the magnetic moment of a unit volume, that 
is, it is found as the total moment o% divided 
by the body volume V: 


M i ey es i 
¥ Vv vo 
where vy = V/S is the volume per atom. 

The word “spontaneous” emphasizes that the 
magnetization M is not induced by an external 
magnetic field but appears because of internal 
forces. A permanent magnet is just such a body 
in which the spontaneous magnetization is non- 
zero. ‘This magnetization produces a magnetic 
field in the medium surrounding the magnet 
(or in the vacuum). 

The unit of spontaneous magnetization in the 
CGS system of units is one gauss (1 G). For in- 
stance, in iron kept at very low temperatures 
M = 1740 G. From this figure we can find the 
magnetic moment p per atom. It is approximately 
2.2 of the spin magnetic moment of the electron. 
The fact that the moment p is found to be of the 
order of the spin moment confirms the correctness 
of our concept of the nature of spontaneous mag- 
netization. 

Thermal motion destroys the magnetic ordering, 
so that there is a critical temperature called 
Curie temperature (or point), above which the 
spontaneous magnetization is zero. For instance, 
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the Curie temperature of iron is 770°C. Iron 
cannot form permanent magnets at higher tem- 
peratures. 


Doped Ferromagnetics 


Now let us consider a material which is a solid 
solution (mixture) of magnetic and nonmagnetic 
(i.e. having a zero magnetic moment) atoms. 
This is a crystal in whose lattice sites magnetic 
and nonmagnetic atoms sit, their arrangement 
being quite random. 

Let us assume that the interaction between 
the magnetic moments of the atoms decreases 
with distance so fast that we have to take into 
account only the interaction between nearest 
neighbors. This means that if two magnetic 
atoms are at neighboring lattice sites, their 
moments are necessarily parallel, but if they are 
separated by at least one nonmagnetic atom, 
their moments’ orientations are arbitrary: they 
“know nothing” about each other. 

The question that we want to pose now is 
whether the spontaneous magnetization arises 
in the presence of nonmagnetic atoms, and how 
many nonmagnetic atoms are needed to destroy 
the spontaneous magnetization. It will be shown 
below that the answer to this question reduces 
to solving the site percolation problem that 
has been formulated in Chapter 4. 

Let us make some definitions. Two magnetic 
moments will be said to be connected with each 
other if they are neighbors or if they are connected 
via a string of neighboring magnetic atoms (Fig. 8). 
The phrase “are neighbors” signifies that the 
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atoms are nearest neighbors. In a square lattice 
shown in Fig. 8 the nearest neighbors are the 
neighbors along the horizontal and vertical 
directions, but not the neighbors along the dia- 
gonals. An ensemble of connected aloms is said 
to form a cluster. This definition carries the 


Fig. 8. Fragment of a plane lattice with magnetic (open 
circles) and nonmagnetic (black circles) atoms. The mag- 
netic atoms form one four-atom cluster, one two-atom 
cluster, and five one-atom clusters. The dashed lines de- 
lineate the cluster boundaries. The moments of different 
clusters can have different orientations. 


following meaning. The magnetic interaction 
makes the connected atoms orient their magnetic 
moments in the same direction. As a result, 
each cluster has a resultant magnetic moment 
proportional to the number of atoms of which 
it consists. Moreover, we agreed that magnetic 
atoms which are not nearest neighbors do not 
interact at all. Consequently, the atoms belong- 
ing to different clusters do not interact with 
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one another, and therefore, the relative orienta- 
tion of the magnetic moments of different clusters 
is arbitrary (see Fig. 8). 

Let us denote by x the fraction of magnetic 
atoms, that is, the ratio of the number of mag- 
netic atoms to the total number of lattice sites. 
By definition, x varies within the interval from 
0 to 4. 

First we assume that the fraction of magnetic 
aloms is very small (z <1). Obviously, in this 
case they are primarily separated (as raisins in 
a bun). A cluster made of two magnetic atoms 
is a rare event, a cluster of three atoms is still 
rarer, and so on. This statement is extremely 
important for what follows and will be proved 
mathematically a little farther below. So far 
we suggest that those who refuse to recognize 
this statement as obvious should take it for 
granted. 

At « <1 the number of clusters is therefore 
approximately equal to the number f° of mag- 
netic atoms, and hence, grows proportionally 
to, as the total number of lattice sites increases. 
ITowever, the magnetic moments of these clusters 
“know nothing” about one another, and hence, 
are randomly oriented with respect to each 
other (see Fig. 8). In order to find the net mag- 
netic moment o% of the system, we need to add 
up the magnetic moments of individual atoms 
by the head-to-tail method. By virtue of random 
directions, these moments cancel out, so that 
the spontaneous magnetization is ultimately 
zero. We have thus found that the spontaneous 
magnetization vanishes at low concentrations of 
magnetic atoms. 
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Formation of an Infinite Cluster 


Now consider the case in which almost all atoms 
are magnetic. Obviously, a small admixture of 
nonmagnetic atoms does not cancel the spontane- 
ous magnetization but only diminishes it. Let 
us discuss this effect in terms of clusters. At 
x = 1 all £ atoms belong to a single cluster. 


Fig. 9. Fragment of a plane lattice with magnetic (open 
circles) and nonmagnetic (black circles) atoms at a high 
concentration of magnetic atoms. With the exception of 
atom B, all magnetic atoms belong to a single cluster and 
have identically oriented magnetic moments, 


If x only slightly differs from unity, some atoms 
drop out of this cluster. This occurs, first, because 
some atoms are replaced by nonmagnetic atoms 
(atoms A in Fig. 9), and second, because some 
magnetic atoms form isolated clusters (atom B 
in Fig. 9) with their own orientation of magnetic 
moment. Nevertheless, when z is close to unity, 
a single cluster survives, permeating the whole 
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lattice however large it may be. This cluster is 
said to be infinite. 

Of course, this concept acquires a rigorous 
meaning only for an infinite system. Let us take 
a large series of samples with a given number of 
magnetic atoms, each containing the same total 
number of atoms, and in each sample let us 
find a cluster with the maximum number of mag- 
netic atoms. Now we average the number of 
magnetic atoms belonging to the maximum cluster 
over all samples of the series and denote the 
result of averaging by SM max- Therefore, M max 
is the average number of atoms in the largest 
cluster. The quantity WS max is a function of 
both g and z. The existence of an infinite cluster 
has the following corollary: at a given value of x 
the ratio SM max/ 4 tends, when ° increases 
infinitely, to a nonzero limit 


lim mas — P (2) 
San I 


The fraction of atoms P (x) belonging to the larg- 
est cluster does not depend on the number 4° 
of atoms if f is sufficiently large but depends 
on z. And as f° goes to infinity, the quantity 
MS max also tends to infinity. It is for this reason 
that we speak of an infinite cluster. 

Only one infinite cluster can exist in a system. 
Assume that not only the average number of 
atoms in the largest cluster is prescribed at 
the given values of .f and x but also the average 
number of atoms in the next-largest cluster. 
Let us denote this last number by Smax. By 
definition, WS frax < Smax- The statement that 
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only one infinite cluster exists in a system signi- 
fies that 


lim 
MW +00 
at any value of z. This means that two clusters 
permeating the system necessarily merge some- 
where¥and transform into a single cluster.* 

We have thus obtained that at a sufficiently 
high concentration z of magnetic atoms, a certain 


P(x) a(x) 
o 


(1) 
1 
0 


Xe OS EA 


Fig. 10. Graphs of functions P (x) and o (2)/o (4). Although 
the two functions vanish at a common point, at the 
critical value x, they behave quite differently, for reasons 
to be discussed in Part III of the book. 


fraction of these atoms belong to one cluster, 
and they all have a common direction of atomic 
magnetic moments. This means that the spontane- 


* Strictly speaking, the statement on the survival of a 
single infinite cluster has not been proved. There are 
arguments supporting the validity of this conjecture to a 
greater or lesser extent, but it would be more correct to 
say that specialists simply take it for granted, 
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ous magnetization is nonzero: 
M=- P(z) 
vo 


Recall now that only small clusters exist at 
a low concentration x of magnetic atoms. In this 
case an increase in the number of sites, f°, results 
only in the growth of the number of small clusters, 
but not in a larger number of particles in each 
cluster. Then 

lime 3 0, that is, P (x) =0 
M -+0 


We thus have to conclude that there is a critical 
concentration x, at which an infinite cluster is 
formed, x satisfying the inequality 0 < te < 1. 
At this very concentration z, the spontaneous 
magnetization appears, and the function P (x) 
becomes distinct from zero (Fig. 10). Consequent- 
ly, a material cannot form a permanent magnet if 
the fraction of nonmagnetic atoms is greater than 
1 — xe (the fraction of magnetic atoms is less 
than ag): 


Exercise 


1. Find the function P (x) for x not very differ- 
ent from unity. 


Sile Percolation Problem Revisited 


What remains to be done now is to say that from 
the point of view of the critical concentration x, 
the problem of network conduction and the prob- 
lem of doped ferromagnetic are identical. 
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The electric conduction problem can be refor- 
mulated just as easily in cluster terms. We then 
only want to replace the concept “nonmagnetic 
atom” in all definitions by the term “blocked site”. 

Figure 8 shows a configuration of magnetic 
(open circles) and nonmagnetic (black circles) 


Fig. 11. The same configuration as in Fig. 8, but with 
nonblocked sites substituted for magnetic atoms, 


atoms. Let us perform the above-described repla- 
cement in this configuration and convert it 
from a doped ferromagnetic to a wire mesh with 
cut-out knots. This is done by removing in 
Fig. 8 the arrows indicating the directions of 
magnetic moments and then by tracing the wires 
so that they connect the lattice sites (Fig. 11). 

Figure 11 clearly shows the main property of 
clusters applied to the network problem. if 
a potential difference is applied to any pair of 
sites within one cluster, a closed circuit with electric 
current is formed. (Of course, this property is 
meaningful only for clusters containing at least 
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two siles.) When a potential difference is applied 
to any pair of sites belonging to different clusters, 
a closed circuit is not formed and electric current 
is zero. If £ < £e, the system contains only clus- 
ters composed of a finite number of sites, so 
that as the size of the system increases, the cur- 
rent between the lateral electrodes will neces- 
sarily be interrupted sooner or later. But if 
> ae, a very large system cannot but include 
into its lateral faces some sites belonging to an 
infinite cluster. This infinite cluster will create 
the electric conductivity o (x) distinct from zero 
and independent of the size of the system. 

Let us return to Fig. 10 which shows the func- 
tions P (x) (the fraction of the sites belonging 
to an infinite cluster) and o (x)/o (1) (6 (1) is 
the electric conductivity at x = 1, i.e. with 
no blocked sites). Both functions vanish at the 
same point that we first identified as the per- 
colation threshold and later the point at which 
an infinite cluster is born. 

We were thus always dealing with a problem 
of percolation theory that is referred to as the 
site percolation problem. If we were interested 
in the value of £e for a “plane ferromagnetic”, 
then we could look up the result of the experiment 
with the wire mesh and say that x, = 0.59. 
However, actual ferromagnetic materials crystal- 
lize into three-dimensional, not plane, lattices. 
An example of a three-dimensional lattice is the 
primitive cubic lattice whose unit cell has already 
been shown in Fig. 7. 

í The problem of wire mesh electric conductance 
is easily generalized to the three-dimensional 
case. Imagine a cube consisting of numerous 
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cells and made by welding together pieces of wire, 
as shown in Fig. 12. We can solder two metal 
plates on the opposite faces of this cube and, as 
in Fig. 1, make an electric circuil and study its 
electric conductance as a function of the number 
of blocked sites. The blocking of each site in ter- 
ires entering 
two-dimensional case, there 


rupts the contact between six w 
this site. As in the 


is a critical Concentration zg of nonblocked sites 
below which the electric conductance vanishes, 

The problem of doped ferromagnetic and the 
related concept of an infinite cluster were equally 
valid for plane and three-dimensional lattices. 
The critical Concentration a, of Magnetic atoms 
at which an infinite cluster appears is at the same 
time the threshold of percolation from one cube 
face to the opposite face in a sufficiently large 
cube. It must be borne in mind that the quantity 
Ta is strongly dependent on the type of the lat- 
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lice. It was found to be 0.59 for a square lattice, 
but ze = 0.31 for a primitive cubic lattice. 


(For details see Chapter 6.) 


Clusters at a Low Concentration 
of Magnetic Atoms ** 


The conclusions drawn in the preceding sections 
were based to a large extent on the statement that 
at a low concentration x of magnetic atoms these 
atoms are usually single, two-atom clusters are 
rare, three-atom clusters are even rarer, and so on. 
Let us prove this statement. 

We introduce a function Pl) giving the 
probability for a randomly chosen atom to belong 
to a cluster consisting of not less than o# atoms. 
This means that the atom chosen at random is 
(a) magnetic, and (b) connected to not less than 
oll — 1 other magnetic atoms. Let us calculate 
the function P g(£) for æ = 4 and ch = 2. 

The function P, (æ) is the probability for a 
randomly chosen atom to be magnetic. This 
probability equals æ (see Exercise 1 to Chapter 2 
where the word “nonblocked” must be replaced 
by “magnetic”, and the word “blocked” by “non- 
magnetic”): 

P, (&) = 2 (1) 

The function P, (x) equals the probability for 
a randomly chosen atom to be magnetic and to 
have another magnetic atom among its nearest 
neighbors. Obviously, these two events are inde- 
pendent, so that the sought probability can be 
represented by the product of the probabilities 
of these two events. The first of them (the prob- 


——— 
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ability for the atom to be m 
and thus 

P, (x) = <W (a) (2) 
where W (x) is the probability for at least one 
magnetic atom to be found among the nearest 
neighbors of the atom. The function W (z) 
depends on what kind of lattice we consider. 
Let us limit the analysis to the square lattice 


agnetic) equals z, 


© 
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Fig. 13. 


in which each atom has four nearest neighbors 
(Fig. 13). We need to find the probability for 
at least one of the atoms 7, 2, 3, 4 to be mag- 
netic. 

This problem is best solved 
event consisting in all four 


netic and the event consis 
of the four 


as follows. The 
atoms being nonmag- 
ting in at least one 
m a complete 
probabilities 


W. As we have said, W +W, 
bility for atom Z to be nonmagi 
The probability for 
magnetic also equals 


= 1. The proba- 
netic equals 1 — z. 
atom 2 or 3 or 4 to be non- 


1 — z. The events consisting 


Ch. 3. Infinite Cluster 65 


in different atoms being nonmagnetic are inde- 
pendent. Consequently, the probability for all 
four atoms to be nonmagnetic equals the product 
of all four probabilities: Wy = (1 — x)’. This 
yields W = 1 Wy = 4 (l x)*. By virtue 
of formula (2), 


P, (z) = z [1 — (1 —2)'] (3) 


If « <1, _the expression for P, (x) can be sim- 
plified by dropping the terms wilh high powers 
of a. Using the binomial formula, we obtain 
1— (1 — x)! = 4r — 62? + 42° — xt 

Note that if « <1, the ratio of each next lerm 
lo ils predecessor is small: 

Gr? 3 i 4 2 

rio w€ i, wat <t, and 


zi z 
E 


Therefore, if x < 1, we can write to a good accu- 
racy that 

1 — (1 — x) x 4x 

whence 

P, (2) = 42? (4) 
A comparison of formulas (1) and (4) shows that 
at x <1 the ratio 

P, (2) 

Sera wad 1 5 
Pay © 4K (5) 
that is, the probability that a randomly chosen 
atom belongs to a cluster of two or more atoms 
5—0240 
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is much less than the probability for it to form 
e-atom cluster. 

a Koska pler derivation of formula (4) for P, (z) 
can be given, immediately taking into account 
the condition z < 1. It is relatively simple with 
this derivation to calculate the functions Poy (x) 
for of > 2 (see Exercise 3). The derivation is as 
follows: a cluster consisting of more than two 
atoms, one of them being atom 0, necessarily 
incorporates either alom 7 or 2 or 2 or 4. The 
probability for atoms 0 and Z to belong to one 
cluster equals the probability of both of them 
being magnetic, and is given by the product 
of the probabilities for each of these atoms to 
be magnetic, that is, t-z = z?. The same can 
be said about the probability for a cluster to be 
formed by atoms 02, 03 or 04. All these pro- 
babilities equal 22. The probability for at least 
one of these events to realize equals the sum of 
the probabilities, that is, 4x°, hence we arrive 
at formula (4). 

This conclusion holds only for z <1, Only 
under this condition can we use the rule for the 
addition of probabilities, Indeed, the addition 
tule holds only for incompatible events. But the 
event consisting in atoms 0 and J being magnetic 
is compatible with the event of atoms 0 and 2 
being magnetic. The simultaneous realization 
of these events means that all three atoms 0, 7, 
and 2 are magnetic, and hence, form a three-atom 
cluster. The probability of the compound event 
equals the product of the probabilities for all 
three atoms to be magnetic, that is, equals 
Trg = 2, 

If z <1, this probability is much less than 
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the calculated probability of the formation of 
a two-atom cluster. Therefore, the probability 
for the events to occur jointly can he neglected, 
and the events can be treated as incompatible. 
This justifies the derivation given above, under 
the condition z <1. 

In fact this means that if x <1, the probability of 
the formation of a three-atom cluster can be neglected 
in calculating P, (x). 

Therefore, in the case of x <1 the function 
P, (£) actually coincides with the probability 
for a site chosen at random to belong to a cluster 
consisting of two (and nol more than two!) 
atoms. Correspondingly, the function Ps (x) 
describes a three-atom cluster. It is proportional 
to z? and is small in comparison with P, (x). 
The general result is that the function Py (x) 


contains the powers of x not smaller than 2%, 
and that for «<1 we should have P y2) < 


K Pan, (2). 

We thus found that if a site chosen at random 
in the case of xv <1 is magnetic, it almost cer- 
tainly will form a one-site cluster. The probabil- 
ity for it to belong to an o/-site cluster sharply 
drops with increasing of. 


Exercises 


2. Find P, (x) for the primitive cubic lattice 
shown in Fig. 12. Find it for an arbitrary lattice 
in which each atom has z nearest neighbors. 

3. Find P, (x) for a square lattice, making use 
of the condition z <1. 

4. Find P, (x) for a square lattice, without 
employing the condition z <1. 


5a 
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Chapter 4 


Solution 

of the Site Percolation Problem 
by Monte Carlo 

Computer Techniques 


The Monte Carlo method is one of the most wide- 
spread methods of solving Percolation theory 
problems. The aim of this chapter is to give 
a general idea of this method, to explain in detail 
how the main element of the method—the gene- 
rator of random numbers—operates, and to give, 
in conclusion, a concrete computer program that 
makes it possible to determine the percolation 


threshold of the site problem. The very first 
question that arises now is: 


Why Monte Carlo? 


“Bul what is zéro? You see th 
curly-headed one, the chief one, 
And why did he Scoop up every 
the table? Such a heap, he took 
What is the meaning of it?” 

“Zéro, Granny, means that the bank wins all. 


If the litte ball falls on zéro, everything on the 
lable goes to the bank...” 


“You don’t say so! And I shall gel nothing?” 
“No, Granny 


, if before this you had staked on 
zéro you would have got thirty-five times what 
you staked.” 


“What! Thirty-five times, and does it often 
turn up? Why don’t they stake on it, the fools?” 


al croupier, the 
showed zéro now? 
thing that was on 
it all for himself, 
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Ee a are thirty-six chances against it, Gran- 

“What nonsense. Potapitch! Potapitch! Stay, 
I’ve money with me—here.” 

She took out of her pocket a tightly packed 
purse, and picked out of it a friedrich d’or. 

“Stake it on the zéro at once.” 

“Granny, zéro has only just turned up,” I 
said, “so now it won’t turn up for a long time. 
You will lose a great deal, wait a little, anyway.” 

“Oh, nonsense; put it down!” 

“As you please, but it may not turn up again 
till the evening. You may go on staking thous- 
ands; it has happened.” 

“Oh, nonsense, nonsense. If you are afraid of 
the wolf you shouldn’t go into the forest. What? 
Have I lost? Stake again!”* 

This excerpt from Dostoevsky’s “Gambler” 
describes the most exciting game of chance of 
the last century: the roulette. It should be noted 
that from the standpoint of probability theory, 
the inexperienced cestatic grandmother reveals 
more common sense than the Gambler who is her 
consultant. The probability of getting zéro will 
not diminish in the least if it turned up on the 
preceding round. There is no sense in wailing, as 
the Gambler advises. This frequently encounter- 
ed misconception is presumably based on the 
misunderstanding of the fact that the probabil- 
ity of getting zéro lwo limes running is small. 
But this does not mean at all that if zéro turned 
up once, the probability to get it the second time 
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on the next throw is less than on the first. Of 
course, the probability remains absolutely the 
same. The city of Monte Carlo in the Monaco 
principality earned its fame as the world capital 
of roulette. It was this city that gave its name 
to one of the most powerful among the modern 
numerical methods in mathematics. 

And what is it that this method and the roulette 
share? It is the fact that the main element of the 
Monte Carlo method is that very revolving 
wheel that decides the fates of people, destroying 
some and rewarding others in numerous casinos 
of Monte Carlo. Actually, mathematicians have 
greatly improved it. There is no revolving wheel, 
but a standard computer program that is called 
the “random-number generator”. But this does not 
alter the principal point. From the mathematical 
viewpoint the wheel of the roulette game is 
nothing but a random-number generator. 


What Is the Monte Carlo Method? 


As a rule, the term “Monte G 
applied to any mathem 
ially based on a r 


arlo method” is 
atical technique essent- 
andom-number generator, 
Usually a modern computer has a stand 
program that generates random numbers 
ributed randomly within the ir 
unity, that is, it “plays oul? 
tinuous random variable th 
probability all values within the interval (0, 1). 

Each time the program is addressed, it out- 
puts one such number with a predetermined 
number of decimal places that depends on the 
computer model, 


ard 
dist- 
iterval from zero to 
the values of a con- 
at assumes with equal 
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The simplest use of the Monte Carlo method 
is for calculating integrals. Imagine, for instance, 
that we need to calculate the volume within a 
closed surface of a complicated shape. Let us 
choose a cube that we know for certain to in- 
clude the whole surface (Fig. 14). Now we obtain 


Fig. 14. To the Monte Carlo determination of a pear’s 
volume. 


from the random-number generator a set of points 
distributed uniformly within the cube. This is 
done as follows. Assume that the cube edge is 
L long, and all the three coordinates of the points 
within the cube vary from zero to L (see Fig. 14). 
Addressing three times the random-number gen- 
erator, we obtain three numbers Yi, Ya, Ya within 
the interval (0, 1). These numbers will give us 
the coordinates of the first point within the cube 
by the formulas X, = Lyn Yı = La Z, = Lys 
Repeating this procedure Q times, we obtain Q 
points that on the average uniformly fill the 
cube. Let Q, be the number of points that hap- 
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pened to fall within the surface. The points being 
distributed uniformly, the number Q, character- 
izes the volume enclosed by the surface. Name- 
ly, if Q is sufficiently large, the sought volume 
equals LQ,/Q. : ? 

A theory is available that makes it possible 
to predict the number of points, Q, required to 
obtain a result with a prescribed accuracy. An 
empirical approach can be used: we repeat the 
experiment several times, each time using a 
new set of random numbers, and then compare 
the results. If the results fluctuate within the 
prescribed accuracy, everything is fine and the 
answer is correct. The Monte Carlo technique 
has very important advantages over conventional 
integration methods in the case of multidimen- 
sional spaces (multiple integrals). 

There are cases in which the Monte Carlo 
method is the only one practicable. Imagine 
that we want to study the behavior of a system 
consisting of an enormous number of particles, 
for example, the behavior of a gas. In principle, 
this problem must he solved by the methods of 
statistical physics; however, these methods are 
found to be ineffective if the interaction between 
Particles is strong (this happens when density 
is high and temperature low). Then the proper- 
ties of the gas are studied by computer simula- 
tion. The number of gas particles participating 
in the simulation is determined by the size of 
computer memory. The memory must store the 
information on the coordinates of all particles. 
The simulation consists in randomly picking 
up one of the particles which then travelsa ran- 
dom distance. (This means that the coordinates 
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of this particle in the computer memory change.) 
Then another particle is chosen randomly, and 
so on. The potential energy of the interaction 
between the particles in the gas depends on their 
mutual arrangement. The energy is calculated 
from the very beginning and then is recalculated 
after each displacement. The probabilities for 
particles to be displaced by a specific distance 
are chosen to agree with the potential energy 
in such a way that the model system “live” on 
the average as the real system does. 

As a result, the computer memory stores what 
resembles “instantaneous photographs” of the 
gas taken at successive time moments. The 
photographs include the coordinates of all gas 
particles and thus make possible the calcula- 
tion of mean thermodynamic characteristics, such 
as pressure, heat capacity, and some others. 

The very procedure of simulation much re- 
sembles a game which is ruled by a rigid code 
that includes calling the roulette, that is, the 
random-number generator. Slight deviations from 
the rules, or a “crooked” roulette, make some 
configurations of atoms in the gas appear more 
often than others; this affects the result of aver- 
aging and invalidates the answer. 

Random-number generators are used not only 
in the Monte Carlo method but also in the so- 
called analogue experiments exemplified by the 
experiment with a wire mesh described in Chapter 
1. We have mentioned earlier that the random 
sequence of blocked sites required for this ex- 
Periment was generated on a computer. In order 
to choose the next sile, it is necessary to address 
the program and obtain a random number y. 
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e must multiply it by the total number „f° of 
ve and then add unity to the product. Then 
we have to take the integral part of f° y F 4. 
This gives an integer that falls in the prescribed 
interval from 1 to. f*. In fact, such numbers may 
appear more than once. But there is no harm in 
it. If it is found that the site with the number 
just generated was blocked earlier, it is neces- 
sary to address the computer for a new random 


number and convert it into the number of the 
site. 


Somewhat later we sh 


all describe a program 
with which we can calcul 


ate percolation threshold 
via the Monte Carlo method, but now we shall 


discuss the principal element of this method, 
namely, the random-number generator. 


How to Think Up a Random Number 


We thus need random 


numbers distributed uni- 
formly in the r 


ange from zero to unity. So far 
the problem is not yet rigorously formulated. 
We need to know how many decimal places are 
required in each number. Let us assume that 
only two decimal places are necessary. Then take 
a telephone directory, open it at an arbitrary 
page, and copy a sequence of the last two figures 
of each telephone number, preceding these numer- 
als with “0.”. This will give an acceptable two- 
digit table of random numbers. And what are 
we to do if ten-digit random numbers are required? 
Looks like a computer is then unavoidable. 

Come to think of it, the very idea that a com- 
puter can generate random numbers may appear 
strange. Indeed, any computer follows the alga- 
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rithm that was fed into it, thatis, it performs pre- 
cisely the actions programmed by a man. How 
then can a chance element be introduced into 
its performance? 

In fact, there is no chance element in the pro- 
gram of the random-number generator. The princi- 
ple of its functioning islas follows. When the pro- 
gram is addressed for the first time, some num- 
ber yọ must be fixed. Then a rigorously defined 
sequence of operations transforms this number 
into a new number 


yı = ® (yo) (1) 


where @ is a specially chosen function or a se- 
quence of operations that transform yo into yj. 
It is this function that determines the algorithm 
of random-number generation. In its turn, the 
number y, is the starting point for generating 
the next number y, by the same recipe: 


Yo = P (yı) (2) 


Obviously, the function @ is such that all 
numbers yy, Yo, «+ +> Yn Satisfy the inequalities 
O0<%m<1. This sequence of numbers repre- 
sents the sought sequence of random numbers. 

It can be readily ascertained that the thus ge- 
nerated sequence cannot be infinite. Indeed, 
any computer operates only with numbers con- 
taining a limited number of decimal places. 
Their number is limited. (There are only 10° 
two-digit numbers and 10” n-digit numbers.) 
Therefore, sooner or later a number yn will 
coincide with an earlier generated number, say, 
Yn-ı. Then the sequence will be repeated: Yn+. 
will coincide with yn-p+p and so on. 
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Consequently, the sequence generated by formulas 
(1) and (2) is inevitably periodic. In view of 
this fact, such numbers are not truly random and 
are referred to as pseudorandom (i.e. “as if ran- 
dom”, or resembling random numbers), 

Nevertheless, they can be used as random 
numbers provided the sequence of numbers re- 
quired to solve a given problem is shorter than 
the period L of the sequence. 

In its turn, the period Z is determined by the 
number of decimal places with which the com- 
puter operates (i.e. the number of memory cells 
allotted to each number) and by the quality 
of the algorithm (i.e. by the properties of the 


function ® that enters into formulas (1) and 
(2)). 


The development of a good random-number 
generator is a very difficult problem. As a rule, 


“off-the-cuff” generators prove to be of poor qual- 
ity. Some specific generators are discussed below. 


The Mid-Square Method 


Historically, this was the first ever 
random-number generation by a computer. It 
was suggested in 1946 by the brilliant mathe- 
matician John von Neumann. The method makes 
it possible to generate random numbers with 
an arbilrary number of decimal places, corro- 
sponding to the capabilities of the computer, 
The method is extremely simple. Assume that 
we need four-digit numbers. Let us choose the 
first number Ai arbitrarily. For instance, 
X, = 8219. Raise it to the second power. This 
gives an eight-digit number 67551961. Extract 


method of 
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the four middle digits: 5519. The second number 
of the sequence is X, = 5519. Now raise 5519 to 
the second power, obtaining 30459361. The 


third random number is then X, = 4593. If 
the first digits among the middle group are zeros, 
the resultant number has a smaller number of 
nonzero digits. For instance, X = 21095649, 
X, = 956. Squaring X3, we can obtain an eight- 
digit number by adding zeros on the left, 
X3 = 00913936, so that X, = 9139, and so 
forth. 

Random numbers yn distributed uniformly 
within the interval from zero to unity are ob- 
tained from the numbers X, via the formula 
Yn = X,/10', where z = 0, 1, 2, 3, ..., so that 
Yo = 0.8219, y, = 0.5519, yp = 0.4593, and so 
on. 

At first glance, the method looks attractive. 
Ilowever, a careful investigation demonstrated 
that this is not true. The main shortcoming of 
the method is that some starting numbers make 
the sequence “go into a cycle”. For instance, 
it was found that in the four-digit class of num- 
bers the sequences often terminate by a cycle 
6100, 2100, 4100, 8100, 6100. The period of this 
cycle is mere 4, and this is obviously unaccept- 
able. 

There is even a number that immediately 
reproduces itself. ‘This is 3792 (3792? = 14379264). 
Zero also reproduces itself, and quite often se- 
quences generated by the mid-square method 
degenerate to zero. Consequently, nowadays 
the mid-square method is only of historical in- 
terest. 
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Exercises 


1. Compose a sequence of four-digit numbers, 
starting with 0085, 0067, 0032. Show that all 
three sequences are monotone decreasing (each 
subsequent number is smaller than the preceding 
number) and rather rapidly degenerate to zero. 

2. Now prove that this constitutes the general 
shortcoming of the mid-square method: if 2n- 
digit numbers X; are used and a number b appear- 
ing in the sequence has zeros for n leading digits, 
the sequence becomes monotone decreasing and 
finally degenerates to zero, 


Linear Congruent Method 


At the present moment this method of generat- 
ing random numbers is r 


egarded as the best. 

The idea is as follows. Four integers are chosen: 

(1) multiplier k, 

(2) shift c, 

(8) modulus m, 

(4) the first number of the se 

The sequence of random 
mined by the formula 


Xn =(kXp +0) mod m 


quence X,. 
numbers is deter- 


(3) 
where the subscript n runs through UAR Pe REE 
The symbol b mod m denotes the remainder after 
dividing b by m. For instance, 

b 2 6 30 3 147 
m 10 10 10 49 12 


bmodm 5 6 0 3 3 


Ch. 4. Site Problem Solved hy Monte Carlo Method 79 


Obviously, b mod m < m. Consequently, all 
numbers of the sequence, Xn, satisfy the inequal- 
ity X, < m. The sequence of the numbers yn 
distributed uniformly within the interval from 
zero to unity is obtained by the formula 


m=i, n=O, 1, 2... (4) 

It is not just any choice of the four starting 
numbers that leads to good results. Note first 
of all that the sequence X, must necessarily be 
periodic, and the period cannot be greater than 
m. Indeed, all X, being integers, with Xn < m, 
the number of different numbers cannot exceed m. 
For this reason, a number that has already oc- 
curred in the sequence will appear at least be- 
ginning withn = m, and the sequence will re- 
peat itself. 

However, it is far from simple to generate a 
sequence with a maximum possible period L = m. 
If the starting numbers are not carefully selected, 
the generated sequences will have, as a rule, 
short periods. 


Exercises 


3. Write the sequence of numbers Xn generated 
by means of formula (3) for k = 3,¢ = 0, X, = 5, 
m = 20. 

4. Write the sequence of numbers Xp generated 
by means of formula (3) for k = 3,¢ = 1, Xo = 
= 5, m = 20. 

5. Write the sequence of numbers Xn generated 
by means of formula (3) fork = 3,c = 2,X = 5, 
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m = 20. Make sure that in all three cases the 
period of the sequences is essentially shorter than 
20. Analyze other examples. , 

The following theorem is true. Ifa sequence is 
generated by means of formula (3) for c #0), its 
period equals m if and only if the following condi- 
tions are met: 

(i) c and m are coprime numbers 
common divisor than 1); 

(ii) b= k— 1 isa multiple of p for any prime p 
which is a divisor of m; 

(iii) b is a multiple of 4 if m is a multiple of 4. 
Unfortunately, the proof of this theorem is too 
complicated to be given here. 

6. Make sure that the conditi 
the theorem given above did not hold in all the 
examples presented in Exercises 3-5, 

7. Make sure that the set of integers k = 11, 
c = 3, m = 5 meets the conditions of the above 
theorem and yields the period L = 5 for arbi- 
trary X,. 

Therefore, a generator with a maximum possi- 
ble period L will be obtained if we take for m the 
largest number with which a given computer can 
operate, and choose the other numbers in accord- 
ance with the theorem given above. 

However, the period is not the only character- 
istic of the quality of a random sequence. For 
instance, let us consider a sequence gencrated 
wilh k = ¢ = 1. The Sequence is 0,1,2,3, .. re 
th 4, 0;. dy. 2 2 Gee =A, Oe a... Its 
period is indeed m, but a 
is absolutely unacceptable. 

A complicated system of tests has been elabo- 
rated to determine the quality of a random- 


(have no other 


ons imposed by 
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number generator. Therefore, only reliable ge- 
nerators can be advised for applications. 

When a random-number generator is being 
chosen, the properties of the computer are im- 
portant not only from the standpoint of the choice 
of the period of maximum possible length. The 
rate at which random numbers are generated 
also depends on the choice of the starting numbers. 
And it is found that different generators prove 
to be optimal for different types of computers. 

Programs involving Monte Carlo computa- 
tions often have to address a random-number ge- 
nerator an enormous number of times (tens and 
hundreds of millions of times). Consequently, 
high speed is one of the most important character- 
istics of a generator. 

The generator recommended for a BESM-6 
computer has k = 517, c=0, m= 2%, and 
odd X,. This set of numbers does not meet the 
conditions of the above theorem (c = 0), and 
the period of this generator is less than m. How- 
ever, another theorem has been proved for gene- 
rators with c = 0, and by virtue of this theorem, 
the period of the recommended generator is 
288 = 2.75-1011, 


Determination of Percolation Threshold 
by Monte Carlo Simulation on a Computer. 
Distribution of Blocked 

and Nonblocked Sites 


Now we shall describe in detail a computer pro- 
gram that determines the percolation threshold 
by Monte Carlo techniques. Note that this pro- 
gram is not unique. In fact, each group of re- 


6—0240 


82 Part I. Site Percolation Problem 


searchers that deals with these problems pre- 
fers to work with their own program that has some 
individual features. This stems from specific 
features of different computers and to some ex- 
tent from the experience accumulated by indi- 
vidual programmers. : 

The problem that we mean is the site problem, 
and for the sake of simplicity, we consider only 
a two-dimensional square lattice. Actually, it 
will become clear later that a generalization of 
the method to an arbitrary lattice of arbitrary 
dimensionality is readily obtainable. 

Let us look at percolation in a square with the 
side containing L sites, so that the total number 
of sites is f° = L?, We assume the distance’ be- 
tween the sites to be unity, and describe the 
sites by their coordinates X and Y. For instance, 
a site with coordinates X = 9, Y = 25 is the 
site in the ninth column from the left and in the 
twenty-fifth row from below. 

In order to study percolation, we must fix 
which of the sites are blocked and which are 
not, and we must be able to vary the number of 
blocked sites in order to exceed the percolation 
threshold. To achieve this, we first assign a speci- 
fic number V to each site. Each site being charac- 
terized by two coordinates X and Y, this is equiv- 
alent to introducing a function of two variables 
V (X, Y) whose arguments X and Y do not run 
through all possible values but are allowed to 
assume only integral values within the interval 
from 1 to L. Programmers call such a function 
a two-dimensional array, and the values taken 
on by this function are said to be the elements of 
this array. For instance, the element V (34, 97) 
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of an array is a number assigned to the site with 
coordinates X = 31, Y = 97. Altogether the 
array Vhas L x L = f° elements, and we have 
to reserve in the computer memory a space 
necessary to store f° numbers. 

The work of the program begins with generat- 
ing this array. Its elements are random numbers 
uniformly distributed between zero and unity. 
A random-number generator outputs a number 
y, and this number is assigned to the array ele- 
ment V (1, 1). This means that this number is 
written into the corresponding memory cell of 
the computer, and that from this moment on the 
computer “memorizes” that V (1, 1) =y. The 
next number output by the generator is assigned 
to the element V (1, 2), and so forth. Thus all 
the elements of the array V are generated. 

Then a second two-dimensional array, that we 
denote by K, is formed. The elements of this 
array are zeros and unities, so that if, for exam- 
ple, K (25, 16) =0, this means that the site 
with coordinates X = 25, Y = 16 is blocked, 
but if K (25, 16) = 4, this site is nonblocked. 
The array Ķ is generated by using the array V 
and a certain number ¢ that falls within the in- 
terval from zero to unity. Varying t, we can 
change the number of blocked sites. 

The array K is generated by obeying the fol- 
lowing rule. Choose a site with coordinates X 
and Y. If V(X, Y)<t, then X (X, Y) = 1, and 
if V (X, Y) >#, then K (X, Y) = 0. The site 
with coordinates X and Y is treated as nonblocked 
in the former case and as blocked in the latter 
case. The quantities V being uniformly distri- 
buted within the interval from zero to unity, 


6e 
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we can assume ¢ being close to zero and obtain 
that almost all sites are blocked. Conversely, 
if ¢ is close to unity, almost all sites are non- 
blocked. When £ = 1/2, the number of blocked 
and nonblocked sites must be nearly 
equal. 

Making use of the distribution of random num- 
bers output by the generator, it is possible to 
relate the quantity £ to the mean fraction x of 
nonblocked sites obtained as a result of the above- 
described procedure. It can be shown (see Exer- 
cise 8) thatt = z. However, this equality is val- 
id if we take a very large number ff of sites or 
if we generate many arrays K with the same t 
and then average the fractions of nonblocked 
sites found in each array. Actually, in each con- 
crete array, z may deviate from ¢ to some extent, 
in both directions, but the equality holds the 
better, the larger Lf is. 

The computer thus stores the array V, from 
which we can generate an array K that describes 
which site is blocked and which is not. The form 
of the array K is dictated by the quantity ¢ which 
is approximately equal to the fraction of non- 
blocked sites generated in this array. Smoothly 
varying t, we can generate the distributions of 
blocked and nonblocked sites with smoothly 
varying concentration of nonblocked sites, z. 


Exercise 


8. Prove that the mean fraction z of nonblocked 
sites in an array K equals t. 
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Search for Percolation Path 


Let us assume that the array V is stored, the 
criterion ¢ is fixed, and the array K with a certain 
fraction of nonblocked sites is found. 

Now the computer knows precisely which site is 
blocked and which is not, and we pass to the 
second stage of the program, namely, the search 
for percolation path. Assume that we look fo 
percolation from left to right. First we replace 
all unities in the leftmost column (X = 1) by 
twos. The replacement consists in erasing a unity 
in the memory cell corresponding to a given ele- 
ment of the array K and in writing a two instead. 
A list of coordinates of the sites equipped with 
twos is compiled in the computer memory. Then 
the computer analyzes each site of this list. 
The computer calculates the nearest neighbors of 
the site under study and requests from the array 
K the information concerning these neighbors. 
If the nearest neighbor contains a unity, it is 
given a two, and its coordinates are added to 
the new list. After the first list has been analyz- 
ed, the computer memory contains a list of “sec- 
ond-generation” twos, that is, the list of unities 
replaced by twos because they were in contact 
with the first-generation twos. 

In order to economize on the computer memory, 
the first list is erased at this stage: it is no longer 
necessary; the appropriate memory cells are 
emptied. The computer starts an analysis of 
the second list and the formation of the third- 
generation twos. When this has been completed, 
the second list is erased and an analysis of the 
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third list is initiated. This analysis generates 
the fourth list, and so on. ; 

The number of twos in the array K grows in 
the course of this procedure. Twos are nonblocked 
sites linked by percolation to a nonblocked site 
of the leftmost column, that is, the twos label 
the percolation path. 

The search for percolation path terminates in 
two cases: 

1. A two appears on the right-hand side of a 
square. The computer recognizes the existence 
of percolation at the given value of t. 

2. There are no twos on the right-hand side of 
a square, and the analysis of the last list has not 
generated new sites marked with 2. This means 


that all paths broke down, and there is no per- 
colation at this ¢, 


Determination of the Threshold 


Assume that the Computer recognized percola- 
tion ata given ¢. Then it diminishes ż, and, mak- 
ing use of the same array V, finds a new array K 
with a reduced number of nonblocked sites. 
The search for percolation paths is then repeated. 
If percolation is found again, ¢ is further reduced, 
until at a certain t percolation ceases. Then the 
interval between this value of t and the minimum 
value at which percolation was found is divided 
in two, and percolation paths are looked for at 
this intermediate value of t. If it is found that 
percolation paths are cut, the interval between 
this last value and the minimum value at which 
percolation was detected is again divided in two. 
If percolation is found, the interval to be divid- 
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Fig. 15. Pattern of the distribution of zeros, unities, and 
twos at the moment of percolation onset. The broken line 
shows the path by which twos “percolated” from the left- 
hand side to the right-hand side of the square. In this 
case the computer did not stop calculations when the 
first “2” appeared on the right-hand side, but continued 
until new twos ceased to appear. 


ed in two is the one between the last value of t 
and the value at which no percolation was found. 

Percolation threshold is thus “bracketed” 
with an arbitrary precision. If we find no percola- 
lion at the first chosen value of ¢, then £ must be 
increased until percolation is detected, and this 
must be followed by “bracketing”. This method 
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enables us to find ż, corresponding to the perco- 
lation threshold, with any prescribed accuracy. 
At this t the computer finds the fraction of non- 
blocked sites, z, which is close to the value of t 
but is not necessarily equal to it. This value of 
x is taken for the percolation threshold found 
in this experiment. Figure 15 shows the first 
percolation path from left to right found in a 
30 X 30 square. 


Identical experiments are then rerun m 


any times 
with different 


sets of random numbers in the 


rmine the mean 
hold for a prescrib- 
is we have to add 
percolation thre- 


value xo(.f*) of percolation thres} 

ed number of sites, S. (For th 

up all the obtained values of 

shold and divide the sum by the number of runs.) 

In order to find the true percolation threshold 

Ta = lim zo (SP), it is necessary to vary the numb- 
N> 


er ofsites, b°, inthe square and determine the func- 
tion ze (4°). This function must then be approxi- 
mated by an analytical expression of the type* 


Xe (M) = ze (00) + S (5) 


that is, we have to ch 
Ta (œ), D, and y such th 
the best way the results 
ee 


oose three quantities, 
at expression (5) fit in 
oblained by means of 


uestion why this 
in the power-law manner, but the 
i iments show that this 
lation theory problems, 
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the computer. If this fit is achieved so that y > 0, 
we can say that it is the quantity x, (co) that 
gives us the limiting value of xe. Indeed, by vir- 
tue of expression (5), 

lim x, (f°) = 2¢ (00) 
MS +00 
The accuracy achieved by this procedure is the 
better, the greater the amount of data necessary 
to establish the dependence of ze (4°). In its turn, 
this amount of data is limited by the speed and 
memory of the computer employed. 


Exercise 


9. Look carefully at Fig. 15 and reconstruct the 
way by which individual groups of twos were 
formed. 


Part II 


Various Problems 
of Percolation Theory 
and Their Applications 


Chapter 5 ! ; 
Problems on Two-Dimensional Lattices 


We Are Planting an Orchard 
(the Bond Problem) 


Imagine that a vast orch 


hard is being designed. 
The fruit trees in t 


he orchard must be planted 
not arbitrarily but in aregular manner. They are 
to be located at the sites of a periodic lattice 
drawn on the ground. Many such lattices can be 
invented, but the following three will be suffic- 
ient for us here: a Square, a triangular, and a 
hexagonal lattices (the latter is often referred 
to as the “honeycomb lattice”). The lattices are 
shown in Fig. 16, Land being expensive, it is 
natural to try and plant the trees as close to one 
another as possible, but this cannot be done for a 
number of reasons, One of the reasons is that 
the designers are afraid of infectious diseases of 
the trees. Let us assume that experts on tree 
infections supplied the following information*®: 


* The author is not responsible for this information, and 
thus prays the reader not to be too serious about practical 
conclusions drawn from the solution of the problem formu- 
lated above. The problem is given mere y to illustrate 
the potential inherent to percolation theory, 
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(a) (b) 


(c) 


Fig. 16. Plane lattices: (a) square; (b) triangular; (c) ho- 
neycomb. 


1. A tree blighted by the disease can infect 
only its nearest neighbors. 

2. Some nearest-neighbor trees inevitably in- 
fect one another when one of the trees has been 
infected. In other cases the disease is not trans- 
ferred (e.g. this may depend on the separation 
between the branches of a given pair of trees). 
A pair of trees in which the infection is inevit- 
able will be referred to as a connected pair. 

3. The experts also supply us with the form of 
the function x (a), mainly, the probability for 
a pair of nearest neighbors chosen at random to 
be connected. This probability depends on the 
distance a between the nearest neighbors in a 
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given lattice. Obviously, the function x (a) 
increases with decreasing argument a: two trees 
infect each other the easier, the smaller the se- 
paration between the trees. 


We are to answer the following question: What 


ase to a concrete number 
of trees in the orchard? 


is convenient to switch to 
the cluster language introduced in the preceding 
chapter. We assume that two neighbor sites 
i ing a connected pair are 
ted by a piece of wire 
connecting the two sites. If two nearest-neighbor 


cted pair, the bond be- 


are linked by an unbroken bond 
linked by an unbroken 
nearest neighbors and are linked 
bonds (e.g. the sites 4 
siles C and D, in Fig. 17 are connected), 


An ensemble of Connected sites will be said 
to form a cluster. In the conte 


problem, the most import 
is thal a diseased lre 
its cluster and none outside this cluster, 

By definition, the fraction of unbroken bonds 
equals z. Further argument runs as in Chapter 3, 


xt of the given 
ant property of a cluster 
e infects all the trees of 
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When x is small, unbroken bonds are mostly 
single, almost all clusters consist of two sites, 
three-site clusters are infrequent, and four-site 
clusters are even less so. When æ is large, there 
exists an infinite cluster of connected sites. When 
x = 1, this cluster comprises all the sites in the 


a alee ay 


b+ 


L, cerei -l 


Fig. 17. Fragment of a square lattice with broken bonds. 
Shown in the figure are three two-atom clusters (1, 2, 3), 
one four-atom cluster (4), one six-atom cluster (5), and 
one ten-atom cluster (6). 


system. As x diminishes, this cluster loses some 
of the sites, and finally, at a certain critical 
value zo, the infinite cluster ceases to exist. 

The infinite cluster is precisely the catastrophe 
from which the orchard must be saved. Let 
P» (£) be the probability for a randomly chosen 
site to belong to the infinite cluster. If £ < Te, 
so that P> (x) = 0, one blighted tree can infect 
only several other trees. And if 22> 2; one 
blighted tree infects an infinite number of trees 
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in an infinite orchard with a probability P» (x). 
Therefore, if £< ze, the pocket of disease ac- 
cidentally brought to the orchard remains local- 
ized in the vicinity of the infected site, but it 
propagates all over the orchard if T >: 

In order to arrive at practical recommenda- 


between the trees exceeds a, otherwise the in- 


fection propagates throughout the orchard. 
Percolation theory 


ciated problem of determining z, as the bond 


square, T for triangular, and 
or honeycomb. Then Zs (H) der 
tion threshold of the sit 
lattice, zp (T) the per 
bond problem o 
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The function PP (x) introduced in this section 
for the bond problem should not be confused 
with the function P (x) defined earlier for the 
site problem. 

The bond problem can be formulated not only 
in cluster terms but also as a problem of percola- 
tion from one side of a square to the opposite 
side. Recall the experiment with a wire mesh 
with which we started the book. Some readers 
may have wanted to ask why it was necessary 
to cut simultaneously the four wires entering 
each site, instead of cutting randomly chosen 
individual wires (bonds). Now it can be clearly 
understood that by cutting the bonds the research- 
ers would find zy (S) instead of zs (S) which was 
actually determined in their experiment. Now 
we can explain why the site problem has been 
initially selected: it will be shown below that 
on a square lattice the bond problem has an 
exact analytical solution which yields that 
zy (S) = 0.5. Consequently, it would not be 
reasonable to carry out such a time- and labor-con- 
suming experiment for the sake of x (S), with 
x, (S) being known only from approximate solu- 
tions. 


Exercise 
1. Find the function P? (x) for 1 — £ &« 1 for 
the three lattices shown in Fig. 16. 


Inequality Relating zp to 2s 


When analyzing the site problem, the function 
P (x) is often replaced by a function Ps (x) 
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related to P (x) by a relation 
P (x) = aPs (z) (1) 


By definition, P (æ) is the probability for a ran- 
domly chosen site to belong to an infinite cluster. 
It can be written as the product of the probabili- 
ties of two independent events. In the language 
of the ferromagnetics problem, the first of these 
events is for a randomly selected site to be magnet- 
ic. The probability of this event equals x (see 
Exercise 1 to Chapter 1). The second event con- 
sists in this site being connected to an infinite 
cluster of magnetic sites, Therefore, the function 
Ps (x) defined by formula (1) is the probability 
fora randomly selected magnetic site to be con- 
nected to an infinite cluster. In other words, 
Ps (x) is the fraction of magnetic sites helonging 
to the infinite cluster, that is, the ratio of the 
number of sites belonging to the infinite cluster 
to the number of magnetic sites. We remind the 
reader that P (x) is the ratio of the number of 
sites belonging to an infinite cluster to the total 
number of sites. Naturally, the function Ps (x) 
grows monotonically with increasing x, equals 
unity for z = 1, and vanishes for US Ay 

The English mathematici 


an Hammersley, who 
was the first to enunciate 


percolation theory, 
proved a theorem that states that 


PS (x) < P» (a) (2) 
Both functions Ps and pb increase monotonically 
with increasing argument x. Therefore (Fig. 18) 
formula (2) implies that 


ThS Ts 


, 


(3) 
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that is, the threshold for the bond problem is 
not greater than that for the site problem on any 
lattice (not necessarily two-dimensional). This 
result can be rewritten as a different inequality: 


1— z> 1 — ts (4) 


that permits the following interpretation. Let 
us assume that we need to block electric current 


pe 
P5 
Xb Xs X 


through a network of wires or block the flow of 
liquid through a network of pipes, and this can 
be done either by blocking network sites or by 
cutting network bonds (wires or pipes) connect- 
ing these sites. Inequality (4) means that the 
system is more readily blocked by cutting out 
its sites than by cutting the bonds. The fraction 
of blocked sites that terminates the flow is less 
than the fraction of broken bonds giving the 
same result. This conclusion appears quite natur- 
al because not one but all bonds entering a site 
are cut when this site is blocked. 


ps p> 


Fig. 18. 


7—-0240 
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Exercise 


2. Find the function Ps (x) for 1 — z & 1 for 
the three lattices shown in Fig. 16. Compare the 
result with that obtained in Exercise 1 and check 
whether inequality (2) holds. 

Clarification. In Exercise 1 to Chapter 3 we 
recommended that you should find the function 
P (x) for the site problem, for 1 — z < 1. How- 
ever, there we meant only the principal term of 
the function, that is, P (x) = x. If this expression 
is substituted into formula (1), we obtain Ps = 14. 
This is the correct result in the sense that 


lim Ps (x)= 1. We suggest that the reader 
amd 


find the small terms that make the function 
Ps (x) differ from unity. Obviously, as «— 4, 
these terms have zero for their limits. Conse- 
quently, the result may be written in the form 


Ps (z) = 1 — A (1 — 3)" 


where A and n are positive numerical coefficients 
depending on the type of the lattice, 


Covering and Containing Lattices 


The site problem is more general than the bond 
problem. The bond problem is reducible to the 
site problem but on a different lattice said to 
cover the lcrmer. A covering lattice is constructed 
by using the following procedure: 

1. Place a site of the covering lattice in the 
middle of cach bond of the initial lattice. 

2. Connect two sites of the covering lattice 
if and only if the bonds of the initial lattice on 
which these sites were placed meet al a site of 
the initial lattice. 
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The result of such a construction is a new per- 
iodic lattice which is said to cover the initial 
lattice. 

Figure 19 shows the covering lattice in the 
case of an initial square lattice. Thin lines trace 
the initial square lattice. Semicircles show the 


Fig. 19. Covering lattice for the square lattice. 


places where the initial lattice has sites. The 
covering lattice consists of thin and solid lines, 
but there are no sites of the covering lattice 
where thin lines cross. The sites lie only at the 
crossings of solid lines and are marked by black 
circles. 

Each bond of the initial lattice joins three other 
bonds at one end and three more bonds at the 
other end. Hence, each site of the covering lat- 
tice must be connected with six other sites. This 
is shown in Fig. 19. Each site is connected with 
four other sites by solid lines, and with two more 
sites by thin lines. 


qr 
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Now assume that a bond problem has been 
formulated on the initial lattice, that is, a cer- 
tain fraction of randomly selected bonds has 

cut. 

ae ree now that ifa bond of the initial lattice 
is broken, the site of the covering lattice placed 
at this bond is blocked. We thus, obtain a site 
problem on the covering lattice. Its sites are 
randomly blocked, and the fraction of blocked 
sites equals the fraction of broken bonds in the 
initial lattice. 

Note that the existence of an infinite cluster 
of connected sites in the bond problem inevitably 
signifies the existence of an infinite cluster made 
of coupled unbroken bonds. Conversely, the ab- 
sence of an infinite cluster of sites signifies that 
bonds do not form an infinite cluster, 

As follows from the method of construction 
of the covering lattice, the existence of an infin- 
ite cluster formed by unbroken bonds in the 
initial lattice signifies that there exists an in- 
finite cluster formed by nonblocked sites in the 
covering lattice, and conversely, if there is no 
infinite cluster of bonds in the initial lattice, 
there is no infinite cluster of sites in the cover- 
ing lattice. Therefore, the percolation threshold 
of the bond problem in the initial lattice equals 
the percolation threshold of the site problem in 
the covering lattice. If the initial lattice is 
denoted by L, and the covering lattice by L 
this statement can be written as formula 


Lp (L) = zs (Leov) 


Introducing the conce 
another lattice, it is pos 


coy? 


(5) 
pt of a lattice containing 
sible to derive a number 
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of inequalities relating percolation thresholds on 
a number of lattices. Assume, for instance, that 
lattice L is obtained from lattice Leone by cross- 
ing out a certain number of bonds. Then lattice 
Leont is said to contain lattice L. 

For instance, let us take a triangular lattice. 
If we erase all the bonds that are marked in 


Fig. 20, The triangular lattice contains the square lattice. 


Fig. 20 by two short lines, it will transform into 
the lattice shown on the right-hand side of the 
drawing. It can be easily seen that from the stand- 
point of the site or bond problems this new lattice 
is equivalent to a square lattice. Indeed, the fact 
that the angles between the bonds of the new 
lattice are not 90° is found to be immaterial if 
we analyze the relations between different sites 
(the new lattice can be just “straightened out”). 
The percolation threshold of the bond problem 
(as well as of the site problem!) on this lattice 
is exactly equal to the percolation threshold on 
the square lattice. For this reason, a triangular 
lattice is said to contain a square lattice. 

Now let us assume that a certain fraction of 
bonds in the containing lattice is broken. The 
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bonds of the containing lattice can be classified 
into the bonds that are common for the contain- 
ing and contained lattices, and the bonds that 
are specific for the containing lattice (these last 
are marked by two short lines in Fig.320). The 
bonds being broken quite randomly, the frac- 
tion of broken bonds in one category of bonds is 
absolutely the same as in the other category, and 
equals the fraction of broken bonds in the whole 
lattice. Therefore, in order to obtain a contained 
lattice with the same fraction of broken bonds, 
it is necessary to break additionally those bonds 
of the containing lattice that are intact but are 
Specific for this lattice, that is, are completely 
absent in the contained lattice, 

This argument demonstrates that the number 
of unbroken bonds leaving each site of the con- 
taining lattice is not less than (is greater than 
or equal to) the number of unbroken bonds leav- 
i f the contained lattice. As 
lity for a randomly selected 
an infinite cluster of the con- 
ot less than that of the contain- 
atement implies the inequality 


PE (E)K P? (2) 


The left-hand side of inequality (6) includes the 
function p» (x) calculated for the contained 
lattice, and the right-hand side includes this 
function for the containing lattice, By analogy 
to inequality (2) employing inequality (3), ine- 
quality (6) implies that 


(6) 


ty (Leont)<ay (L) (7) 
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Ht J 


a ig e 


| 
Fig. 22. Transformation of one unit cell of the lattice 


shown on the right-hand side of Fig. 21 into a honeycomb 
unit cell. 


that is, the containing lattice has a lower 2p 
than the contained lattice. 

It has been already mentioned that a triangular 
lattice contains a square lattice. For this reason, 


zy (T) < ty (8) © 


Now assume that some bonds are deleted in 
a square lattice, as shown in Fig. 24. This yields 
the lattice shown on the right-hand side of the 
drawing. Look attentively at this drawing. It 
is equivalent to a hexagonal lattice. Just pull 
it slightly upwards, somewhat warping the bonds, 
and its cells (Fig. 22) will transform into the 
“honeycomb” lattice shown in Fig. 16. Conse- 
quently, a square lattice contains a hexagonal 
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lattice, and hence, 


z» (S) < z» (H) 9) 
Inequalities (8) and (9) imply that 
zp (T) < z» (H) (10) 


Let us turn now to the site problem. Assume 
that the same sites are blocked in the contain- 


bonds can 
only facilitate its formation. But if it is known 
exists in the containing 
on of sites, no conclusion 
can be drawn from this as to the existence of an 


infinite cluster in the contained lattice. The 
removal of some bonds between sites in the tran- 
sition from the containing to the contained lattice 
can be “lethal” for the infinite cluster, Conse- 
quently, the percolation threshold of the site 
problem in the contained lattice cannot be small- 
er than that in the Containing lattice: 


Ts (Leont)< x, (L) 


Therefore, we 
for the triangul 


(11) 
can write a chain of inequalities 
ar, Square, and hexagonal lattices: 
zs (T) < zs (S) < z, (H) (12) 


quite identical to that for the bond problem. 
Let us turn again to Fig, 19 where the cover- 
ing lattice for the Square lattice is shown. Imag- 
ine that the bonds traced by thin lines have been 
deleted, Obviously, this gives us a square lattice 


Ch. 5. Problems on Two-Dimensional Lattices 105 


merely rotated through 45°, which, of course, 
is quite immaterial for percolation theory prob- 
lems. 

The covering lattice for a square lattice thus 
contains a square lattice. By virtue of formula 
(5), 

Ts (Loov) = Th (S) (13) 


where Loy is interpreted as the covering lattice 
shown in Fig. 19. However, inequality (11) and 
the fact that this covering lattice contains a square 
lattice imply that 


Ts (Leov) < ws (S) (14) 
Inequalities (13) and (14) imply that 
T (S) < as (S) (15) 


We have thus obtained that the bond problem 
threshold on a square lattice is lower than the 
site problem threshold. Inequality (15) is a parti- 
cular case of Hammersley’s general theorem that 
is embodied in formula (3) and was given earl- 
ier without proof. This theorem was not used in 
deriving (15), and we thus can say that the argu- 
ments given above prove this theorem in the case 
of a square lattice. 


“White” Percolation and 
“Black” Percolation 


Now let us look at the bond problem from a slight- 
ly different viewpoint. So far we normally have 
said that there are intact and broken bonds distri- 
buted randomly over a lattice and applied the 
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term “cluster” to a ensemble of sites connected 
en bonds. 
by ga aeai can be reformulated in a more 
symmetric manner. Let us rename broken bonds 
“black”, and unbroken bonds “white”. An ensem- 
ble of sites connected by white bonds will be re- 
ferred to as a white cluster, and an ensemble of 
sites connected by black bonds as a black cluster 
(in our former terminology it wasa “white” cluster 
that we called a cluster). Let the frac- 
tion ofwhite bonds be denoted, as before, by x. 
The fraction of black bonds will be denoted by 
q. Each bond has to be either black or white, 
so that q = 1 — z. 

After this reformulation we can speak of per- 
colation through both white bonds and black 
bonds. 

When the concentration x of white bonds is low, 
there is no infinite white cluster, but there is 
an infinite black cluster, that is, an infinite clust- 
er of sites connected through black bonds. Con- 
versely, when the concentration q of black bonds 
is low (i.e. at x nearly equal to unity), there is 
an infinite white cluster but not an infinite 
black cluster. 

When x varies from zero to unity, two events 
take place: the black infinite cluster vanishes, 
to be replaced by a white cluster, or, to say the 
same in different terms, percolation through black 
bonds vanishes, to be replaced by percolation 
through white bonds. But what is the sequence 
in which these events take place? 

The white and black bonds differ only in labels. 
It is obvious, therefore, that the critical concen- 
tration qp at which percolation appears through 
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black bonds equals the concentration x, at 
which percolation appears through white bonds. 

Therefore, as x increases, percolation through 
white bonds appears when z = £p, and _percola- 
tion through black bonds stops when z= 
= 1 — qp = 1 — xp. The sequence in which 


Il- l 
IN, 
as 
0 SEn 0.5 Xp Doe: 
(u) 


Fig, 23. 


these events take place depends on the sign of 
the difference x, — 0.5. 

If it is necessary for percolation through white 
bonds that the bonds be more than 50% white 
(this means that the number of white bonds must 
exceed that of black bonds at the percolation 
threshold), then at increasing x first percolation 
through black bonds ceases and then percolation 
through white bonds appears (Fig. 23a). In region 
I in Fig. 23a there is percolation only through 
black bonds, in region JJI only through white 
bonds, and in region IJ there is no percolation 
through either white or black bonds. 
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If zp < 0.5, first percolation through white 
bonds appears and only then percolation through 
black bonds ceases. In region J in Fig. 23d there 
is percolation through black bonds, in region III 
it proceeds through white bonds, and in region 
II through both black and white bonds. 

The same symmetric reformulation can be 
found for the site problem. We remind the reader 
that in this case all bonds are unbroken, but 
sites can be of two sorts. In the wire-mesh prob- 
lem we called them nonblocked and blocked. In 
the ferromagnetics problem these were magnetic 
and nonmagnetic sites. And now we introduce, as 
in the bond problem, a universal notation: 
nonblocked, or magnetic, sites are said to be 
white, while blocked, or nonmagnetic, siles are 
said to be black. White sites are said to be connect- 
ed if they are nearest neighbors or are connected 
by a string of white sites that are formed by pairs 
of nearest neighbors. Black sites can be connect- 
ed in the same symmetric manner. 

We are justified in speaking of percolation 
through white and black sites. If x, > 0.5, there 
is a region of z values with no percolation through 
either white or black sites (l—a,<2< z). 
If z3< 0.5, the region ts < x< 1 — zr, is 
the region in which percolation proceeds both 
through white and black sitos. 

The symmetric approach proves constructive 
because sometimes it is possible to conclude from 
an analysis of the lattice that it cannot accommo- 
date either percolation through white or percola- 
tion through black, or that one of these must ne- 
cessarily be present. This information is suffic- 
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ient to make readily certain conclusions on per- 
colation threshold. 

For instance, let us consider the site problem 
on a triangular lattice. Assume that the lattice 
percolates through white sites. You will readily 
see that this excludes percolation through black 
sites. Let us assume that we study percolation 
of current from left to right on a very large net- 
work, as in the experiment whose description 
opened this book. But now the network is not 
a square but a triangular lattice. It is easy to 
understand that the structure of any triangular 
lattice is such that a rightward percolation path 
through white sites excludes any possibility for 
downward percolation through black sites. In- 
deed, black sites cannot “penetrate” through the 
broken line traversing the entire network from 
left to right and connecting white sites. 

As we explained in Chapter 3, an infinite cluster 
ensures percolation in any direction if the size 
of the system is sufficiently large. It then follows 
that a triangular lattice cannot contain, for the 
same value of z, two coexisting infinite clusters, 
one of white and another of black sites, that is, 
percolation through white sites and percolation 
through black sites are incompatible. Hence, 


s (I) > 0.5 
The same conclusion follows for square lattices: 
za (S)> 0.5 


A theorem has been proved for triangular lat- 
tices, stating that zs (T) = 0.5. The proof of 
this theorem cannot be given here, but the main 
points can be grasped more or less easily, making 
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use of the concepts of percolation through white 
and percolation through black. Drawing various 
configurations of black and white sites, it can 
be noticed that blocked rightward percolation 
through white sites necessarily implies downward 
percolation through black sites (square lattices 
do not possess this property!). 

A triangular lattice thus cannot simultaneously 
have percolation through white sites and percolation 
through black sites, but there must always be perco- 
lation through one sort of sites. This means that 
region II in Fig, 23 (a, in this figure must now 
be replaced by 2) degenerates to a point, that 
is; 2a (T) = 0.5, 

In the cases of bor 
analysis can be f. 
lattices, 


nd problems of this type an 
acilitated by operating with dual 


Dual Lattices 


Only plane lattices can be dual. By definition, 
a plane lattice is a lattice that can be placed on 
a plane, with bonds intersecting only at points 
where lattice sites are located. For instance, all 
the lattices shown in Fig. 16 are plane, while 
the covering lattice in Fig. 19 is not plane be- 
cause its bonds intersect at points shown by arcs, 
but there are no sites a 


t these crossing points. 
(These arcs act as “bridges” that uncouple right- 
ward and downward roads.) 


Each plane lattice divides the plane into cells. 
Lattice Là is said to be dual to lattice L if each 
bond of L4 intersects one and only one bond be- 
longing to L, and vice versa, each bond of L 
intersects one and only one bond belonging to L4, 
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Fig. 24. Construction demonstrating that square lattices 
are self-dual. The solid lines and black circles show the 
bonds and sites of the initial lattice; the dashed lines 
and open circles show the bonds and sites of the dual 
attice. 


Fig, 25. Construction demonstrating that a triangular 
lattice is dual to a honeycomb lattice (and vice versa). 
Notations are the same as in Fig. 24. 
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Furthermore, each cell of lattice L must con- 
tain only one site of lattice Zà (and vice versa). 

You see from the definition that duality is a 
mutual property: if Lê is dual to L, then L is dual 
to L4. Figures 24 and 25 show that a square lattice 
is dual to a square lattice, and a triangular and 
a hexagonal lattices are mutually dual, 

Let us return now to the bond problem. We 
assume that if a bond in the initial lattice is 
white (unbroken), the bond of the dual lattice 
intersecting the former is black (broken). There- 
fore, if the concentration of white bonds in the 
initial lattice is z, the concentration of white 
bonds in the dual lattice is q=1—r. 

In what follows it will be more convenient to 
interpret percolation threshold as that concen- 
tration of white bonds at which conductivity 
from left to right first appears (or disappears) 
in a very large wire mesh with electrodes soldered 
to it (see Fig. 1). This formulation of the problem 
is identical to that discussed at the very begin- 
ning of the book, but now the disrupted elements 
of the net are the bonds and not the sites of the 
lattice. 

Assume now that a rightward percolation path 
exists through white bonds of the initial lattice, 
It is easy to understand that this signifies no 
downward percolation through white bonds in 
the dual lattice. Indeed, by definition, a white 
bond of the initial lattice is intersected only 
by a black bond of the dual lattice. H 


Tence, if 
the initial lattice contains a nowhere disconnected 


broken line composed of white bonds and tra- 
versing the entire network from left to right, this 
means that the white bonds of the dual lattice 
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cannot at any point “break through it” and form 
a downward path. 

It a lattice is sufficiently large, the existence 
of a rightward path on the initial lattice signifies 
that the fraction of white bonds exceeds the thre- 
shold value: 


«> sp (L) (16) 


Conversely, the absence of a downward path 
through white bonds of the dual lattice indicates 
that the fraction of white bonds, q = 1 — x, 
in the dual lattice is less than the threshold 
fraction for the dual lattice, that is, 1 — z < 
< zp (L8), or 

xz > 1 — x, (Lì) (17) 
According to all this, all the values of x satis- 
fying inequality (16) also satisfy inequality (17). 
Hence, xp (L)> 1 — zp (L8), or 

£» (L) + z» (L4) > 1 (18) 


On a square lattice L = L4, so that inequal- 
ity (18) implies that 
ay (8) 0.5 (19) 
and making use of inequality (15), we addition- 
ally obtain 
zs (S) > 0.5 (20) 
that is, we come to the same conclusion as in 
the preceding section. Recall that the experi- 


ment with a square lattice described in Chapter 
1 gave a, (S) = 0.59, which does not contradict 


inequality (20). 
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It has been rigorously proved for the square 
and triangular lattices that inequality sign in 
formula (18) must be replaced by equality sign, that 
is, 


zn (L) + xy (L3) = 1 (21) 
This immediately leads to two new results: 
zy (S) = 0.5 (22) 

and 
zy (T) + x (A) = 1 (23) 


We shall not give here the rigorous proof of 
formula (24): it would require that we introduce 
a number of new concepts that will not be nec- 
essary in the further exposition. However, this 
formula can be given a sufficiently clear interpre- 
tation. Tracing out various configurations of 
white and black bonds, it can be discerned that in 
the case of square and triangular lattices the 
absence of rightward percolation through white 
bonds in the initial lattice always signifies that 
there exists a downward percolation path through 
white bonds in the dual lattice. Let us take this 
statement for granted. Assume that the fraction 
of white bonds is such that there is no rightward 
percolation on the initial lattice. For a suffi- 
ciently large network this means that 


&< Ly (L) (24) 
Correspondingly, the dual lattice contains a 
a downward percolation path. The fraction of 


white bonds in this lattice is 1 — z. Therefore, 
{= a2 a, (L3), or 


z< 1 — xy (LÀ) (25) 
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The above arguments signify that any value of 
w satisfying inequality (24) must also satisfy 
inequality (25). This is possible if 1 — sp (L4)> 


2 xy (L), or 
x (L) + ty (L9) < 1 (26) 


Inequalities (18) and (26) yield equality (21). 

Formula (23) alone is not sufficient for find- 
ing separately zp (T) and x, (H). However, if 
we make use of the star-delta transformation 
known in the theory of electric circuits, we arrive 
at a relation between percolation thresholds of the 
bond problem on the triangular and hexagonal 
lattices. As a result, each of these thresholds is 
calculated: 


zy (T) = 2 sin (7/18) ~ 0.347 296 
(27) 
zp (H) = 1 — ay (T) = 1 — 2 sin (7/18) ~ 


æ 0.652 704 


Exercise 


3. Look at Fig. 53 (see p. 252) illustrating 
the star-delta transformation. If you understand 
by looking at this drawing how to write the rela- 
tion mentioned above and arrive at formulas 
(27), you have done very well. If you do not go 
that far, do not despair because the problem is 
not that simple. The derivation of formulas (27) 
first obtained in 1963 by the English mathema- 
ticians Sykes and Essam was an important event 
in percolation theory. Now read attentively the 
text placed alongside the figure and you are likely 


g* 
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to feel pleased with the beautiful application of 
probability theory yielding such a nontrivial 
result. 


Results for Plane Lattices 


In conclusion we must give a summary table for 
percolation thresholds for plane lattices (Table 1). 


Table 1 
Percolation Threshold for Plane Lattices 
Lattice type Xb Xs 
SSS ell 


Triangular 0 
Square 0 
Hexagonal 0 


Only two figures in this table, n 


and a; (H), were obtained by appro 
ods. All the 


amely, x, (S) 
ximate meth- 


dimensional lattices is much less 
» hot a single exact solution 

You need not regard this 
o universal methods exist 
ms of percolation theory 
exact solution mentioned 
above looks like a miracle. Rather, we should 


be surprised how many exact solutions have 
already been devised. 
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Exercise 


4. Let us revisit the orchard mentioned at the 
beginning of this chapter. Assume that the dis- 
tance separating the trees is chosen so as to 
satisfy the requirement that the fraction of con- 
nected pairs equals the threshold value. Let 
the function a (x) giving the distance between 
neighboring trees as a function of the fraction 
of connected pairs, 2, be known. Naturally, the 
ereater x is, the shorter a is because the trees 
infect one another the easier, the closer to one 
another they are planted. If « = zp, the sepa- 
ration between trees equals a (xp). Find the 
area per tree for the thus chosen separation, for 
three different lattices. The lattice with the 
least area per tree is the most profitable. Could 
you say in advance which of the lattices gives 
the least area per tree, starting only with the 
information that the function a (x) monotoni- 
cally decreases with increasing <£? 


Directed Percolation 


Assume now that the trees planted at the sites 
of a plane lattice form not an orchard but a for- 
est, and that this forest catches fire. Some near- 
est-neighbor trees have intertwined branches and 
readily pass on the fire. In accordance with our 
general terminology, we shall say that the sites 
containing such trees are connected by white 
bonds. Other nearest-neighbor trees do nob light 
each other up. We shall say that the correspond- 
ing sites are connected by black bonds. The white 
and black bonds are distributed randomly over 
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the lattice, with the fraction of white bonds being 
equal to z. . my 

Our problem is to find the critical value ty 
such that at £z < x, the pocket of fire remains 
localized, while at z > z, the fire spreads over 
the whole forest. 

Obviously, this is simply one more example of 
the bond problem. The value of xy can be found 
by looking it up in Table 4 (see p. 116). 

Assume now that the forest fire is 
by a strong wind, so that the fl 
only along the wind. 
interesting problem, 
rected percolation. 

We postulate that the 
siles of a square lattice 
along the diagonal of tl 
the wind direction is i 


accompanied 
ames propagate 
This results in a new and 
called the problem of di- 


forest was planted at the 
, and the wind is blowing 
1e square cells. In Fig, 26 
ndicated by the arrow al 
the top of the figure, and the lattice has been 
rotated through 45° with respect to the standard 
arrangement of square lattices in this book. The 
problem is now reformulated as follows. Each 
white bond is converted into a vector whose 
arrow is placed so that the projection of the vec- 
lorontothe wind direction is Positive. As before, 
black bonds do not let the fire pass whatever 
the direction, while white bonds let it pass only 
in the direction of the arrow. We want lo find the 
critical fraction of white bonds beginning with 
which one tree in flames in an infinitely large 
forest can start a fire Spreading lo an infinitely 
large distance, 

The solid lines with arrows in Fig, 26 stand 
lor white bonds: two percolation paths are shown: 
Z and 2. Path 7 traces the locus of the fire, while 
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Wind 


vig. 26. Directed white bonds are shown by the solid 
lines with arrows; black bonds are shown by the thin lines. 
Flame propagates from left to right along path 7, and 
cannot propagate along path 2. Two segments of path 2 
on which the fire would have to move against the wind 
are marked by two dashes. 


path 2 does not correspond to the conditions of 
directed percolation: at two points marked by 
dashes the propagation goes counter to a while 
bond, that is, against the wind. 

We thus find that while in nondirected perco- 
lation white bonds worked in both directions, 
in directed percolation they work only in one 
direction. Hence, the critical fraction of white 
bonds in directed percolation, 2p, cannot be 
smaller than in conventional percolation, that 
is, ap > Th 

At the present time, approximate solutions 
have been obtained for a number of problems 
with directed percolation. Thus, in the above- 
described bond problem on a square lattice 
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zp = 0.63 or 0.64 (the results yielded by several 
methods slightly differ from one another). Recall 
that in the nondirected bond problem on a square 
lattice we had 2, = 0.5. 

Quite a few physical problems are reducible 
to directed percolation. As an example, we can 
cite the motion of electrons in a strong electric 
field in a randomly inhomogeneous medium. 
These are media with properties that vary ran- 
domly from point to point, so that a moving 
electron runs into obstacles that it has to bypass. 
At the same time, the uniform electric field acts as 
a wind that drives electrons in one direction only. 

Directed percolation also realizes in the prob- 
lem of the electric conductance of a wire nel 
(see Chapter 1) if we assume that each bond be- 
tween the sites of the net contains a diode (a 
rectifying element) that lets the current flow 
in one direction only. Then the fraction of broken 
bonds at which the current through the net drops 
to zero corresponds to the threshold of directed 
percolation. A combined problem, with diodes 


present not in each bond of the net, has also been 
analyzed, 


Chapter 6 


Three-Dimensional Lattices 
and Approximate Evaluation 
of Percolation Thresholds 


As was shown in Chapter 5, percolation threshold 
essentially depends on lattice type. The purpose 
of the present chapter is to explain qualita- 
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tively what exactly are the properties of lattices 
that are important for percolation thresholds. 
Having understood this aspect, it is possible to 
predict percolation thresholds (to within 10%) 
without solving the problem. The ability to come 
up with such predictions (be the prediction 
accurate or not) appears to be valuable because, 
first, the number of different lattices is quite 
substantial, and second, the computation of a 
threshold takes (on top of a high skill) about half 
an hour of computation time of the most advanced 
computers. 

Moreover, percolation problems are not re- 
stricted to lattice problems only. We shall be 
able to see that applications mostly require 
solving problems that are formulated not on 
lattices. It is found that the ideas of approximate 
evaluation presented in this chapter are quite 
fruitful for nonlattice problems. It was on the 
basis of these ideas that a number of _ percola- 
tion thresholds for nonlattice problems had been 
predicted with high accuracy long before these 
problems were solved with computers. 

In order to collect the necessary experience, 
we have to go beyond the framework of plane 
lattices discussed in the preceding chapter and 
turn to three-dimensional lattices. 


Three-Dimensional Lattices 

The simplest of the three-dimensional (or simply 
3D) lattices is the simple cubic laltice. 1t is shown 
in Fig, 12. Its unit cell is a cube shown in Fig. 27. 
The vectors aj, äp, a3 are called translation vec- 
tors. Elongating each of these translation vectors 
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an integral number of times (n,, ns, n) and then 
adding up the resultant vectors, we can obtain 
a vector Ra, n, n, beginning at the origin of 
coordinates and ending at any site of the simple 
cubic lattice: 


Rau, na, ms = 74a + ngas + Ngaz 


In cubic lattices the three vectors aj, ay, ay are 
of identical length a, so that the numbers n, 


£ 
/ 


Fig. 27, Unit cell of a simple cubic lattice. 


a; 


Na, Ng are simply three Cartesian coordinates of 
lattice sites in units of a. The simple cubic lat- 
liceissaid to be generated by translation (parallel 
displacement) of the cubic unit cell by vectors 
that are integral multiples of äl, ay, ay. 

The most important characteristic of a lattice 
is the number of nearest neighbors (also called 
the coordination number), denoted by z. In the 
simple cubic lattice z = 6, 

Such alkali halides as NaCl (common salt), 
KCI (rock salt), LiF, Nal, and a number of 
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others crystallize into the simple cubic lattice. 
The alkali metal ions (e.g. Na*) alternate in 
this lattice with the halogen ions (e.g. Cl-). 

Body-centered cubic lattice (abbreviated to bec). 
This lattice can be composed of two simple 


(a) 


Fig. 28. (a) Body-centered cubic lattice; (b) unit cell of 
the bee lattice. 


cubic lattices by displacing them relative to 
each other by one half of the body diagonal of 
the unit cell. (By definition, the body diagonal 
of a cube is the line connecting the opposite 
vertices of the cube and passing through the cen- 
tral point.) 

A bee lattice is shown in Fig. 28a. lts unit 
cell is shown in Fig. 28b. The entire hee lattice 
can be generated by translating this unit cell 
by vectors aj, äs, ay; however, it will not be suffici- 
ent to place the vector Rn, na n, dl the origin if 
we want lo generate all sites of the lattice. We 
also need to place the origin of the vector al 
the central atom. Only then will the two simple 
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cubic sublattices of the bee lattice be translated 
together. 

The nearest neighbors of each site of a bec lat- 
tice lie in the direction of the body diagonal of 
the unit cell (see Fig, 28b). Each site has eight 
such neighbors, so that z = 8. The distance to 
a nearest neighbor equals one half of the body 


Fig. 29. Unit cell of a face-centered cubic lattice. 


diagonal, that is, /Ba/2, where tho letter a 
denotes, as before, the length of the edge of the 
unit cell. 


Univalent alkali metals, such as lithium, sod- 
ium, potassium, rubidium, cesium, bivalent 
barium, and a number of other substances 
crystallize into bee lattices. 

Face-centered cubic lattice (abbreviated to fec). 
A unit cell of an fce lattice is shown in Fig. 29. 
In comparison with a simple cubic lattice, it 
contains additional sites placed at the centers 
of cach face. Tn order to translate this cube, the 
origin of the vector Ras, na n, must be placed 
at the origin of coordinates and also at the centers 
of the three nonopposing faces of the cubic unit 
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cell. The nearest neighbors of each site of an 
fee lattice lie in the direction of face diagonals. 
The distance to a nearest neighbor is a/ V 2. Each 
of the three mutually perpendicular plane in- 
tersecting at a given site contains four nearest 


(a) (b) 


Fig. 30. (a) Diamond lattice; (b) tetrahedral structure of 
bonds in the diamond lattice. 


neighbors of this site, so that the number of 
nearest neighbors, z, of a site is 12. 

Substances that crystallize into fcc lattices 
are such metals as copper, silver, gold, alumin- 
ium, lead. 

The last lattice that we introduce into our ana- 
lysis is the diamond lattice. It is shown in Fig. 30a. 
You obtain it if you imagine two face-centered 
cubic lattices displaced with respect to each other 
along the body diagonal of the cubic unit cell 
by one fourth of its length. 

Substances that crystallize into diamond-type 
lattices are elements of the fourth group of the 
periodic table: carbon (diamond), as well as the 
two most important semiconductor elements: 
germanium and silicon. All these elements are 
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telravalent, and the atoms in their lattices are 
bonded by covalent forces. To simplify, we can 
imagine that each atom has four “arms” corre- 
sponding to four valence electrons. An atom in 
a lattice then holds hands with its four nearest 
neighbors. The diamond lattice is perfectly 
suited to this type of bonding. Each site of this 
lattice lies at the center of a regular tetrahedron 


formed by other sites (Fig. 30b). The number of 
nearest neighbors is z = 4. 


Percolation Thresholds for 3D Lattices 


The bond and site problems are formulated for 
3D lattices exactly as for plane lattices. As there, 
we again assume that bonds connect only near- 
est-neighbor sites. 

Table 2 summarizes percolation thresholds 
of the site and bond problems in the 3D lattices 
described above. As we have already mentioned, 
no exact solution was obtained in the 3D case. 
All the results listed in Table 2 were obtained 
by various approximate methods that, as a rule, 


Table 2 
Percolation Threshold for 3D Lattices 
Lattice type | Xp | £s 
Simple cubic 0.25 0.31 
Body-centered cubic 0.18 0.25 
Face-centered cubic 0.12 0.20 
Diamond 0.39 0.43 
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employ computers. Naturally, certain slight 
discrepancies are found between the results 
published in the scientific literature. ‘Table 2 
is a selection of results that are the most reliable, 
from our point of view. 

Now the problem is to try and understand, by 
looking at this table and Table 4 (see p. 116), 
which summarizes the results for plane lattices, 
why some lattices have relatively high percola- 
tion thresholds, and why other lattices have 
relatively low thresholds. We begin with the 


bond problem. 


Factors Determining Percolation Threshold 
in the Bond Problem 


If all bonds are intact, each site is connected with 
z other sites, where the number z of nearest neigh- 
bors varies considerably among lattice types. 
At a fixed fraction æ of intact bonds, each site 
is connected, on the average, with zx other sites. 
Let us try and test the following hypothesis: 
Is it realistic for the quantity zx that gives the 
average number of sites with which each site 
is connected to carry sufficient information for 
deciding whether percolation is present or absent 
in the lattice? Can it be that we need no other 
information on the properties of a lattice but its 
coordination number z, and that percolation 
sets in in all lattices at the same value of zg? 
It is rather clear that this hypothesis cannot be 
very accurate. But is it approximately valid? 

This is very easy to prove or disprove. Let 
us find the product zzp for all lattices with known 
percolation thresholds of the bond problem. If 
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this product is universal, that is, identical or 
at least nearly identical for all lattices, the pro- 
posed hypothesis holds or holds approximately. 


Table 3 
The Product zzp for Different Lattices 
Lattice type z Xh 2Xh 

Plane lattices 
Square 4 ]|0.5 2.0 
Triangular 6 | 0.35.| 2.4 
Honeycomb 3 | 0.65 | 2.0 

3D lattices 
Simple cubic 6 | 0.251] 4.5 
Body-centered cubic 8 10.48] 1.4 
Face-centered cubic 12 | 0.42] 1.4 
Diamond 4 | 0.39] 1.6 


The required data are summ 


arized in Table 3, 
You see that in plane | 


attices we have 


Zy = 2 (1) 
with an error less than 10%, and in 3D lattices 
Sy = 4.5 (2) 


an number of 
t of percola- 
s hold approx- 
ach of the quan- 
Ties by at least a factor of two 
group of plane lattices and in that 
of 3D lattices, the accuracy with which the pro- 
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duct zx is constant within each group must be 
considered remarkably high. 

We conclude that percolation threshold of the 
hond problem can be approximately evaluated 
if we know the number of nearest neighbors, and 
use formula (1) for plane lattices and formula 
(2) for 3D lattices. Percolation threshold of the 
bond problem is most sensitive to the number of 
nearest neighbors, but is much less sensitive to 
all other properties of lattices (e.g. to the number 
of next-to-nearest neighbors, i.e. the sites 
that are second in their closeness to a given site 
after ils nearest neighbors). 

We have thus arrived at a very simple and at 
the same time a relatively accurate method of 
evaluating percolation thresholds of the bond 
problem applicable to arbitrary lattices. 


How to Evaluate Percolation Threshold 
in the Site Problem 


Let us analyze now a similar scheme for the site 
problem. It will be natural to test first the fami- 
liar version, that is, check the variation of the 
quantity zrg among different lattices. It will be 
readily found that the product varies almost as 
each of the quantities z and zs does independent- 
ly. One should not be surprised: the product 
z2,, has a well-defined physical meaning for the 
bond problem, namely, that of the mean number 
of intact bonds per site. In the case of the site 
problem, a bond is effective if it connects two 
white sites, but is ineffective in all other cases. 
Consequently, the product zr; is virtually 
meaningless. 
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Fig. 31. Construction of tangent circles in a honeycomb 
lattice. The lattice itself is shown in Fig. 16c. The circles 
have radii equal to one half of the distance to the nearest 
neighbor. The white sites correspond to white circles, and 
the black sites to black circles, The percolation paths 
through white circles are traced with the solid lines, 


In 1970 the American physicists I. Scher and 
R. Zallen proposed a different method for evalu- 
ating percolation threshold in the site problem. 
Their idea was to put each site in correspondence 
with a specific part of space. Then they went on 
to say that percolation through white sites appears 
when the fraction of Space occupied by these 
sites exceeds a certain critical value slightly 
dependent on the lattice type. 

Imagine around each lattice site a sphere (or a 
circle if we deal with a plane lattice) of radius 
equal to half the distance to the nearest neigh- 
bor. The spheres (circles) constructed around 
neighboring sites are then tangent to one another 
(Fig. 31). A white site is assigned a white sphere, 
and a black site a black sphere. If two white sites 
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are connected, there is a path between them 
through tangent white spheres (see Fig. 31). There- 
fore, the onset of percolation means the formation 
of infinitely long paths through tangent while 
spheres. 

Now we assume that percolation sets in when 
the fraction of the total volume (or of the total 
area) occupied by while spheres (or circles in a 
plane lattice) exceeds a certain critical value in- 
dependent of lattice type. In order to test this 
hypothesis, we need to calculate the fraction of 
Volume occupied by white spheres for x = 
for different lattices with known values of £s A 
then compare the results. 

First we need to calculate the fraction of volume 
filled up by white spheres for « = 1, that is, 
in the case when all spheres are white. This 
quantity is denoted by letter f and referred to as 
filling factor. The filling factor equals the fraction 
of volume occupied by the spheres constructed 
around each lattice site and having a radius equal 
to half the distance to the nearest neighbor. Fil- 
ling factor essentially depends on lattice type 
and must be calculated for each specific lattice. 

The fraction of volume occupied by white spheres 
for a<1 is found by multiplying the filling 
factor by the fraction of while spheres, that is, 
by x. The fraction of volume filled up by white 
spheres is thus fz. At percolation threshold it 
equals fxs. If the hypothesis on the universality 
of the fraction of volume at which percolation 
sels in is correct, the product fr, must be identi- 
cal for all lattices. 

The filling factors for a number of lattices are 
listed in the second column of Table 4. In order 


ge 
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Table 4 
The Product fzs for Different Lattices 
Lattice type Í Xs tt, 

Plane lattices 
Square 0.79 | 0.59 | 0.47 
Triangular 0.91 | 0.5 0.46 
Honeycomb 0.61 | 0.7 | 0.43 

3D lattices 
Simple cubic 0.52 | 0.31 | 0.16 
Body-centered 0.68 | 0.25 | 0.17 
Face-centered 0.74 | 0.20 | 0.15 
Diamond 0.34 | 0.43 | 0.15 


to illustrate how they were calculated, let us de- 
termine f for the honeycomb lattice given in 
Fig. 34. It was shown in Exercise 4 to Chapter 5 
that the area per site of the honeycomb lattice is 
(3V 3/4) a2, where a is the length of the side of 
the unit cell. This result has the following mean- 
ing: let us draw on a plane, on which the lattice 
was constructed, a square, a rectangle, a circle, 
or any other geometric figure, but necessarily 
such that its size be much greater than the distance 
between the neighboring sites of the lattice. 
Divide its area by the number of lattice sites that 
fell within this figure. The area per site is the 
limit to which this ratio tends when the figures 
become infinitely large. 

The fraction of area occupied by circles equals 
the limit of the ratio of the area occupied by the 
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circles to the area of the large figure. The area 
occupied by the circles equals the product of the 
number of sites belonging to the large figure by 
the area of a single circle. In other words, the 
quantity f equals the ratio of the area of one cir- 
cle to the area per site. 

The circles drawn in Fig. 31 have the radius 
a/2, and hence, their area is na?/4. This gives 


_ na 33a n a 
f= 7 a Tea 


The filling factors for other lattices are calculated 
in a similar manner; you can see from Table 4 
that f varies within a wide range. 

The last column of Table 4 lists the products 
fxs. You can see that the hypothesis that fa, is 
independent of lattice type does not hold too 
well. However, its variations are small both 
within the group of plane lattices and within the 
group of 3D lattices. Hence, to within 10 to 15% 
the following formulas hold: 


fiey 0.5 (3) 


for plane lattices, and 
fr, = 0.16 (4) 
for 3D lattices. 

The calculation of the filling factor f being rel- 
atively simple, formulas (3) and (4) make it 
possible to evaluate percolation threshold of the 
site problem for arbitrary lattices. 

It is readily understood that the critical fraction 
of volume occupied by white spheres, at which 
percolation appears, decreases monotonically 
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with increasing dimensionality of space. In one- 
dimensional space, that is, in a linear string of 
sites, percolation through white sites vanishes 
at no matter how small concentration of black 
sites. Even a single black site blocks the perco- 
lation path because no bypass is possible. ‘The 
possibility to go around black sites appears in 
plane (two-dimensional) lattices, and such pos- 
sibilities become more numerous in 3D lattices 
because detours are nol restricted to a single 
plane. 

The critical volume concept proves fruitful not 
only in lattice problems. In Chapter 9 we shall 
encounter a problem in which white and black 
balls are not located at lattice sites but are ran- 
domly poured into a jar. We shall be interested 
in percolation through white balls in contact with 
one another. This percolation is also found to arise 
when the volume occupied by white balls comes 
to about 0.16 of the total volume. This result 
changes only slightly if the balls have unequal 
radii. 

Chapter 10 treats the problem of space painted 
in a random manner by white and black paints. 
It is found that percolation through regions of one 
of these colors appears in the 2D case when the 
fraction of surface area painted white (or black) 
is exactly equal to 0.5, and in the 3D case when 


the fraction of volume painted white (or black) 
1S approximately equal to 0.46. 


Exercise 


1. Check whether 


, the filling factors f listed 
in Table 4 were corre 


ctly calculated. 
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Chapter 7 


Ferromagnetics 
with Long-Range Interaction. 
Phe Sphere Problem 


So far, when considering site and bond problems, 
we have always assumed that each site can be di- 
rectly connected only with its nearest neighbors, 
and that the connections between farther removed 
sites are established through strings of sites each 
of which is connected with its nearest neighbors. 
In this chapter the site problem is generalized to 
the case in which the sites that are not nearest 
neighbors are directly connected with each other. 
This problem may be of practical importance, and 
so it may prove useful to know that it has been 
sufficiently well elaborated. 

If the number of sites with which a given sile 
is connected becomes very large, the site problem 
transforms to a quite new problem referred to as 
the sphere problem. In fact, this problem plays an 
especially important role in percolation theory: 
it is used to understand the transition to metallic 
conductivity taking place in semiconductors as 
the impurity concentration in them increases; its 
solution is at the basis of the theory of hopping 
conduction in semiconductors, which is an impor- 
tant and interesting phenomenon observed at very 
low temperatures. Consequently, numerous scien- 
tists worked with the sphere problem, and intere- 
sting results have been obtained for this and 
similar problems. 

The sphere problem is also interesting because 
this is the first nonlattice problem that we en- 
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counter in this book. The random elements with 
which it operates are defined not at sites of a 
periodic lattice. 


Ferromagnetics with Long-Range Interaction 


Let us return to the problem of a ferromagnetic 
with nonmagnetic impurity atoms, discussed in 
Chapter 3. There it was assumed that magnetic 
atoms orient their magnetic moments parallel to 
each other only if they are nearest neighbors or 
if they are connected through a string of magnet- 
ic atoms that form Successive pairs of nearest 
neighbors. But if all nearest neighbors of a mag- 
netic atom were nonmagnetic atoms, the magnet- 
ic moment of this atom was assumed to be arbi- 
trarily oriented. 

This model was based on the fact that the inter- 
action between magnetic moments resulting in 
the parallel arrangement decreases very rapidly 
with distance; in fact, the decrease is so steep 
that the magnetic moments that are not nearest 
neighbors simply “ignore” one another, that is, do 
not interact at all. 

In crystallography, the group of all nearest- 
neighbor atoms of a given atom is usually called 
ils first coordination group, and it has already 
been mentioned that the number z of nearest 
neighbors is called the coordination number, 
The ensemble of equivalent atoms that are 
next-lo-nearest neighbors of the given atom 
is called the second coordination group, and 
so on. As an example, let us consider the 
simple cubic lattice (see Fig. 12). The first coor- 
dination group in this lattice Consists of six atoms 
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located at the edges of the cube projecting from 
the initial atom. The second coordination group is 
formed by 12 atoms located at the diagonals of the 
cube’s faces passing through the initial atom. 
And finally, the third coordination group consists 
of eight atoms located at the body diagonals of 
the cube passing through the initial atom. 

In Chapter 3 we assumed that the interaction 
between atomic magnetic moments vanishes be- 
yond the first coordination group. Under this 
assumption, the calculation of the critical con- 
centration of magnetic atoms at which spontane- 
ous magnetization appeared (or vanished) re- 
duced to the site problem with bonds only between 
nearest neighbors. 

The assumption of short-range interaction does 
not always hold, so that it is justifiable to discuss 
a problem in which the interaction between mag- 
netic moments covers several coordination groups, 
and to determine the critical fraction of magnetic 
atoms resulting in spontaneous magnetization. 

This problem reduces to the site problem in 
which bonds are formed not only between nea- 
rest neighbors. In fact, its formulation contains 
nothing new. Sites may be white and black (mag- 
netic and nonmagnetic). Two white sites are re- 
garded as connected if the bonds connect the 
coordination groups in which they are located. 
If site A is connected with site B, and site B with 
site C, then A is connected with C. The ensemble 
of connected sites forms a cluster. The percola- 
tion threshold is defined as the fraction of white 
sites at which an infinite cluster is born. 

When bonds are made to extend to farther coor- 
dination groups, the percolation threshold z must 
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Table 5 


Percolation Threshold in the Site Problem with 
Bonded Nonnearest Neighbors* 


Lattice type Z Xg 


Plane lattices 
Honeycomb, 1 3 0.700 


2.10 
Square, 4 4 0.590 | 2.36 
Triangular, 1 6 0.500 | 3.00 
Square, 1, 2 8 0.410 | 3.28 
Triangular, 1, 2 12 0.295 | 3.54 
Honeycomb, 1, 2, 3 12 0.300 | 3.60 
Square, 1, 2, 3 12 0.292 | 3.50 
Triangular, 1, 2, 3 18 0.225 | 4.05 


3D lattices 


Diamond 4 0.425 | 1.70 
Simple cubic, 1 6 0.307 | 1.84 
Body-centered, 4 8 0.243 | 1.94 
Face-centered, 4 12 | 0.1495 | 2.34 
Body-centered, 1, 2 14 0.475 

Simple cubic, 1, 2 18 0.137 | 2.47 
Face-centered, 1, 2 18 0.136 | 2.45 
Simple cubic, 4, 2,3 26 0.097 | 2.52 
Body-centered, 1, 2, 3 26 0.095 | 2.47 
Face-centered, 4, 2, 3 42 0.061 | 2.56 


— 


* The listed thresholds for a single coordination group are 
not always exactly equal to the data given in the preceding 
tables. The reason for this is that those tables gave the data 
that we considered the m s 


ost reliable, while the present table 
compares the data calculated by the same method. 


E a U 
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obviously diminish. The greater the number of 
bonds originating at a given white site, the great- 
er the probability that at least one of these 
bonds will lead to another while site, and cor- 
respondingly, the smaller the number of white 
siles necessary to ensure percolation. 

Clearly, this problem is nota bit simpler (rath- 
er, it will be more complex) than the standard 
sile problem. However, a number of problems of 
this type was solved by various approximate 
techniques, and the results of one of them are 
listed in Table 5. The first column of the table 
gives the lattice type and the number of coordi- 
nation groups to which the bonds were extended. 
The second column shows the number of sites, 
Z, with which each site is connected, that is, the 
total number of sites located in the considered 
coordination groups (in the case of a single coor- 
dination group it coincides with the coordination 
number 2). 

The last column of the table gives the product 
Zx;. This product, as we indicated in the preced- 
ing section, is strongly dependent on lattice 
type in the case of the site problem with bonds 
only in the first coordination group. However, as 
we find in Table 5, at high values of Z it changes 
the slower, the greater Z is. This is especially 
clear for three-dimensional lattices where large 
values of Z are used. Obviously, the product Ze, 
tends to a number of the order of 2.6-2.7 inde- 
pendent of lattice type. 

In the 3D case the number 2.7 is now consid- 
ered the most likely value (to within 0.1) of 
the quantity Be defined as the limit of Zz, for 


140 Part II. Various Problems of Percolation Theory 


very large Z: 


B. = lim Zz, (1) 
Z-r00 

It will be explained in the subsequent sections 
why this limit exists and why it is independent 
of lattice type, being only a function of the di- 
mensionality of space, that is, of whether the lat- 
tice under consideration is plane or three-dimen- 
sional. 

In order to understand this feature, it is neces- 
sary to study the so-called sphere problem. 


Exercise 


1. Find the location of 42 sites belonging to the 
first three coordination groups of an fee lattice, 


The Sphere (Circle) Problem 


Now we shall consider a different problem or, 
rather, a problem that seems to be different at 
first glance. Assume that a plane is covered with 
circles of identical radius R whose centers are 
spread on the plane quite randomly and, on the 
average, uniformly. This means that both coor- 
dinates of centers of circles are random numbers 
distributed uniformly in the range from zero to 
L, where L is a very large length (in comparison 
with R) characterizing the size of the system under 
discussion. The important distinctive feature of 
this problem is that there is no limitation on 
the extent to which the circles overlap. The 


mean number of centers of circles per unit area is 
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N. In other words, N is the concentration of 
centers of circles. 

Two circles are said to be mutually connected 
if the center of one of them lies within the other; 
such circles are sometimes said to be inclusive. 
If circle A is connected with circle B, and B is 


Fig. 32. Percolation paths through inclusive circles are 
traced by the broken lines. The dots mark the centers of 
circles. 


connected with C, then A is connected with C. 
Therefore, too far removed circles can be connect- 
ed via a string of inclusive circles (Fig. 32). 
The problem now is to find the critical concen- 
tration V at which percolation sets in through 
inclusive circles, that is, paths are formed tra- 
versing the whole system and composed of inclu- 
sive circles. (In other words, an infinite cluster is 
formed out of mutually connected circles.) 
The problem involves two parameters: the con- 
centration M and radius R. (We have left out 
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the size of the system, L, but it is clear that if 
the system is sufficiently large, the critical value 
N is almost independent of L.) It is not difficult, 
however, to ascertain that the presence or absence 
of percolation depends not on two parameters bul 
only on one, namely, on the dimensionless pro- 
duct NR. (The dimensionality of concentration 
in a plane problem is cm-®.) Let us choose for 
this parameter the mean number of centers of 
circles within a circle. It equals 


B = nNR? 


The following arguments will immediately 
show that percolation sets in at a certain value 
of the parameter B, regardless of the specific val- 
ues of V and R making up B. Imagine looking at 
a plane with the circles drawn on it. Let us mag- 
nify this picture severalfold, for instance, through 
a projector. This will be a transformation vary- 
ing NV and R but leaving B invariant because the 
mean number of centers of circles within a circle 
will not be changed by this magnification. 

It is also easy to understand that this tr 
mation does not affect percolation. If there was 
no percolation in the initial picture, there will 
not be any in the magnified picture, and con- 
versely, if there was percolation through inclusive 
circles in the initial picture, this percolation will 
not vanish because of magnification. 

The transformation changing N and R but not 
B thus does not affect percolation. Consequently, 
it depends only on the value of 


the parameter 2 
whether there will be or not be percolation in a 


system. Percolation is present when B is large, 
and absent when it is small. 


anslor- 
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This new problem that has just been formulat- 
ed above is referred to as the circle problem. Its 
three-dimensional analogue is called the sphere 
problem. The sphere problem is formulated as 
follows. The coordinates of centers of spheres of 
radius R are generated in the three-dimensional 
space by a suitable random-number generator. 
Two spheres are said to be connected (or inclu- 
sive) if the center of one sphere lies within another 
sphere. We need to determine the critical con- 
centration of centers at which percolation sets in 
through inclusive spheres. 

It will be readily understood that, as in the 
plane problem, the presence or absence of perco- 
lation is delermined exclusively by the value 
of the parameter B which is the mean number of 
centers of spheres within a sphere: 


B = (4/3) 2NR? 


where V is the mean number of centers of spheres 
per unit volume (the dimensionality of volume 
concentration is cm~). 

We have already mentioned that the sphere 
problem has important applications to the theory 
of electric conduction in semiconductors at low 
temperatures. For this reason, it was analyzed 
by a number of authors resorting to most various 
techniques. By today’s data, the critical value Ba 
at which percolation through spheres sets in is 
2.7 + 0.4. The circle problem was not studied as 
thoroughly, so that a substantial spread is found 
among the results published by different authors. 
It appears that Be = 4.1 + 0.4. 
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The Circle (Sphere) Problem 
Is the Limiting Case of the Site Problem 


Let us return to the site problem in which not only 
nearest neighbors are connected, and explain 
why the limit on the right-hand side of formu- 
la (1) exists, why it is independent of lattice 
type, and finally, why we denote it by the same 
symbol Be as the threshold values of the circle 
and sphere problems. 

We begin with slightly reformulating the site 
problem. To be specific, we first discuss plane 
lattices. Generalization to 3D lattices proves to 
be quite simple. 

-et us surround each while site with 
radius R that was chosen to be greater than the 
distance from this site to the sites of the farthest 
coordination group with which the given site is 
connected, but smaller than th 
sites of the next coordin 
that two circles are connected if the white sites 
around which tl 


This means that two circles 


a circle of 


are connected if the 
tin the other, that is, 


an infinite cluster of connect- 
quivalent to the onset of per- 
clusive circles. In other words, 
gh the circles of radius R drawn 
around white sites sets in at a critical fraction 
of white sites, x, (Fig. 33). 

The meaning of the quantity Zx is easily un- 
derstood. The quantity Z is the number of sites 
(both black and white) that lie within a circle. 
The product Zz is the mean nuniber of centers of 
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other circles that lie within a circle (or the mean 
number of white sites within a circle). The quan- 
lily Zr, is the mean number of centers of circles 
within a circle at which percolation sets in, 
that is, at which infinite paths can be traced 
through inclusive circles. 

It is now obvious that the product Zz, has the 
same meaning as the quantity Be in the circle 


1 SEEE T E k aN: 
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Fig. 33. Percolation path through inclusive circles on a 
square lattice. The interaction is taken into account at 
a distance of up to three distances to the nearest neigh- 
bors. The percolation path is traced by the broken line. 


problem. Quite likely, now it may be more dif- 
ficult for the reader to point out the difference 
between the site and circle problems than to no- 
tice the similarity. But the difference is there, 
and it is an important one. The fact is that in 
the circle problem the points that can serve as the 
centers of circles can lie anywhere on the plane, 
while in the site problem these are only the sites 


10 —0240 
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of the lattice in question (see Fig. 33). If the 
number of sites within a circle is small, the differ- 
ence between the two problems becomes very 
important. It is then natural that the critical 
value of Zr; is a function of lattice Lype. The to- 
tal number of sites lying within a circle is Z. As 
we see from Table 5, Zz, indeed varies among 
lattices when Z is not too large. 

However, if Z is large, the difference between 
the two problems fades away. Imagine that we 
have started with the circle problem, but then 
shifted the center of each circle to the lattice 
site nearest to it. This is already a site problem. 
If the number of lattice sites within a circle is 
very large, this shift will not, with a high prob- 
ability, connect the circles that were not con- 
nected before the shift, and vice versa. This argu- 
ment shows that the site and circle problems 
become equivalent for Z — œ. The value of B, 
determined by formula (1) is independent of lat- 
tice type and coincides with the quantity Be de- 
fined in the circle problem. 

This chain of argumentis can be carried over 
entirely to the 3D case. The quantity Be defined 
by formula (1) coincides for 3D lattices with Be 
defined in the sphere problem. 

To recapitulate, the site problem on any plane 
lattice is reducible at large Z to the circle prob- 
lem, and the site problem on a 3D lattice is re- 
ducible to the sphere problem. The limit in formu- 
la (1) is, therefore, indeed independent of lattice 
type but depends on the dimensionality of the 
space on which the problem has been formulated. 


(The values of B, are different for the circle and 
sphere problems!) 
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Chapter 8 


Electric Conduction 
of Impurity Semiconductors. 
The Sphere Problem 


Percolation theory proved extremely useful for 
the understanding of the processes in doped se- 
miconductors (called extrinsic semiconductors to 
distinguish them from pure ones, called intrinsic 
semiconductors). It can rightly be said that at the 
present time impurity semiconductors represent 
one of the main fields to which percolation theory 
is applied. Several chapters of this book deal 
with impurity semiconductors. The present chap- 
ter begins with an exposition of current notions 
of semiconductor science, briefly interrupting the 
story of percolation theory. 


Intrinsic Semiconductors 


Let us start to discuss pure semiconductors, choos- 
ing as a basis semiconducting elements of the 
fourth group of the periodic table, such as ger- 
manium and silicon. These elements crystallize 
into diamond lattice (see Fig. 30a, b). The four 
electrons making up the outer shell of each atom 
form bonds with the four nearest neighbors. The 
electrons move in such a way that their density 
is distributed nonuniformly around the atoms, 
being concentrated in four “strands” stretching 
from the centers of tetrahedrons, in which the 
atoms are located, to the vertices of these tetra- 
hedrons. It is these strands that bind up the atoms 


10* 
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of the crystal and do not let them fly apart. 

It can be imagined, by way of illustration, thal 
each atom has four hands and holds hands with 
its four neighbors. Each eleetron is then strongly 
hound lo a neighbor, so that a moderately strong 
electric field applied to the semiconductor can- 
nol generate an electric current: the number of 
electrons is exactly equal to that necessary lo 
form the bonds, while the energy required to re- 
lease a bonding clectron is quite high. This 
energy, or rather, the minimum work that has to 
be done to transfer an electron from the bound 
stale to the state in which it can move freely 
through the crystal is called the energy gap width 
(or simply gap width) and is denoted by Bi 

The gap width is an extremely important char- 
acteristic of a semiconductor, determining to a 
large extent all its electrical properties. Mor in- 
stance, let us consider electric conduction. The 
electric conduction in pure (intrinsic) semiconduc- 
tors exists because a certain f raction of bonds 
are broken. A bond breakdown releases an elec- 
tron ready to carry electric current and a “crip- 
pled”, “three-armed” atom. This “cripple” is also 
a charge carrier, only its ign is opposite to that 
of the electron, Indeed, an clectric field can 
transfer an electron from a neighbor to the “three- 
armed” atom, changing its energy only slightly. 
As a result, another atom becomes “three-armed”. 
This process will continue, and you easily un- 
derstand that in contrast to electrons moving in 
the field, say, from left to right, a “three-armed” 
atom moves from right to left. (It is important 
to realize that the atoms themselves do not 
move, In fact, a moving electron turns the atom 
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which it abandons into a “three-armed” atom.) For 
this reason, a “three-armed” atom can also be 
regarded as a charge carrier bul, in contre st to the 
electron, it must be assigned the positive, not 
negative, charge. 

_ The microscopic picture is somewhat different 
in other semiconductors consisting of atoms with 
a different valence, but in all cases the breakdown 
of a bond produces two charge carriers with oppo- 
site signs. One of them (the negative) is called 
the electron, and the other (the positive) is cal- 
led the hole. A “three-armed” atom is thus a par- 
ticular case of hole. 

The energy required to brea 
supplied from the energy of thermal motion of 
the atoms. As it is proved in statist cal physics, 
at high temperatures the mean energy of thermal 
motion per atom that performs sma l-amplitude 
vibrations around the equilibrium position is 
BET, where T is temperature in K, and & is the 
Boltzmann constant ( = 1.388-10- erg/K). 

As a rule, the gap width Æg is measured in elec- 
tron volts (eV). One electron volt is the work 
done by an electron traversing a potential differ- 
once of 1V: 1 eV =1.6-10- J = 1.6-10-" erg 
In germanium Eg = 0.7 eV, and in silicon £y 
— 14 eV. At room temperature (300 K) the 
energy 3kT' is only 0.08 eV. This is much less 
than necessary for a bond breakdown both in 
germanium and silicon. 

However, thermal motion is chaotic, and at 
some random moments the energy of motion of an 
atom may be very high. This happens extremely 
rarely to cach individual atom, but the total 
number of atoms is very large. There are roughly 


k down a bond is 
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10% atoms in 1 cm3. (The reader may remember 
that the number of atoms per mole equals Avo- 
gadro’s number, 6-105.) f 

For this reason, the concentration of electrons 
and holes released by broken bonds is much smal- 
ler than the concentration of atoms; nevertheless, 
it is so high that it can result in an appreciable 
electrice conductivity in a certain temperature 
range. Calculations show that in germanium the 
concentration of electrons is about 101 em-3 at 
T = 300 K. In silicon, where the gap is wider, 
the concentration at the same temperature is 
much lower (101° cm~*). As temperature is red uced, 
the concentration of electrons decreases sharply, 
and the resistance of an intrinsic semiconductor 
undergoes a corresponding sharp enhancement. 

The number of charge carriers at a given tem- 
perature is the smaller, the wider the energy gap 
is. The difference between a semiconductor and 
a dielectric is merely the width of the gap. Asa 


rule, materials with Eg of about 5 eV or more are 
referred lo as die 


lectries. They contain practical- 
ly no charge carriers at room temperature, 
While the difference between diclectrics and 
semiconductors is more quantitative than quali- 
tative, the difference between metals and diclec- 
trics is a principal one. The energy gap of a me- 
tal is zero, and the concentration of charge car- 
riers is high even in the immediate vicinity of 
the absolute zero of the temperature scale, 


Impurity Semiconductors 


Let us assume that germanium or silicon is doped 
with an impurity element from the fifth group of 
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the periodic table, such as phosphorus, antimony, 
or arsenic. An atom of these elements has five 
electrons on the outer shell. If such an atom sub- 
slilutes, for example, a germanium atom located 
at the center of a tetrahedron, four of its elec- 
trons go to form the bonds with theyfour neigh- 
bors and the fifth electron remains dangling. 
Chis electron is bound to its alom because if the 
electron moves away, the atom becomes positive- 
ly charged and attracts the electron back. There- 
fore, the electron can be removed to a considerable 
distance away from the atom only by doing a work 
against attractive forces. This work is what we call 
the energy of bonding of the electron to an atom. 

As we will show below, the bonding energy of 
this extra electron is relatively low, so that the 
distance between this electron and the impurity 
atom is rather large (in comparison with the period 
of the lattice). As a result, the structure of the 
impurity atom resembles that of the simplest 
among atoms, namely, the hydrogen atom. We 
remind the reader that the hydrogen atom con- 
sists of a positively charged heavy nucleus and a 
light negative electron, the size of the nucleus 
being negligibly small in comparison with the 
distance from the nucleus to the electron. 

An impurity alom has a similar structure. The 
unit that acts as the nucleus of a hydrogen atom 
is not just the nucleus of the impurity atom but 
includes the inner-shell electrons and the four 
electrons that form the bonds. We shall be able 
to show now that the size of this complex is small 
compared with the distance to the extra elec- 
tron; the charge of the complex is positive and 
equals in magnitude the electron charge. 
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The electron in a hydrogen atom is known to be 
at a distance of the Bohr radius ay from the 
nucleus: 

B , 240-8 { 
an = sy & 0.53-40-8 cm (1) 
Here A = 1.05-10-27 erg-s is Planck’s constant 
(divided by 2x), m = 9.8.19-28 g is the electron 
mass, and e = 4.8-10- CGSB units is the elec- 
tron charge. 

The fundamental constant & was introduced hy 
the German physicist. Max Planck in 1904 in 
connection with the hypothesis on the quantum 
nature of electromagnetic waves. This constant 
enters into the quantum-mechanical equation 
describing the motion of an electron around the 
nucleus. The only other constants in this equation 
are e and m. The Bohr radius ay is the only pos- 
sible quantity with the dimensionality of length 
that can be Composed of e, m, and hi. 

The bonding energy of the electron in the hyd- 
rogen alom equals 


Ey = a 13.6 eV (2) 


The following interpretation of this formula by 
the order of magnitude is possible. The positively 
charged nucleus creates a potential e/rat a dis- 
tance r. At a distance of order ap it equals elap, 
and at infinity it is zero. Consequently, the work 
that has to be done to transfer an electron located 
at a distance of order ag from the nucleus lo an 
infinilely Jaree distance equals, by the order of 
magnitude, e?/ap (the potential difference times 
the electron charge), 
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Let us return to the impurity atom. The for- 
mulas written for the hydrogen atom have to be 
modified in order to take into account that the 
extra electron moves not in a vacuum but in a 
semiconductor crystal. This modifies the form of 
Coulomb’s law. The force acting on an electron 
separated by a distance r from the nucleus is 
now e/er? (in CGSE units), where e is the dielec- 
trie constant, or permittivity, of the crystal, 
and not e?/r?. This complication can be taken 
into account by replacing the quantity e by 
ele in the expressions for ag and Æg. Further- 
more, it should be borne in mind that the mass 
describing the motion of an electron through the 
crystal does not equal that of a free electron, m. 

The point is that the nuclei of the semiconduc- 
tor’s atoms and the inner-shell electrons produce 
in the crystal a periodic electric potential. One 
of the most interesting conclusions of the quantum 
theory of solid state is that the electron very 
nearly “overlooks” this potential if it is exactly 
periodic. The qualification “very nearly” consists 
in that the neglection of the periodic potential 
in the equations of motion of the electron must be 
simultaneously accompanied by replacing the elec- 
tron mass m by a mass m* which is a function 
of crystal properties. The quantity m* is called 
the effective mass. 

An impurity atom produces a nonperiodic po- 
lential that can by no means be omitted. But, 
when the motion of the electron in this potential 
is described, we can neglect the periodic potential 
of the erystal provided we replace m by m*. 

The equation describing the motion of the extra 
electron in a crystal around a charged impurity 
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atom thus differs from the equation of motion of 
an electron in the hydrogen atom by the re- 
placements è — e/s and m —> m*. . 

Let us denote the characteristic distance to 
which the extra electron is removed from the im- 
purity atom by aj, and the bonding energy of 
this electron by £$. Making use of formulas (1) 
and (2), we obtain 


ay, = 0.53.1078 (5) e [em] (3) 
À 1 
B$= 13.6 (2$) 4 lev] (4) 


As a rule, the effective masses in semiconduc- 
tors are substantially smaller than the mass of 
the free electron, and the dielectric constant 
equals 10-20 (e.g. in germanium m* œ 0.14 m 
and e = 416). Consequently, the characteristic 
distance at which the extra electron is found 
equals, in typical semiconductors, from several 
tens to several hundreds of angstroms (1 À = 
= 10-8 em), which is much greater than the in- 
leratomic distance (e.g. in germanium this spa- 
cing is 2.45 A). 

We have thus obtained that the e 
in an impurity atom is removed to a considerable 
distance from the atom and is kept at this distance 
by the attractive forces caused by the nueleus. 
The fact that the distance aj, is large is already 
a sufficient indication of smallness of the bonding 
energy “ii. Indeed, the data cited above show 
that in typical semiconductors the bonding energy 
E% is from several hundredths to several thous- 
andths of an electron volt; hence, it is much 


xtra electron 
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smaller than the gap width Æg (e.g. in germanium 
Eg = 0.7 eV). 

This is natural because it is easier to break 
loose the extra electron off the impurity atom 
than the electron bonding the host atoms. The 
thermal energy kT reaches the level of 0.01 eV 
at a temperature of the order of 100K. As 
a rule, at this temperature a considerable frac- 
tion of extra electrons separate from their impu- 
rity atoms and take part in the transfer of elec- 
tric current. Therefore, impurity atoms belong- 
ing to the fifth group of the periodic table 
let go off their extra electrons rather easily. 
For this reason, they were given the name 
donors. 

Suppose that the impurity atoms belong to an 
clement of the third group of the periodic table, 
for example, boron, aluminium, gallium, or others. 
These atoms have three electrons on the outer 
shell, so that they lack one electron for forming 
the bonds with their four neighbors. This elec- 
tron is readily borrowed from the neighboring 
atoms of the semiconductor, but then one of the 
neighbors becomes “three-armed” or, in other 
words, a hole appears in the neighborhood of the 
impurity atom. The impurity atom captures a 
fourth electron and thereby becomes negatively 
charged. The hole is thus attracted to this atom 
by electric forces, and a work has to be done to 
pull the hole away. This work is called the bond- 
ing energy between the hole and the impurity 
atom. 

The calculation of the bonding energy of the 
hole again leads to the problem of a hydrogen-like 
atom, but in this last case it is a positively charged 
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hole that moves around a negatively charged 
nucleus. The role of the nucleus is now played by 
an impurity atom together with the captured ad- 
ditional electron. This additional electron turns 
into a strand forming an interatomic bond, so 
that the region in which it moves does not exceed 
one interatomic distance. At the same time, a 
hole is bonded much less strongly to an impu- 
rity atom. A distance from the atom to the hole, 
and the bonding energy, are given by formulas (3) 
and (4). We only have to take into account that 
the effective mass of the hole that must be used 
in this case differs, generally speaking, from the 
electron effective mass. As a rule, this mass is 
also considerably smaller than the free electron 
mass, so that a hydrogen-like atom with a hole 
also has a size of the order of tens of angstroms, 
and its bonding energy is of the order of several 
hundredths of one electron volt. 

At temperatures of about 100 K the thermal 
motion breaks the bonds between holes and 
impurity atoms, after which the holes begin an 
“independent way of life” and, when an electric 
field is applied, participate in the electric cur- 
rent. 

Impurity atoms from the third group of the 
periodic table thus readily accept an electron 
and form a hole; hence, the term acceptors for 
such dopants, 

Let us summarize this section. Mobile charge 
carriers are formed in semiconductors only at the 
expense of the thermal motion energy. They can 
be formed when bonds of the lattice are broken. 
This requires that a work be done equal to the 
energy gap width Eg. The outcome is the simul- 
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taneous appearance of an electron and a hole. An 
electron or a hole can also be born individually 
when their bond to an impurity atom is broken. 
The bonding energy between an electron or a hole 
and an impurity atom being much smaller than 
the gap width, the probability for the electron 
to break away from the impurity alom is much 
higher than the probability of breaking a lattice 
bond. On the other hand, the concentration of 
impurity atoms is usually lower by many orders 
of magnitude than the concentration of lattice 
atoms. Consequently, when temperature is in- 
creased, first of all charge carriers on impurity 
atoms are separated, the concentration of intrin- 
sic carriers (electrons and holes born of broken 
lattice bonds) being negligibly small. This tem- 
perature range is called the range of extrinsic 
conduction. However, as the temperature is 
further increased, the concentration of intrinsic 
charge carriers becomes comparable with that of 
impurity atoms and then grows larger than it. 
This range is called the range of intrinsic conduc- 
tion. 

The conclusion which is the most important for 
further exposition is that at a very low tempera- 
ture at which the thermal motion energy kT is 
small in comparison with the bonding energy of 
electrons to impurity atoms, the semiconductor 
has neither extrinsic nor intrinsic charge carriers. 
The bonds are intact where they ought to be, and 
impurity electrons and holes are localized in the 
vicinity of their atoms. In this range, electric 
conduction of the semiconductor drops very 
sharply to zero as its temperature is lowered. 
This behavior distinguishes a semiconductor 
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from a metal in which the concentration of mobile 
charge carriers remains high no matter how 
small the temperature is. 


Transition to Metallic Electric Conduction 
at Increased Impurity Concentrations 


This principal difference between semiconductors 
and metals suddenly disappears when impurity 
concentration is increased. If concentration goes 
beyond a certain critical value Ne, electric 
conductivity remains relatively high and weakly 
dependent on temperature no matter how low the 
temperature is. This electric conduction is refer- 
red to as metallic. This does not mean that the 
electric conductivity of a semiconductor becomes 
comparable with that of good metals. Far from 
it: the electric conductivity of a semiconductor 
always stays many orders of magnitude lower, 
The term only reflects the behavior of electric 
conductivity at low temperatures. The transition 
to metallic electric conduction that takes place at 
an increased impurity concentration is called the 
metal-insulator transition, or the Mott transition 
(after the famous British physicist Sir Nevill 
Mott). 

Experiments have demonstrated that the criti- 
cal impurity concentration Ne at which the Mott 
transition takes place varies quite considerably 
among different semiconductors. A reliable esti- 
mate of the quantity Ne is given by a relation 


Nean ~ 0.02 (5) 


For instance, in germanium Ne = 1017, and in 
Silicon M, = 3-1018 cm, 


on 
D 
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Even now we lack a satisfactory mathematical 
theory of the Mott transition in semiconductors. 
In fact, this is one of the most complicated prob- 
lems in solid state theory. The qualitative picture 
looks as follows. 

The structure of impurity atoms resembles that 
of univalent elements of the first group of the 
periodic table (H, Li, Na, K). Like these ele- 
ments, an impurity atom has a single electron on 
its outer shell. When crystallized, elements of 
the first group form good metals. Hydrogen forms 
a molecular crystal which is an insulator. How- 
ever, there is every reason to believe that this 
insulator also converts to a metal at a sufficient- 
ly high pressure. (The research in this field is 
conducted nowadays with such vigor that it could 
make the topic for a separate book.) 

It appears as nearly obvious that if hydrogen- 
like impurity atoms are distributed in a semi- 
conductor with a sufficiently high density, they 
will also form a metallic system. 

It might seem strange at first glance that this 
argument has a bearing on the Mott transition. 
Indeed, the concentration of atoms in a metallic 
sodium is of the order of 107? cm™, that is, great- 
er than the concentration of impurities in ger- 
manium, at which the Mott transition occurs, 
by a factor of about 100 000. Obviously, sodium 
cannot be a metal at such low concentra- 
tions. 

However, it is not just the concentration but 
the degree to which the electron shells of neigh- 
boring atoms overlap that determines whether a 
system is metallic or not. If the atoms are at 
such a large distance that the regions of space in 
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which their outer-shell electrons move are tar 
from one another, these are simply isolated atoms. 
But if the nuclei of the atoms are SULTS 
close for these regions to overlap, the atoms lose 
their individuality. Their outer-shell electrons 
“cannot know” to which of the nuclei they be- 
long. A system of a large number of atoms Appin- 
priates, or collectivizes, electrons. These electrons 
form a separate system capable ol conducting an 
electric current. This material is a metal. 

The degree of overlapping is determined for 
hydrogen-like atoms by a dimensionless parame- 
ter Nai’. (The quantity (4/3) naj" is the mean 
number of atomic nuclei within a sphere of radius 
equal to the effective Bohr radius.) It has already 
been mentioned that since an impurity atom is 
located inside a semiconductor, its radius is ano- 
malously large. For this reason, condition (5) is 
met already at N = 1017 cm- (for germanium). 

It must be clear by now from the arguments 
given above why the Mott transition oecurs at ap- 
proximately the same value of the parameter 
Nai in all conductors, although critical con- 
centrations N, may vary quite considerably. 
Indeed, it is this parameter that determines the 
overlapping of neighboring atoms! 

The following question can be posed: What 
should be the concentration of atoms in a hydrogen 
crystal for the overlapping of its atoms to be the 
same as the overlapping of impurity atoms in a 
semiconductor at the concentration corresponding 
to the Mott transition? The overlapping is given 
by formula (5), but now it is natural to substitute 
into it the Bohr radius found from formula (1). 

This will give an enormously high concentration: 


aw 


wv 
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10% em™. Hence, the overlapping corresponding 
to the Mott transition is extremely large. 
There is no ground, therefore, to be surprised 


semiconductors with a high impurity concen- 
tration. 


The Mott Transition and Sphere Problem 


Percolation theory offers a simplified description 
of the Mott transition which takes into account 
that impurities are located in a semiconductor 
randomly and can form tighter or looser ensem- 
bles. Assume that impurity atoms in a certain re- 
gion are so close to one another that their electron 
shells strongly overlap and the outer-shell elec- 
trons become “collective property”. This region 
is a piece of metal: a potential difference applied 
to this region would cause an electric cur- 
rent. 

The presence of such regions is not sufficient, 
however, for a large sample to behave as a me- 
tal. If metallic regions are infrequent, they are 
not in contact, forming only isolated metallic 
islands in the material that at low temperatures 
behaves as a dielectric. Taken as a whole, the 
combination is equally dielectric. 

As impurity concentration increases, the frac- 
tion of space occupied by metallic regions grows, 
and at a certain critical concentration N, the 
metallic regions form a connected system of “lakes 
and canals” that permeates the whole semiconduc- 
tor crystal. Beginning with this concentration, 


11-0240 


162 Part II. Various Problems of Percolation Theory 


the electric conduction of a large sample be- 
comes metallic. 

Obviously, the concepts presented above must 
be formulated mathematically in terms of perco- 
lation theory, although the formulation will nol 
come easily. The main difficully is that we do 
not know the impurily concentration at which a 
region can be regarded as metallic. ; 

The simplest model of the Mott transition 
enunciated in the sixties reads as follows. Ima- 
gine that each impurity atom is a metallic ball 


(a) th) 
Fig. 34. (a) Overlapping spheres representing aloms with 
collectivized electrons; (b) a string of overlapping spheres 


forming a metallic channel through which electric current 
can flow. 


of radius ry. The balls can penetrate into one 
another (Fig. 34a), that is, the regions of space 
that they occupy can overlap. The balls form 
strings (Fig. 34b) and regions of more complex 
shape that are regarded as metallic by definition. 
We want to find the concentration Ne of balls 
beginning with which the metal regions make the 
whole large sample electrically conductive. 

But how to choose the radius ro? A hydrogen- 
like impurity atom does not have well-defined 
boundaries. The probability of finding an electron 
ala distance r from the nucleus diminishes with 
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increasing r by the law exp (—2r/ag), where 
e =: 2.718 is the base of natural logarithm. The 
probability of finding an electron at r = aj/2 is 
less by a factor of 2.7 than that atr = 0; atr 
: it is less by a factor of 7.4 than that at r 
= (); atr = 1.5 aj it is smaller by a factor of 20. 
The effective radius of the alom must definitely 
be proportional to the length aj: 


ro= qah (6) 


This is an important statement. Tt signifies 
that the numerical coefficient g must be fairly 
universal: it varies rather weakly in going from 
one semiconductor to another, while the length 
aj varies quite considerably. 

It is very difficult to calculate g on the basis of 
physical arguments. A simpler approach is to find 
the critical concentration Ne corresponding to 
the Mott transition in terms of the model for- 
mulated above. This will give us We as a func- 
tion of the ball of radius ry. Then, making use 
of formula (6), we must express Ne in terms of q 
and aj. According to experimental data, Ne is 
given by formula (5). Comparing the theoretical 
expression with formula (5), we can find q, that 
is, determine it “as if? from experimental data. 

Let us start the realization of this program. 
We have to solve the following problem of perco- 
lation theory. Spheres of radius rọ, whose cen- 
lers are distributed in space randomly and on the 
average uniformly, are drawn in the three-di- 
mensional space. The mean number of centers of 
spheres per unit volume is N. Two spheres are 
regarded as connected if they overlap (see 
Fig. 34b). We need to find the critical concentra- 


iix 
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tion Ne at which percolation sets in through 
overlapping spheres, that is, paths appear that 
traverse the whole system and are composed of 
overlapping spheres (see Fig. 34b). 

This problem differs from the sphere problem 
formulated in the preceding chapter in that the 
connected spheres in the sphere problem were not 


Fig. 35. String of inclusive circles of radius 2r is shown 
by the dashed lines. The centers of circles are marked by 
the crosses. The percolation path through these centers 
is shown by the broken line, The percolation path through 


overlapping circles of radius rọ follows the same 
centers. 


merely overlapping but inclusive, that is, their 
centers were at a distance smaller than To, not 
than 2r,. However, this difference is not too 
significant, and the results of one problem are 
easily transferred to the other. Indeed, if an in- 
finite path exists through inclusive spheres of 
radius 2r, at a certain concentration of centers 
of spheres, there also exists a path through over- 
lapping spheres of radius To This is just the 
same path, that is, the path through the same 
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centers. This argument is illustrated in Fig. 35 
for the plane problem. (The reader will easily 
visualize the corresponding drawing for the 3D 
case.) If there is no path through inclusive spheres 
of radius 2r, there is no path through overlap- 
ping spheres of radius rọ. Hence, the critical con- 
centration corresponding to the percolation thre- 
shold through inclusive sphere of radius 2r equals 
the critical concentration corresponding to the 
percolation threshold through overlapping spheres 
of radius ro. 

According to the results given in the preceding 
chapter, the critical concentration for percolation 
through inclusive spheres of radius 2ry is deter- 
mined by the condition 


(4/3) Ne (2r)? = Be ~ 2.7 (7) 


Substituting (6) into (7), we obtain 


Neap = (8) 


The coefficient gean now be found by comparing 
(8) with expression (5) obtained from experimen- 
tal data. In order for these expressions to coin- 
cide, we have to set q = 1.6. Therefore, the effec- 
tive radius rą equals 1.6ag. The probability to 
find an electron at the point at a distance of 
r = r, from the nucleus is less by a factor of 24 
than the probability to find it at the pointr = 0. 

The most important feature of this distribution 
of the metal-insulator transition lies not so much 
in calculating the quantity q that determines the 
effective atomic’ radius but, of course, in gaining 
an insight, through percolation theory, into the 
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internal organization of semiconductors al impu- 
rity concentrations close to Ne. , 
If the suggested description is correcl, a semi- 
conductor contains a system of metallic channels 
that run across the matrix. Electric current flows 
through these channels as through wires. In 
Part ITI of this book it will be shown that the 
fraction of volume of a semiconductor occupied by 
these channels is very small if impurity concen- 
tration is close to Ne. This entails specific prop- 
erties of the electric conduction and other im- 
portant characteristics of semiconductors, 


Exercise 


1. Indium antimonide 


(InSb) is a semiconduc- 
tor with a very 


narrow gap (0.18 eV at room 
temperature), Effective masses are also small in 


such semiconductors, Assuming that the effective 
electron mass m* is 0.015 m and the dielectric 
constant is e = 18, calculate the effective Bohr 
radius ağ and the critical concentration Ne cor- 
responding to the Mott metal-insulator transition. 


Chapter 9 


Various Generalizations 
of the Sphere Problem 


Inclusive Figures of Arbitrary Shape 


We have already mentioned that the motion of 
electrons through a crystal is described by the 
effective mass m*; this m* may considerably 
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differ from the free electron mass, that is, the 
mass of an electron in an empty space. It was 
also found that the effective mass may differ from 
the free electron mass not only in magnitude. 
The point is that not all directions in a crystal 
are equivalent. For instance, the motion of an 
electron along the edges of a cube may differ 
from the motion along body diagonals and from 
the motion along face diagonals of the cube. 
Consequently, the effective mass is not necessari- 
ly identical in different directions. As a result, 
the region in which an electron moves in the 
neighborhood of a donor atom will not be spheri- 
cal. Tt may be ellipsoidal or even more complex. 

These arguments explain why the sphere prob- 
lem was generalized to the case of figures of arbi- 
trary shape. The new problem is formulated as 
follows: the sites (centers) are distributed in 
space randomly and on the average uniformly. 
The site concentration is V. The sites are sur- 
rounded by identical closed surfaces of arbitrary 
shape. 

The surfaces surrounding different sites are iden- 
tical not only in shape but also in surface orien- 
tation in space. If, for instance, the surfaces are 
fish-shaped, the tails of all the fishes must face in 


the same direction. 
The volume within one surface is V. Two sites 


are said to be connected if one of them lies with- 
in the surface drawn around the other site (inclu- 
sive surfaces). We have to find the critical con- 
centration Ne at which percolation through con- 
nected sites sets in. 

As in the sphere problem, the presence or ab- 
sence of percolation is determined only by the 
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value of the parameter B which equals the pro- 
duct VN. This is the mean number of sites within 
the volume bounded by one surface. This number 
will not change if we enlarge or contract the 
scales along all the three directions, that is, if we 
multiply or divide the coordinates of all sites in 
the system and of all points on the surfaces by 
the same factor. It is also clear that this rescal- 
ing will not suppress percolation if it existed; 
neither will it bring it to existence. Therefore, 
percolation “does not respond” to those changes 
in N andY that leave B unaltered. Consequently, 
it will be more convenient, as in the sphere prob- 
lem, to argue not in terms of the critical concen- 


tration M, but in terms of the critical parameter 
o 


Be = N,V (1) 


around 
volume V 
the critical 
arge, while the 
ed. It depends 


If we keep the shape of the surface drawn 
each site constant but increase the 
bounded by this surface, say, twice, 
concentration NV, will be half as 1 
parameter B, will not be effect 
only on the surface shape. 
Formula (1) generalizes the formula 


Bo = (4/3) nN, R (2) 


that we used earlier for the sphere problem. 

In the general case, the critical value of B, at 
which percolation sets in must not be equal to 
the value of B. = 2.7 obtained for the sphere 
problem. At present the value of Be for different 
figures is a subject of intensive study. 
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The Ellipsoid Problem 


The Soviet scientists Ya. G. Sinai and B. I. 
Shklovskii proved that certain distinct surfaces 
have the same value of Be. Thus, the ellipsoid and 
sphere (and the ellipse and circle in the plane 
case) are such surfaces. 

First we remind the reader of the definitions of 
ellipse and ellipsoid. The ellipse is a closed curve 


Bi 


Fig. 36. Ellipse Æ obtained from circle O by extending 
it along the z-axis. 


drawn on a plane obtained from a circle by extend- 
ing (or contracting) it along one of its axes 
(Fig. 36). In order to realize this extension or 
contraction, it is necessary to transform each 
point M of the circle with z-, y-coordinates into 
the point M’ with coordinates x! = kr, y =Y, 
where jt is the extension factor (k > 1 corresponds 
to extension, and k < 1 to contraction). 

The surface obtained from a sphere by exten- 
sion (or contraction) along one of its axes is 
called the ellipsoid of revolution. In order to realize 
this extension or contraction along the z-axis 
it is necessary to transform each point on the 
sphere with z-, y-, z-coordinates into the point 
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with coordinates 2’ = z, y’ = J, 2° = 
figure thus obtained is shown in Fig. 

called the ellipsoid of revolution because a rotla- 
tion by any angle around the z-axis (called the 


Fig. 37. Ellipsoid of revolution, 


axis of revolution) transforms this solid into 
itself. If it is cut by a plane drawn through the 
z-axis the section is an ellipse. The section by a 
plane drawn perpendicularly to the z-axis is a 
circle. 

The general-type ellipsoid is obtained from an 
ellipsoid of revolution by extending (or contract- 
ing) it along one of its axes perpendicular to 
the axis of revolution. The new extension factor 
is not necessarily equal to X. The section of the 
gencral-type ellipsoid by any plane parallel to 
the zOy-, xOz-, yOz-planes is an ellipse, 


Ch. 9. Generalizations of the Sphere Problem 174 


Now we shall prove that ellipsoids and spheres 
have identical By. Let the sites be randomly 
distributed in space with the mean concentration 
N. Draw a sphere of radius R around each site. 
If B = (4/3) xR? > BY" there is percolation, 
otherwise, there is no percolation through con- 
tacting spheres. Here BS" is the critical value 


of B for the sphere problem (B® = 2.7 + 0.1). 

, Let us apply the extension y’ = y, 2 = kt, 
z' = kz along the z- and a-axes in order to 
transform the coordinates of both the sites and 
points on the spheres. (If we dealt with a plane 
problem, the picture could be modelled by mark- 
ing the sites and circles on a rubber band which 
is then extended in one direction. Likewise, we 
could imagine a three-dimensional “rubber space” 
in which we have marked the sites and spherical 
surfaces. Then this “rubber space” is extended 
in two directions.) 

If the 2-coordinates of sites are random numbers 
uniformly distributed within the interval from 
0 to L, where L is the size of the system, the new 
x'-coordinates obtained by multiplying x by ky 
are also random numbers uniformly distributed 
within the range from 0 to kZ. The z'-coordinates 
are uniformly distributed within the interval 
from 0 to kL, while the y-coordinates remain 
unaltered. The concentration N’ of sites is not 
equal to V. see is 

All spheres now transformed into ellipsoids 
with volume V’. (It can be shown that V’ = 
= kka. (4/3) nR. This is irrelevant for the 
derivalion to follow.) As a result, the sphere 
problem converted to the ellipsoid problem, Here 
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we can introduce a quantity B’ = N’V’ and 
find the critical value of B’ at which percolation 
through ellipsoids sets in. Let us denote it by B% 

The rest of the proof separates into the following 
steps: 

f B = P'. All sites that lay within a certain 
sphere prior to transformation are found inside 
the ellipsoid obtained from this sphere after 
the transformation. Indeed, when the sphere is 
extended (or contracted), the points internal to 
the sphere remain internal, and those external to 
it remain external at each step of extension. 
Hence, the mean number of sites, B, which fell 
inside one sphere prior to transformation equals 
the mean number of sites, B’, which fell inside 
one ellipsoid after the transformation. 


2. If B>B, then B>B” (3) 
If B< Be, then BaB (4) 
Indeed, if two spheres were connected prior to 

transformation, that is, the center of one of them 


was inside another, the two ellipsoids formed 
out of these spheres are also connected because 


surfaces. If two spheres were not connected, the 
ellipsoids obtained from them are not connected 
either. This yields that if p > Be", that is, 
if infinite percolation paths exist through con- 
nected spheres, infinite percolation paths also 
exist through connected ellipsoids, and this sig- 
nifies, in turn, that B’ is greater than the threshold 
value B8", This gives us condition (3) because 


B' =B. It B< Be", that is, no percolation 
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exists through connected spheres, there is no 
percolation through connected ellipsoids, that is, 
B' < Be". This gives us condition (4). 

3. Since conditions (3) and (4) must hold for 
any value of B, they imply that Bg" = get. 
which was to be proved. 


Other Surfaces 


The ellipsoid problem is very important in semi- 
conductor physics, but in some cases the ellip- 
soid is too simple a surface. Physicists were there- 
fore much interested in how the critical value of 
B, depends on the shape of the surface when exact 
relations cannot be obtained. To gain this knowl- 
edge, Monte Carlo computer simulations were 
carried out, and the following surprising and 
extremely important fact was found out: Be is 
only very weakly dependent on the surface 
shape. 

The shapes that were analyzed were a cube and 
a tetrahedron. It was found that the critical 
values for these solids do not differ, within the 
error of calculations (+0.41), either from each 
other or from the value of B, for the sphere prob- 
lem. Attempts to find a solid with least resem- 
blance to a sphere led to a “three-dimensional 
cross”, that is, a figure formed by three elongated 
parallelepipeds intersecting at the origin of 
coordinates (Fig. 38). It was found that B, for 
this surface is only 20% less than that for the 
sphere. 

The same conclusion was obtained for plane 
figures. An analysis demonstrated that B, for 
squares differs from Be for circles by mere 2%. 
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Consequently, Be is, lo a satisfactory accuracy, 
universal within a class of figures of identical 
dimensionality (i.e. either plane or three-dimen- 
sional). The reason for this observation may be 
a certain relation similar to that proved in the 
preceding section, so that exact equalities in 
fact hold where computer simulation discovers 
only a slight deviation or no deviation al all, 


Fig. 38. “Three-dimensional cross”, 


within computational errors. (Small deviations 
found by computer calculations must be looked 
al critically because it is not always easy to carry 
oul a correct evaluation of results.) Unfortunately, 
this is all we know at present. 


Another Experiment at the House Kitchen. 
The Hard-Sphere Problem 


Tn 1974 three students of the Harvard University 
in the USA, J.B. Fitzpatrick, R.B. Malt, and 
F. Spaepen, conducted the following simple 
experiment. Five thousand small balls, part 
of them made of aluminium and another part of 
plastic, were placed in a jar. Before starting the 
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measurements, the balls were thoroughly stirred 
and then the jar was intensely shaken to make 
the packing of balls maximally tight. A metal 
foil electrode was placed at the bottom of the jar, 
and another electrode was placed on top of the 
balls, afler which the system was compressed by 
a load of 15 kg. 

Aluminium is a highly conductive metal, and 
the plastic was an insulator. The aim of the 
experiment was to determine the critical fraction 
x, of aluminium balls at which electric current 
flows between the electrodes, that is, paths appear 
through the aluminium balls in contact. 16 was 
found that a, & 0.25. 

In addition, the experiment made it possible 
to study the electric conductivity of the system 
as a function of z at x > £e. 

Here we encounter another problem of per- 
colation theory, namely, the hard-sphere problem. 
The closest relation to this new problem is per- 
haps the site problem. Let us recall our approach 
to an approximate evaluation of percolation 
threshold in the site problem used in Chapter 6. 
A sphere is constructed around each lattice site, 
its radius being equal to half the distance to the 
nearest neighbor. The spheres constructed around 
white sites are said to be white, and those around 
black sites are said to be black. Percolation 
through white sites is equivalent to the existence 
of percolation paths through tangent white 
spheres (see Fig. 31). 

The difference between the site problem and 
the new problem lies in that the centers of spheres 
in the site problem are the sites of a regular lat- 
tice, while in the new problem the centers can be 
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located at any point. It will be shown below 
that this difference is not very important. 

The difference between the new problem and 
the sphere problem discussed in the preceding 
chapters and the preceding sections of this chap- 
ter lies in that the spheres that we operate with 
in the new problem are assumed to be hard and 
thus cannot overlap. This is an essential distinc- 
tion. 

It was shown in Chapter 6 that in the sile 
problem percolation through white sites sets in 
when the fraction of space filled by the spheres 
constructed around these sites is approximately 
0.16. It was found that this number is almost 
independent of lattice type. It is then natural 
to conjecture that if it is nearly independent of 
lattice type, it should not strongly depend on 
whether the lattice is or is not there. If this hypo- 
thesis is correct, the fraction of volume filled 
with metallic balls at which percolation appears 
through these balls must be approximately 0.16. 

Let us denote by f, as in Chapter 6, the filling 
factor, that is, the fraction of volume occupied 
Jointly by both aluminium and plastic balls. 
By definition, the quantity x is the ratio of the 
number of aluminium balls to the total number 
of balls. Hence, the fraction of volume occupied 
by aluminium balls equals fz. If this fraction 
of volume equals 0.16 at percolation threshold, 
the critical value Te can be found from the con- 
dition fx, = 0.16. 

The filling factor of a system of tightly packed 
but randomly distributed balls is well known. 
Such a system is familiar to mankind from the 
times of antiquity. If it was necessary to measure 
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a certain amount of grain or some other free- 
flowing material, it was poured into a special 
vessel (a measure), shaken, and compacted. In 
modern science this system is a model of atomic 
arrangement in amorphous metals. Amorphous 
metals, also referred to as metallic glasses, are 
materials with metallic electric conductance but 
having no crystal structure. It was found that 
the atomic arrangement of amorphous metals 
strongly resembles the arrangement of tightly 
packed incompressible balls. For this reason, 
the properties of tightly packed randomly dis- 
tributed balls were studied very thoroughly 
(mostly in computer models). It was found that 
the fraction of volume filled with the balls is 
j = 0.637. 

Let us return to the problem of percolation 
through metallic balls. Having determined Te 
by means of the formula xe = 0.16/f, we obtain 
ze = 0.25, in complete agreement with the 
result of the three students. 

Experiments on determining the percolation 
threshold were often repeated, cach time with 
better instruments. Experiments were conducted 
in which balls in the mixture were of different 
radii. The radii of both metallic and dielectric 
balls were varied in the mixture in a wide range. 
It was found that in this case the critical volume 
fraction fđe is approximately 0.17, being almost 
the same, to within experimental accuracy, as 
in the case of identical balls. 

Percolation threshold in the hard-sphere prob- 
lem can thus be evaluated rather easily after 
we assume that the critical volume fraction filled 
with metallic balls is approximately 0.16. Note 
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that this fraction is much greater in the problem 
of overlapping spheres. 

The hard-sphere problem proved very impor- 
tant for applications as well. It accounts for the 
theory of heterogencous materials manufactured 
from an insulator with tiny metallic inclusions. 
Such materials are now an imporlant object of 
study. They are prepared and utilized both as 
thin films and as bulk samples. The electrical 
properties of such materials in the vicinity of 
percolation threshold are fascinating. For in- 
stance, the capacitance of a capacitor filled with 
such a material grows to infinity when the frac- 
tion of volume filled with the metal tends to per- 
colation threshold. This phenomenon stems from 
an enormously high mutual capacitance of large 
metallic clusters. The description of electrical 
properties of heterogeneous materials is growing 


nowadays into an independent domain of per- 
colation theory. 
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Chapter 10 


Percolation Level 


“The Flood” 


Six days and [six] nights 
Blows the flood wind, as the south-storm 
sweeps the land. 
When the seventh day arrived, 
The flood (-carrying) south-storm 
subsided in the battle, 
Which it had fought like an army. 
The sea grew quiet, the tempest 
was still, the flood ceased. 
I looked at the weather: stillness 
had set in, 
And all of mankind had returned to clay. 
The landscape was as level 
as a flat roof. 
I opened a hatch, and light fell upon 
my face. 


In: Ancient Near Eastern Terts 
Relating to the Old Testament, 
Ed. by J. B. Pritchard, 
Princeton Univ. Press, 1955, 
p. 94 
(Accadian Myths and Epics, 
the Epic of 
Gilgamesh, 

Gilgamesh, Tablet XI, trans- 
lator 
E. A. Spaiser) 


This description of the Flood is found in the 
Babylonian Epic of Gilgamesh, that most an- 
cient work of literature which dates back roughly 
to the second millennium BC. 


1Q* 
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When the waters began to subside, the sum- 
mits of the highest mountains appeared above 
the surface. The water subsided lower and lower 
until it reached its normal level. Just imagine 
this vast picture: a huge system of mountain 


Fig. 39. Voyage during the Flood. 


ridges, such as the Himalayas, gradually emerged 
from the water. First the highest peaks broke 
the surface and formed islands, then the zone of 
alpine meadows was liberated, and finally, water 
dropped to foothills, 
Suppose that we want to 
To what level must the wa 


the last waterway across the whole system of 
ridges to disappear (Fig. 39)? Clearly, such a 
way exists as long as a certain fraction of moun- 


answer the question: 
ter drop in order for 
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Fig. 40. Map of a mountain system flooded with water. 
The submerged regions are shaded. 
lainous passes remain under water, but then it 
dries up. 

Figure 40 illustrates a geographic map of a 
mountain system on which the thin lines trace 
the contour lines (curves of identical height). 
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The thicker lines mark the contours corresponding 
to water level. They separate water from dry 
land. The submerged regions are shaded in the 
drawing. Figure 40a shows the very beginning of 
the Flood: the lakes do not communicate with 
one another. Then the water level climbed, and 
Fig. 40b now shows the waterways across the 
system. And in Fig. 40c only individual summits 
are above the water, 

The level of water at which this waterway 
appears (or disappears) is called percolation level. 
The problem of finding this waterway is a plane 
problem of percolation theory. It can be refor- 
mulated. Let us assume that a plane is randomly 
painted by white and black paints. Let the frac- 
tion of area painted white be x. At small values 
of x, while spots form isolated islands, while at 
w nearly equal to unity, the isolated spols are 
black. We want to find the critical value of x 
at which a noninterrupted path appears (or dis- 
appears) across the whole system, going only 
through white regions, 

Likewise, we can also formul 
ing 3D problem, filling a large volume with a 
white and a black substances. Then we have to 
vary the fraction of volume occupied by one of 
these substances until percolation sets in. 


ate the correspond- 


How to Construct a Random Function ** 
This new problem is nota | 
us analyze its mathematical formulation in the 
plane case. We have to define on the whole plane 
a random function V (X, Y), where X and Y 
are coordinates. In the Flood problem this func- 


altice problem. Let 
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tion is the altitude above sea level of the point 
on the earth’s surface with the given coordinates. 
We shall confine ourselves only to the case of 
Gaussian random functions (named after the great 
mathematician Karl Friedrich Gauss). 

The simplest method of constructing a Gaus- 
sian random function is to ascribe to each point 
in space a random number not related in any 
way to a neighboring random number. Such 
random functions were given the name “white 
noise”. The values it takes at neighboring points 
differ abruptly, that is, this function is discon- 
tinuous. To obtain a continuous function, the 
white noise has to be “smoothed”. This proce- 
dure consists in ascribing to each point of space 
a quantity equal to the average of the values that 
the “white noise” function assumes in a certain 
region around this point. The continuous random 
function is then formed by these values. The 
function varies only slightly from this point to 
the neighboring point because of the small size 
of the spatial region over which the white noise 
is averaged. Let us denote by rọ the size of the 
region over which averaging is carried out. This 
size is referred to as the correlation radius of a 
random function. Its main property is that the 
function changes only slightly when the argu- 
ment of the function changes by an increment 
small compared with ro 

The quantities V can be described by a dis- 
tribution function that we denote by { (V). By 
definition, the probability for the function 
V (X, Y) to have at a point taken randomly in 
space a value within a small interval from V, 
to V, + AV equals f (V,) AV. 
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The function V (X, Y) can always be con- 
structed in such a way that, averaged over all 
points of a plane, it gives zero. To achieve this, 
we need to use in the construction of the “white 
noise” function the random numbers distributed 
symmetrically with respect to zero. It can be 
shown that the function V (X, Y) obtained from 
such “white noise” has the Gaussian distribution 
(see Chapter 2) of the type 

4 . y2 1 
10 = Spee ( aa) (4) 


Note that the probability density (1) is sym- 
metric relative to the positive and negative val- 
ues of V. In “mountain lingo” this means that 
peaks and valleys are encountered with equal 
probability, 

Now we are ready to enunciate the percolation 
problem. The contours shown in Fig. 40 are de- 
termined by the condition V (X, Y) = const. In 
order to prescribe the water level, we need to 
introduce a quantity ¢ that can vary from —oo 
to -+-co. The regions on the plane where V (X, Y)< 
<? will be referred to as white (covered with 
water), and those where V (X, Y) >t as black 
(protruding above water). The thick contours in 
Fig. 40 correspond to V (X, FaR 

The percolation level is defined as the critical 
value t, at which white regions form paths 
going uninterruptedly across the whole system. 
It is also possible to speak of the critical fraction 
Te ol space filled with while regions al the moment 
when percolation sets in. This fraction of space 
equals the Probability for a continuous random 
variable V to assume a value in the range —coo < 
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< V < te- By virtue of the definition of distribu- 
tion functions, 
t 


w= | s(v)av (2) 


Formula (2) relates the critical fraction ze of 
space to the percolation level tg. 


Analogy to the Site Problem ** 


Imagine that a plane lattice (its type is immate- 
rial) with small period is superposed over the 
geographic map shown in Fig. 40. Let the lattice 
period be substantially smaller than the charac- 
teristic extensions of land and water regions. 
Let us “paint” the lattice sites in water-covered 
areas white, and those in land areas black. As in 
the familiar site problem, we again consider 
white sites connected if they are nearest neigh- 
bors. The critical fraction of area under water 
equals the critical fraction of while sites at which 
percolation through white sets in. 

However, the new problem is not identical to 
the site problem. Remember Chapter 4 where 
the construction of a system of white and black 
sites has been described in detail. (Sites were 
labelled in that chapter as blocked and non- 
blocked.) It was an important aspect of that 
construction that each site became black or 
while as prescribed by a random-number genc- 
ation to the color of the sites in 
the nearest neighborhood. As a result, the white 
and black sites were thoroughly mixed. The 
situation is quite different in the new problem. 


valor, in no rel 
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The lattice period being very small, black and 
white sites form large groups. Neighbors of a 
white site are almost certainly white, and those 
of a black site are most probably black. (The 
size of blocks is determined by the correlation 
radius ry.) , 

It must be clear that to introduce a lattice into 
the percolation level problem is a purely formal 
operation. This is in fact a nonlattice problem, 
and its solutions must be independent of both 
the lattice period (provided it is sufficiently 
small) and lattice type. However, this method 
allows us to make use of a well-developed appa- 
ratus of lattice problems. 

In Chapter 4 we gave an algorithm of the Monte 
Carlo solution of a site problem with a computer. 
This algorithm is entirely transferable to the 
percolation level problem. An array K (X, Y) 
consisting of zeros and unities is generated 
by addressing a random function V (X, Y) by 
the method described in Chapter 4. White sites 
correspond to unities, and black siles to zeros. 
Then percolation paths are searched for, and the 
critical fraction of space at which percolation 
sets in is determined by the same lechniques, 


Percolation Levels 
in Plane and Three-Dimensional P 


A plane problem has an exact solution if the 
properties of a random function V (X, Y) are 
on the average symmetric with respect to V — 0. 
The Gaussian functions described 
among the functions with 

In order to arrive 


roblems ** 


above are 
such properties. 
at the exact solution, we 
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must use the symmetric formulation of the 
percolation problem given in Chapter 5. As be- 
fore, we refer to the regions in which V (X, Y) < t 
as while, and to those in which V (X, Y) > 
as black. In addition to percolation level through 
white regions, tẹ, we can introduce percolation 
level through black regions, te. The further proof 
goes by the following steps: 

1. By virtue of the symmetry of the function 


V (X, Y), we have te = —t. Indeed, replace V 
by —V at each point of the plane. This will 
give a function V’ = —V that on the average 


has the same properties, so that percolation 
levels calculated by V’ must be equal to those 
calculated by the original function V. The 
inequality V < to yields the inequality —V > 
> —t,, that is, V’ > —te. At t = t, the regions 
that are white with respect to the function V 
(V < te) form an infinite cluster. These same 
regions of space are black with respect to the 
function V’ and t = —tg. Indeed, here V’ > 
> —t,. Therefore, the value ¢ = —t, is per- 
colation level through black regions for the 
function V’. But it has already been said that 
the functions V and V’ should have identical 
percolation levels. Consequently, percolation level 
through black regions for the function V is 
k = —ty. 
° 2. if k < 0, then, as ¢ increases, first percola- 
tion through white regions appears (at t = 
= te <Q) and then (att = —te > 0) percolation 
through black regions disappears. In the range 
tp << t< —t, there is percolation both through 
white and black regions. If t œQ, first per- 
colation through black regions disappears, and 
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only then percolation through while regions 
sets in. And in the range —t,<l< te per- 
colation is impossible. 

3. In plane problems, percolation through 
black regions excludes percolation through white 
regions, and vice versa. Indeed, if we can sail 
across a mountain ridge from west to east, ib 
means that we cannot cross it by land from 
north to south. This excludes the case be <0. 
On the other hand, the absence of percolation 
through white regions necessarily means the 
presence of southward percolation through black 
regions. (You can make sure of that by analyzing 
pictures like Fig. 40.) Therefore, the case te > 0 
is also excluded. This leaves us wilh the single 
possibility: te = 0. This is the sought result: 
percolation level stands at zero. 

Formula (2) makes it Possible to calculate the 
critical fraction ze of area. As follows from the 
normalization condition for distribution func- 
tion (see formula (1) of Chapter 2), at ta = 00 
the right-hand side of formula (2) equals unity. 
The symmetric nature of function Í (y) implies 
that at te = 0 the fraction of area z, = 0.5. 

In three-dimensional cases, eastward percola- 
tion through white regions does not preclude 
southward percolation through black regions 
because percolation channels can be easily un- 
coupled. (Remember highway bypasses arranged 
al different levels.) For this reason, in the 3D 
case, le < O and, by virtue of formula (2); Bs =< 
< 0.5. Monte Carlo calculations by the above- 
described procedure for Gaussian random func- 
lions demonstrated that the 3D crilical fraction 


Ze = 0.16 + 0.01. 
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The quantity £e can be approximately evalu- 
ated by the method we used to evaluate per- 
colation thresholds of the site problem in Chap- 
ter 6. The method consists in surrounding each 
site by a sphere (or circle in the plane case) of 
radius equal to the distance to the nearest neigh- 
bor. A sphere around a white site is regarded as 
white, and a sphere around a black site is regar- 
ded as black. It was found that percolation 
through white spheres in contact appears when 
the fraction of space filled with these spheres is 
approximately the same for all lattices. It is 
natural to assume that this fraction must not be 
very different from the value Te that we en- 
counter in the percolation level problem. Accord- 
ing to formulas (3) and (4) of Chapter 6, in the 
2D case the fraction of area occupied by while 
circles is approximately 0.5, and in the 3D case 
the fraction of volume occupied by white spheres 
is approximately 0.16. The two estimates thus 
coincide strikingly well with the solution of the 
percolation level problem. We can expect that 
z, = 0.16 is a fairly good estimate for non- 


Gaussian random functions. 


Impurity Compensation in Semiconductors 


Percolation level plays a veryzimportant role in 
the theory of extrinsic semiconductors. Assume 
that a semiconductor was doped by both donor 
and acceptor impurities in identical amounts. 
Donor impurity atoms have an extra electron 
on the outer electron shell, while, contrary to 
this, acceptor impurity atoms lack one electron. 
As a result, donors eagerly donate their extra 
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electrons to acceptors and thereby become posi- 
lively charged. Acceptors accept electrons and 
become negative. (This phenomenon is called 
impurity compensation.) Impurities in semicon- 
ductors being distributed chaotically, a random 
system of positively and negatively located char- 
ges is formed. Bach charge produces an electric 
potential --e/er, where e is the absolute magnitude 
of the electron charge, e is the dielectric constant, 
and r is the distance from the charge. The sign of 
the potential depends on the sign of the charge. 
The potential of any point in space is the 
sum of the potentials produced by all impurity 
atoms. The impurities being distributed ran- 
domly through the semiconductor, their poten- 
tial is also a random function. 
If the concentration of donors is slightly higher 
than that of acceplors, some electrons remain on 
the donors. If the donor-electron bonding energy 
is relatively low, the thermal motion of atoms 
easily ionizes the donors. In principle, the liber- 
ated electrons can participate in current transfer, 
but this transfer is inhibited by an electric poten- 
tial produced by charged impurity atoms within 
the semiconductor. The product of electric poten- 
tial by the electron charge is the potential energy 
of interaction between an electron and the electric 
field of impurities. Potential energy essentially 
changes the character of motion of electrons. 


Motion of a Particle 
with Nonzero Potential Energy 


The total energy Æ of a P 


article consists of its 
kinetic energy mv?/2, wher 


e m is the mass of the 
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particle and v is its velocity, and potential energy 
V (r) which is a function of coordinates of the 
point at which the particle is located: 


b=" + (1) (3) 


where r is a vector from the origin of coordinates 
to the location of the particle. 

The trajectory of a particle is described by a 
function r (t), where ¢ is time. The fundamental 
principle of mechanics is that the total energy E 
of a particle does not change in the course of 
motion. This means that a change in coordinates 
of the particle entails a change in its velocity 
so as to compensate for the change in the poten- 
tial energy V (r). The law of energy conservation 
and the fact that the kinetic energy mv*/2 must be 
a positive quantity impose important restrictions 
on the motion of the particle. 

Let us assume for the sake of simplicity that V 
is a function of only one coordinate, namely, X, 
and that velocity is also directed only along the 
X-axis. Let V (X) have the form shown in Fig. 4A, 
The type of motion of the particle is determined 
by its total energy. If the total energy is E,, 
the kinetic energy is positive only in the range 
X, < X < Xj. Points X, and Xj are the bounds 
of motion. Velocity vanishes at these points, hence 
the term turning points. A particle with energy 
E, is trapped into the potential well between the 
turning points and cannot leave this well. Its 
range of motion is bounded on both sides. A par- 
ticle with energy /, cannot penetrate deeper 
on the left than point X>, but its range of motion 
is not bounded on the right. 
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The motion of a particle can thus be bounded or 
unbounded depending on the relative values of 
the total and potential energies. 

Note that quantum mechanics allows a particle 
to penetrate into regions with a negative kinetic 


v 


Fig. 41. Potential energy V of an electron as a function 
of the X-coordinate (the solid curve). 


energy. However, if V (r) is a sufficiently smooth 
function (and in semiconductors with high con- 
centration of compensated impurities it is always 


emooti, this penetration plays only a minor 
role. 


Motion of an Electron 
in the Field of Impurities 


The potential energy of an electron, stemming 
from the potentials of randomly distributed im- 
purity atoms, is a random function of coordinates. 
If the total energy Æ of the electron is low, it 
can move only in a bounded region of space 
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everywhere surrounded by turning points. If the 
energy is sufficiently high, the electron is allowed 
to move throughout the whole space. This is the 
necessary condition for participating in the 
transfer of electric charge (in electric current). 

The reader must understand the difference 
between the one-dimensional problem discussed 
in the preceding section and two-dimensional or 
three-dimensional problems. If a particle can 
move only along a straight line, then the motion 
can be unbounded in both directions only if the 
energy of the particle is greater than all the values 
assumed by the potential energy along this straight 
line. Hence, a stringent condition Æ > Vex 
must. be satisfied (Vmax is the maximum value of 
potential energy). 

No such condition is required in the two- and 
three-dimensional cases. A particle can bypass 
the regions in which its motion is forbidden. It is 
only necessary that the regions where the motion 
is allowed form a system of “lakes and canals” 
through which the particle could travel to infinity 
in an infinite system. 

Now it is obvious that the problem that emerges 
is one of finding percolation level. Let V (r) be 
a random function describing the electron poten- 
tial energy. We fix the total energy Æ and refer 
to spatial regions where E > V (r) (the positive 
kinetic energy) as white, and to those where 
E< V (r) as black. We want to find percolation 
level, that is, the critical value Æe at which 
percolation through white regions sets in. 

Only those electrons whose energy exceeds X'e 
are free and take part in current transfer (some- 
limes this energy is called the mobility edge). 


13-0240 
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At low temperatures, Æe is substantially greater 
than the energy of thermal atomic motion kT. 
Consequently, the probability for an electron 
to reach the level of Æe is low. Correspondingly, 
the concentration of electrons capable of charge 
transfer is small. As temperature increases, this 
concentration grows steeply, thereby increasing 
the electric conductivity of the system. The per- 
colation threshold Æ, thus determines the tem- 
perature dependence of the electric conductivity 
of impurity-compensated semiconductors. 

At very low temperatures, electrons accumulate 
in the deepest potential wells and cannot take 
part in charge transfer. Consequently, at low 
temperatures, impurity-compensated semicon- 
ductors turn into insulators. 

It was mentioned in Chapter 8 that at a suf- 
ficiently high concentration of impurities of 
one sort (e.g. donors), the semiconductor’s con- 
duction becomes metallic, only marginally de- 
pending on temperature down to absolute zero 
(the Mott transition). Impurity compensation 
(e.g. by adding acceptors) enhances the random 
potential energy and switches off the metallic 
conduction. A detailed theory of this phenomenon 
has been developed on the basis of percolation 
level concepts. 


Part IIT 


Critical Behavior 

of Various Quantities 

Near Percolation Threshold. 
Infinite Cluster Geometry. 


In this part we discuss those aspects of percolation 
theory that appear as the most interesting from 
physics’ standpoint: the behavior of various 
quantities in the immediate neighborhood of 
percolation threshold. It was said in the previous 
parts that such physical quantities as the spon- 
taneous magnetization of a doped ferromagnetic 
or the electric conductivity of a network with 
blocked sites vanish ab the threshold point. In 
the present part we discuss the laws that describe 
their behavior in the vicinity of percolation 
threshold. The derivation of these laws requires 
that the geometrical properties of infinite clusters 


be understood. 


Chapter 11** 
The Bethe Lattice 


Tt was shawn in Chapter 5 that exact solutions 
can be found in some plane percolation threshold 
problems. However, it was never said that it 


43% 
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was also possible to find the function P (x), 
that is, the probability for a site to belong to 
an infinite cluster. No exact expressions for this 
function (or for the electric conductivity of a 
network) are known at present cither for plane 
or for three-dimensional problems. The only 
exception is the Bethe lattice which, as we shall 
show below, must be classified as a lattice in the 
space of infinitely high dimensionality. In whal 
follows we will give the formulation and solution 
of the site problem on the Bethe lattice. 


Rumors 


“What nonsensical stories they do spread aboul 
the town. What are things coming to when you 
can hardly turn round before there is some 
scandal going about you, and not a word of sense 
in it either...”*—in these words Nikolai Gogol 
described how a preposterous gossip spread by 
two ladies ruined Chichikov’s promising fraud. 
“This enterprise they contrived to carry out in 
just a trifle over half an hour. The town was posi- 
tively stirred up; everything was in a ferment; 
and if there were but anybody that could make 
out anything!”** 

And indeed, rumors do spread unbelievably 
fast. But this speed ceases to be surprising if 
we analyze the mathemati 


cal model suggested 
below. 


* Dead Souls 
London, 1922; 
Garnett. 

** Dead Souls, a poem by N. G 
New York, 1944; transl 
erney. 


» 4 poem by N. Gogol. Chatto & Windus, 
translated from the Russian by Constance 


ogol. The Heritage Press, 
ated from the Russian by B., G. Gu- 
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tb) 
attice for g = 3. The open circles 
O, and the black circles re- 


present those of category Ts 


Fig. 42. The Bethe 1 
represent people of category 


The model is shown in Fig. 42a. Assume that 


a “lady agreeable in every respect”, denoted in the 
figure by circle A, passed the news lo three of her 
acquaintances, B1, By, By. Each of these threo 
passed the message on to three of her (or his) 
rumor-mongering acquaintances, SO that this 
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“second-hand” information was received by nine 
people marked by circles C. Each of these nine 
people again transmitted information to their 
three contacts, making the news known to another 
27 people. It is readily found that the “tenth- 
hand” information reaches 3! = 59 049 people! 
Assuming that it takes each gossip twenty min- 
utes to pass the information to his (her) three 
listeners, we find that the whole chain takes 
200 min = 3 h 20 min. 

Of course, this model strongly simplifies the real 
process. It assumes that all people have an iden- 
tical number of acquaintances. Furthermore, it 
assumes that each person receives information 
from only one informant. This means that only 
one line enters each circle (see Fig. 42a), By vir- 
tuc of this property, the model resembles a tree 
that branches out infinitely in all directions. 
Kach circle can be regarded as a trunk of its 
own tree, and the trees that grew from, for in- 
stance, circles B,, B}, Bs have no common circles. 
The same can be said about the trees stemming 
from circles C, and so on. 

In scientific literature, this model is called 
just that: a tree. It is also called the Bethe lattice, 
after the famous physicist Hans Bethe. The 
circles shown in Fig. 42a are the sites of this 
lattice. The number of lines originating from 
each sile of the Bethe lattice can be arbitrary 
(but identical for all sites). Let us denote this 


number by q. Figure 42a shows the lattice for 
2 


= 3. 
Now recall that most people have their own 


point of view and do not participate in spreading 
rumors. Let us divide all people into two cate- 
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gories: category O shown by the open circles 
consists of people who transmit the received in- 
formation to the next trio (Fig. 42b). These circles 
give rise to q arrows. Category T shown by the 
black circles consists of people who do not par- 
ticipate in spreading the rumors. No arrows 
originate from the black circles (sce Fig. 42b). 

The introduction of the black circles strongly 
affects the propagation of a rumor. Consider the 
configuration shown in Fig. 42b. Among the 
three circles B,, Ba, Bs, only one is white and 
transmits the rumor further on. Circles Ci, Cs, 
C, would be glad to gossip but received nothing 
from B, and By. Among C4, Cs, Cg, only Ce belongs 
to O. Therefore, instead of nine recipients of the 
second-hand information, only three were reached 
by the rumor, and only C, will transmit it along 
the line. 

Imagine that the system under discussion is 
not bounded and thus has an infinite number of 
circles. Then the following question can be 
posed: Will a rumor originating from point A 
die out after a finite number of transfers, or 
will it stray to an infinite distance from A and 
become known to an infinite number of persons 
in an infinite system? We see from Fig. 42b that 
this depends on the relative number of open and 
black circles and on configurations that arise 
around a sile. 

In fact, we speak of a sile problem of per- 
colation theory, only formulated on the Bethe 
lattice. Let the fraction of people belonging to O 
be x. This means that a person selected at random 
will be found to belong to category O with 
probability z and to category T with probability 
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1 — z. The question we have to answer is as 
follows: What is the probability P (x) for a rumor 
transmitted to a randomly selected person to 
become known to an infinite number of people? 
Obviously, this probability is zero if x is small, 
but becomes greater than zero beginning with a 
certain critical value £ = zp. 


Solution of the Site Problem 
on the Bethe Lattice 


It will be convenient to introduce, instead of 
P (x), the probability for a rumor transmitted 
to a randomly chosen person not to become known 
to an infinite number of persons. Let us denote 
this probability by Q (x). Obviously, 


Q (2) =1 — P (2) (1) 


because these two events form a complete system. 

An algebraic equation can be written for 
Q (£). The following line of reasoning must be 
used. The spreading of a rumor can be stopped 
for two incompatible reasons. The first is that the 
person chosen at random belongs to category T. 
The second reason is that although this person 
belongs to category O and passes the rumor to q 
people, all channels leading from these people 
will be cut at different stages. The probability Q 
is thus the sum of the probabilities of two in- 
compatible events. The probability for a random- 
ly selected person to belong to category T is 
1 — zx. The probability that the person belongs 
to category O but that the rumor propagation 
will be interrupted at later stages is denoted, for 
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the time being, by W’. Then 
C=1—e+4 W 


u Now let us take a look at W’. The outcome 
lat it represents is the result of two simulta- 
neously realized events: (i) the person chosen 
randomly belongs to O (the probability of this 
event is x), (ii) all q channels leading from the 
acquaintances of the person chosen at random 
interrupt at some step. Obviously, these two 
events are independent. Hence, the probability 
W’ equals the product of the probabilities: WwW = 


= xW (x), so that 
Q =1 — z + zW (2) (2) 
y for all q channels 


where W (x) is the probabilit. 
to be interrupted at some step (of course, this 
can happen at different steps in different. chan- 
nels). 
Let us consider one of the q channels that ori- 
ginate from one of the acquaintances of the 
person chosen at random. The event consisting 
in this channel being interrupted at some step 
is equivalent to the statement that the rumor 
communicated to this acquaintance does not 
reach an infinitely large number of people. By 
definition, the} probability of this event is Q (2). 
Tt is very important for what follows that the 
trees rooted at q acquaintances of a randomly 
selected person have no common circles. Hence, 
a certain configuration of open and black circles 
found in one tree does not affect at all the pro- 
bability of any possible configuration of circles 
in other trees. (Obviously, this last statement 
would be incorrect if the trees had common cir- 


202 Part III. Quantities Near Percolation Threshold 


cles.) Consequently, the events consisting in the 
interruption of the rumor in one or another chan- 
nel are independent. f 
The probability for all q channels to be inter- 
rupted is therefore equal to the product of the 


probabilities for each of the q channels to be 
interrupted: 


W (2) = IQ (It (8) 
Substituting formula (3) into (2), we obtain 
an equation for Q (x): 


Q (z) =1 — z + x [Q (x) (4) 


Note that the decisive factor that enables us 
to reduce the problem to algebraic equation (4) 
is the independence of channels. This property 
is an exclusive property of the Bethe lattice, 
so that the method we use is not successful when 
applied to usual lattices, although it is often 
employed to obtain an approximate solution. 

Let us analyze equation (4). It is meaningful 
for all æ within the interval OS z< 1. Re- 
writing it in terms of P (z) = 1 — Q (x), we find 
[(1—P (x)? £4 P (z)-—x=0 (5) 


One of the solutions of (5) is P (x) = 0 forall a, 
although for q > 4 equation (5) is nonlinear and 
has other solutions. Thus, P (1) = 1 is also a 
solution for 2 = 1, and this second solution is 
physically meaningful because the probability P 
must be unity and not zero if all circles are 
open. 


For q = 2 the solutions of equation (5) are 
easily found. There are two solutions in this 
case: P (z) =0 and P (x) = 2 — 14/2. If z> 
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> 1/2, the second solution has a physical mean- 
ing. If z< 1/2, the solution becomes negative 
and thus meaningless. 

Therefore, for q = 2 we have the following 
solution: 


0 for 0Kz < 1/2 


2—1/z for 1/2<r<1 (6) 


P (x)= 


on this case the percolation threshold v, equals 
A similar solution exists for all q > 1, although 
the percolation threshold ze is a function of q. 
In the general case it is possible to find x, and 
the form of P (x) for x close to To, assuming 
from the start that P (x) <1 (this is always 
true in the neighborhood of percolation thre- 
shold). ‘The term (1 — P)’ in equation (5) can be 
expanded via the binomial theorem: 


(1— Py 1—qP 4 E59 pr... (7) 


Each subsequent term here is much less than the 
preceding term because P <A. Therefore, we 
substitule (7) into equation (5), assuming that 
only the three terms written above are significant 
on the right-hand side of (7). This gives 


q (a—1) pe = 
ya? =qPr—P 


Assuming P Æ 0, we divide both sides of this 
equality by P and obtain 


_ (z—1/4):2 8 
P= (8) 
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Formula (8) coincides with (6) for q = 2. This 
solution vanishes at £ = 4/q, which implies that 
t = 1/g. Solution (8) is meaningful if 1> 4. 
æ > 1/q, and only if z is very close to 1/q. There- 
fore, we can set x = 1/q in the denominator of 
(8). Finally, we obtain 


P (a) = ih, (9) 


Formula (9) describes the probability P (x) in 
the vicinity of percolation threshold. 


Discussion 


If qg = 1, the function P (x) = 0 for all æ within 
the range 0 < z< 1, and equation (5) has no 
other solutions within this range. If q = 1 and 
v= 1, equation (5) holds for any value of P. 

If q = 1, the Bethe lattice transforms into a 
linear chain of lattice sites, And no matter how 
small number of black sites interrupts percola- 
tion through white in such a chain. It is therefore 
natural that whatever the value of æ within the 
range UKS z < 1, the interruption of the rumor 
propagation is unavoidable, that is, P (x) = 0; 
In a linear chain ig =f. 

We have shown in the preceding section that 
in the general case te = 1/q. This result could 
be predicted in advance, Each person in the 
problem under discussion transmits the rumor 
lo q of her (his) acquaintances, The mean number 
of people of category O among these acquaintances 
equals gz. Hence, each transfer of information on 
the average creates gx sources instead of one. 
The quantity qz is therefore its branching ratio. 


= 


Ch. 41. The Bethe Lattice 205 


In order for the process to continue, the branching 
ratio must be greater than unity. Hence, the 
critical concentration t is obtained from the 
condition gag = 1, that is, ze = 1/4. 

Recall that the condition necessary to sustain 
the uranium fission chain reaction is writlen 
in the same form. In fact, the process of rumor 
spreading is also a chain reaction and is des- 
cribed in the same terms as a nuclear explosion. 

It is of interest to compare the obtained value 
of ze with the results derived for lattice prob- 
lems in spaces wilh a greater dimensionality. 
Percolation threshold of the site problem was 
approximately computed for the so-called kyper- 
lattices. These are lattices of the same type as 
the square and simple cubic ones bul in a space 
with a greater number of dimensions. The coor- 
dination number z(the number of nearest neigh- 
bors) is found in such lattices from the formula 
z = 2d, where d is the dimensionality of space 
(it gives = 4 for d = 2, and z = ô ford = 3). 
The percolation threshold xs was computed for 
d = 4, 5, 6. It was found that the formula 


n= (1452) T (10) 


provides a good fit to the obtained results. 

The second term in parentheses can be neglect- 
ed if d is sufficiently large, so that t; = 
= 1/(z — 1). But recall that z — 1 = for the 
Bethe lattice. Indeed, only one bond enters 
each site of this lattice, and q bonds originale 
from il. 

We see thal percolation threshold of the Bethe 
lattice (£e = 1/9) is the same as that of a hyper- 
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lattice wilh a very large dimensionality. Hence, 
the Bethe lattice as if represents an infinite- 
dimensional space. 

The Bethe lattice is the only system for which 
it proved possible to find the exact form of P (x) 
at percolation threshold. It was found (see (9)) 
that in this case P (x) tends to zero linearly as 
z> te: P Œ (x — ze). We shall sce later that 
this is a specific property of the Bethe lattice as 
well as of all lattices of high dimensionality. 


Exercise 


1. Analyze the bond problem on the Bethe lat- 
tice. Consider all sites to be identical, and the 
bonds to he intact or broken. Let the fraction 
of intact bonds be x. Find the function P (x) 
defined as in the text. 


Chapter 12 


Structure of Infinite Clusters 


The Shkloyskii-de Gennes Model 


Now we shall consider, for the sake of definiteness, 
the site problem and assume that the concentra- 
tion of nonblocked sites is slightly above the 
threshold value, so that an infinite cluster exists. 
It consists of infinite strings (chains) of inter- 
connected sites. If all connected sites of an 
infinite cluster are Joined by straight segments, 
we obtain a set of intersecting broken lines (see 
Fig. 15 showing one such line). 
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By definition, the structure of an infinite 
cluster is its geometry on a scale much greater 
than the lattice period. On this seale, the sharp 
knees oceurring at individual lattice sites are 
not resolved lo the eye, and strings appear as 
smoothly bent curves. 

Figure 43 shows a small fragment of an infinite 
cluster. The cluster does not end at points A 


D 


Fig. 43. Fragment of an infinite cluster with dead-ends 


and B but stretches leftwards and rightwards to 
infinity. Now let us introduce the following clas- 
sification of points and lines in an infinite cluster: 
the elements of an infinite cluster belong to either 
ils backbone or dead-ends. 

A point is said to belong to the backbone of an 
infinite cluster if there are at least two paths that 
emerge from this point in opposite directions 
and lead to infinity. Point C in Fig. 43 is one of 
such points. Going to the left or to the right of 
this point, we can move infinitely far from it. 
If only one path leaving a point leads to infinity, 
this point belongs to a dead-end. For instance, 
ards from point D in Fig. 43 


only the motion upwa 
leads to infinity. The downward motion ends at 
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a cul-de-sac. Point D is thus said Lo lie on a dead- 
end. 

Let us imagine that all dead-ends are deleted; 
now we can try to discern the structure of the 
backbone of an infinite cluster. The simplest 
backbone model was suggested independently by 
the Soviet physicist B.I. Shklovskii and the 


ia 


Fig. 44. Backbone of an infinite cluster. 


French physicist P.G. de Gennes. In a plane 
problem, this model resembles a very large, old, 
and fairly worn-out fishnet. The net lost its 
regular periodicity, its cords have slack, some 
nots are broken, while others slipped off to wrong 
locations, but nevertheless it’s still a net “on 
the average” (Fig. 44). 

The characteristic linear size R of a cell of 
this net is called the correlation radius of an in- 
finite cluster. It grows drastically as we approach 
percolation threshold: 

l 
lz—ze |” g 


where J is a length that in order of magnitude 
equals the lattice period, and v is a positive 
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quant ily called the exponent of correlation radius. 
The network thus becomes less and less dense as 
percolation threshold is approached. 

A correlation radius approaching infinity is 
the general property of all critical phenomena. 


LoL 
1 
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Fig. 45. Correlation radius as a function of z. The graph 
shows the width 6 of the critical region for the L X L 
square (see the next section). 


The power-law behavior (1) is not a rigorously 
proved property; nevertheless, it forms the basis 
for the modern theories of critical phenomena 
and appears to be well supported by experimental 
data. 

A correlation radius remains meaningful at 
£< te, that is, below the threshold. In this 
region it describes the maximum size of finite 
clusters. If £ —> Te on the side of smaller values 
@< He). a correlation radius also tends to 
infinity as in (1). This means that when percola- 
tion threshold is approached from below, finite 


1414—0240 
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clusters grow infinitely and merge into an infinite 
cluster at £z = zte. The dependence of R on x 
thus has a shape schematically traced in Fig. 45. 

For three-dimensional problems the Shklovskii 
de Gennes model is formulated similarly. It 
resembles a badly damaged wire skeleton of a 
three-dimensional lattice whose mean cell size 
is given by formula (1). But it must be borne 
in mind that the numerical values of the expo- 
nents of correlation radius are different for plane 
and 3D problems. 

Let us consider now the corollaries derived from 
the network-like structure of infinite clusters. 


Role of the System’s Size 


We emphasized in Chapters 1, 2, and 3 that the 
concept of percolation threshold is truly meaning- 
ful only in an infinite system. In a finite system, 
percolation threshold varies among samples, that 
is, it is a random variable. However, the values 
assumed by this random variable fall, with very 
high probability, into a certain range of width 
ô (S°) called the critical region. As the number of 
sites in the system grows, the width of this region 
undergoes a power-law decrease (see formula (8) 
of Chapter 1), so that percolation threshold 
acquires a clear-cut sense as.f*—» oo and con- 
verts from a random variable into a certain 
quantity. 

This information was practically quoted, with 
no attempt to derive it, in the opening chapters 
of the book. The concept of correlation radius 
makes it possible to understand these concepts 
and to derive formula (8) of Chapter 4. 
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For definiteness, let us discuss the wire-mesh 
experiment on a grid L X L in size, the schema- 
lic of this experiment being shown in Fig. 1. 
Assume that numerous experiments using various 
random sequences of blocked sites have been 


Fig. 46. Square superposed on an infinite network, with 
R<L and z> to. 


conducted, and the result is a set of percolation 
threshold values. Recall that the configurations 
of blocked sites are very dissimilar in different 


experimental runs. 
It is convenient to pursue the following line of 


reasoning. Imagine an infinite network (a wire 
mesh), with a prescribed bed fraction x of non- 
blocked sites. Imagine now that an L X L square 
is superposed on the mesh at different places, and 
percolation is analyzed from the left- to the 
right-hand side of this square through the non- 
blocked sites that fell within the square (Fig. 46). 


14% 
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Placing the square over the different parts of 
the infinite wire mesh, it is possible to sample 
one by one the results of successive experiments 
with a finite wire mesh. 

Percolation sets in in an infinite network exact- 
ly at z = ze, but as we shall see below, this does 
not mean at all that percolation is always realized 
in an L X L square when z > ag. 

If £> ze, an infinite cluster exists in the 
infinite system. Let us picture its backbone as 
a fishnet (Fig. 44). The parameter extremely 
important hereafter is the ratio of correlation 
radius R to the square size L. First we assume 
that Z is much greater than R. Then (see Fig. 46) 
the square covers a large number of unit cells 
of the infinite cluster netting that ensures per- 
colation between the opposite sides of the square. 
These cells may be different in size, and the net- 
ting of the infinite cluster may contain large 
holes, but the probability for the cluster to con- 
tain a square-size hole is negligibly small if on 
the average the square contains very many cells. 
The following conclusion is therefore drawn: 


If z > zo, percolation threshold of a square 
cannot lie within the range of z satisfying (1) 
the strong inequality L > R (z). This region 
must lie above the threshold. 


According to formula (1), a correlation radius 
grows infinitely when z tends to Ze, so that at 
some value of x the radius inevitably becomes 
comparable with L. Now we cannot say anything 
definite about percolation across this square. 
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Everything depends on the specificities of the 
configuration of blocked sites within the square. 

Now let £ < Te, with the correlation radius 
much smaller than L. At £ < 2, the correlation 
radius is the maximum length of finite clusters. 
If R < L, no cluster exists long enough to con- 
nect the opposite sides of the square. Hence, 
another definite conclusion is made: 


Ifa< ze, percolation threshold of a square 
also cannot lie within the range of z satisfy- II 
ing the strong inequality L > R (z). This u) 
region must lie below the threshold. 


If £ < To but z is very close to ze, the cor- 
relation radius exceeds L. In this case we cannot 
say anything definite about percolation across 
the square. An infinitely large system contains 
finite clusters whose size is greater than L, but 
these clusters have holes of a similar size, so 
that everything depends on a specific configuration 
of blocked sites within the square. 

Now we can evaluate the size of the critical 
region in which the values of percolation thresh- 
old for an L X L square can lie. According to 
the corollaries (I) and (II), this region must obey 
the condition L < R. Figure 45 shows that this 
region is the narrower and squeezes up to per- 
colation threshold for an infinite system the 
closer, the larger L is. The width 6 of the region 
obeys the condition R (ô) = L. From formula 


(1) we find 1/8 = L, or 


ĝa 


g e) 
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Percolation thresholds of L x L squares are 
distributed uniformly within the critical region, 
that is, for | z — ze | <8 (see Fig. 5 where the 
distribution function of percolation thresholds 
is plotted). Nothing singles out the point x =a 
within this region. Indeed, this is the point at 
which percolation sets in in an infinile system. 
Bul it is impossible to establish whether per- 
colation is achieved or not if we work with finite- 
size squares. If L < R, itis impossible to estab- 
lish, by placing the Square over different regions 
of an infinite network, whether the network con- 
tains only finite clusters or whether these clusters 
have already merged and formed an infinite clus- 
ler. A study of percolation across a Jinite-size 
square only allows us to determine the width of the 
critical region. 

In this section we discussed only plane prob- 
lems. In fact, the arguments completely transfer 
to three-dimensional problems. Formula (2) de- 
lermines the width of the critical region in 3D 
problems equally well. A slight difference appears 
if the width § is written not in terms of the 
system's size L butof the total number of siles, f°. 
The point is that, = (L/a)", where a is the 

dimensionality of 


lattice period, and d is the 
space, Consequently, we have by (2) 


bf )= <= (3) 


where C is a numerical coefficient that ¢ 
found by equally simple re 
(d = 2), formula (3) coin 
of Chapter 4. It was by 
that the exponent 


annot be 
asoning. In plane cases 
cides with formula (8) 
means of this formula 
of correlation radius of the 
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plane problem was established for the first time 
from the dependence & (.f*) determined numeri- 
cally on a computer. It was found that v, = 1.33. 
(Here and below the subscript 2 indicates that 
we deal with the exponent of a two-dimensional 
system.) The exponent v of three-dimensional 
problems is different: Va = 0.8-0.9. (The sub- 
script 3 denotes that the exponent refers to 3D 


problems.) 


Electric Conduction 
Near Percolation Threshold 


Let us consider specifically two- or three-dimen- 
sional networks with blocked sites. As we said 
in the opening chapters, such networks are 
electrically conductive at x > Ze, but at the 
percolation threshold £e their conductivity va- 
nishes. Experimental data and the data obtained 
by computer simulation show that the specific 
electric conductivity of networks tends to zero 


in a power-law fashion, 
(4) 


= wi it 
a (x) = % (© — te) 
where the factor o equals, by the order of mag- 
nitude, the specific electric conductivity of the 


network without blocked sites. The quantity t 
is called the critical exponent of electric conduc- 
tion; it constitutes the subject of a very careful 
analysis, mostly by numerical computer simula- 
lion (e.g. one of the more recent computations 
involved a square lattice with 800 x 800 sites). 
Tt was established that t; = 1.3 in two-dimen- 
sional networks, and ts = 1.6-4.7 in three-di- 


mensional networks, 
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The network model of an infinite cluster makes 
it possible to derive formula (4) and to relate 
its exponent ¢ to the exponent of correlation 
radius. Electric current flows only through an 
infinite cluster, namely, through its backbone. 
No current flows in dead-ends that couple to the 
backbone only at one end. If electric current 
were sufficiently strong, so as to heat the wire lo 
a glow, the backbone of an infinite cluster would 
be visually discernable in the dark, as a rate of 
illuminated channels against a dark background. 
Far from threshold the whole network glows more 
or less uniformly, with the distance between glow- 
ing channels increasing as we approach the 

the glow fades out when the 
: the current ceases to flow. 


(see Fig. 46). The 
ance of a square with 


\ l ver of wires intersect- 
ing this square is 1/R, where R is the spacing 


between wires, given by formula (1). Let us denote 
the resistance of one wire with unit length by po 


The circuit is the parallel Connection of all wires. 
Therefore, the resislivily is 


(5) 
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and the electric conductivity is 
o=pt=p,'Rt (6) 


Substituting (1), we obtain 
0 = Gg (€ — 2o)” (7) 
where o, = ppl. 

In the three-dimensional cases we need to 
calculate the resistivity of a wire skeleton that 
models, for instance, the simple cubic lattice 
with period R (lattice type affects only the 
numerical coefficient). The resistivity equals the 
resistance of a cubic cell with unit-length edge. 
The number of wires connected in parallel pass- 
ing through a face of such a cube equals 1/R°. 
The resistivity therefore equals 


R= iy = pR? (8) 
and the electric conductivity is 
o= p R= 45 (t— =e)” (9) 


where o, = pō l™. 

In order to avoid confusion, note that the 
electric conductivity € has different dimension- 
alities in the two- and three-dimensional cases: 
it is measured in ohm~ in the 2D case, and in 
ohm-!-cm7! in the 3D case. 

The coefficients Ca and Gs correspond, by, the 
order of magnitude, to the electric conductivities 
of two- and three-dimensional networks with no 


blocked sites. Indeed, as we see from formulas 
6) and (9), the electric conductivity © (2) reduces 
to o, OF Ogy respectively, if R = L that is, if the 
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network of an infinite cluster coincides with the 
initial network on which the problem has been 
formulated. The factor oy in formula (4) is thus 
equal to o, in the 2D case, and to og in the 3D 
case. 

A comparison of formulas (6) and (7) with 
formula (4) shows that ¿ = v in the 2D case, and 
t = 2v in the 3D case. Using ve = 1.3 and 


Fig. 47. Backbone of an infinite 


f cluster, with “sinuosity” 
taken into account. 


va = 0.8-0.9, we find t, = 1.3 and t = 1.6-1.8, 
quite close to the data given above. This agree- 


ment strongly supports the Shklovskii-de Gennes 
model. 


The backbone model of 
be generalized as follows. Imagine that the wires 
forming the backbone are very sinuous (Fig. 47). 
The distance between their intersection points 
is, as before, R (x) and is given by formula (1). 
However, if the wire segment between two inter- 


an infinite cluster can 
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section points is straightened out, its length 
will prove considerably greater than R. We 
denote this length by £ and write it in the form 
p l 
(£— te)? 


(2> te) 


where >v. As 2p, the ratio £/R = 
= (£ — £) tends to infinity. 

The introduction of the length £ is justified 
only if £ > R, that is, € > v. However, it has 
been proved that ¢ = 1. Hence, the suggested 
generalization of the Shklovskii-de Gennes model 
(in fact, this generalization has been carried by 
the authors of the model) is meaningful only if 
v< 1. Since v = 1.3 in the 2D case, no “sinuosi- 
ty” is expected there. But in the 3D case v < 1, 
and there is a good reason to believe that “sinuosi- 
ty” indeed characterizes the backbone of an infi- 


nile cluster. 


Exercise 

1. Express the critical exponent £ of electric 
conductivity in terms of the exponents € and v 
in the 3D case, making use of the generalized 


model. 


Function P (x) Near Percolation Threshold. 
Role Played by Dead-Ends 
The function P (x), which is the fraction of sites 


belonging to an infinite cluster, vanishes at « = 
= £e, together with electric conduction. An 


analysis demonstrated that in the immediate 
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vicinity of the threshold this function varies as 
P (x) = D (x — zo)? (10) 


where D is a numerical coefficient of the order 
of unity, and B is another critical exponent. It 
has been found that B, = 0.14 (2D problems), 
and Ba = 0.4 (3D problems). These results were 
mostly obtained in computer simulation. 

All sites of an infinite cluster, both those be- 
longing to the backbone and to dead-ends, con- 
tribute to P (x). The infinite cluster model makes 
it possible to find which of these subsets has a 
greater number of sites. First assume that dead- 
ends are completely absent and calculate the 
contribution of the backbone of the infinite clus- 
ter to P (z). 

The number of sites belonging to the backbone 
per cell of an infinite cluster in the 2D case is of 
the order of R/a, where a is the lattice period 
(as in the preceding section, this is also an 
evaluation that is not meant to find numerical 
coefficients). The cell area is about R?, and 
hence, the total number of all sites in a cell is 
about NR?/a?. Therefore, the fraction of sites 
belonging Lo the backbone of an infinite cluster is 


Pry (2) ~~ (@—a)” (11) 


Here the symbol ~ stands for equality to within 
an order of magnitude (neglecting numerical 
coefficients of the order of unity). 

In the 3D case the number of sites of the back- 
bone per each cell of an infinite cluster is also 
Ria, but the total number of sites per cell is 
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a \2 
Pal) ~ (Fe) ~ (ee) (12) 

Comparing formulas (11) and (12) with for- 
mulas (7) and (9), we obtain that the fraction 
of sites belonging to the backbone of an infinite 
cluster coincides, to within an order of magni- 
tude, wilh the function o (x)/o = (£ — ap)'. 

Comparing (14) and (12) with formula (10), we 
notice that 


Pow (2) 
P (2) 


in the 2D case, and 


P = 
a ~ (t—25)"" Bs 


in the 3D case. 
= 4.3 and vy + 0.9. Hence, 


Recall that vz 
va — Pa = 1.2 and 2v3 — Ba ~ 1.4, Therefore, 


the ratio Pyy (x)/P (x) rapidly tends to zero 
as £ —> To both in the 2D and 3D cases. This 
means that the sites forming the backbone of an 
infinite cluster comprise an infinitesimal fraction 
of the total number of sites belonging to the 
infinite cluster. The “mass” of an infinite cluster 
predominantly resides in its dead-ends and is com- 
pletely useless from the standpoint of electric 
conduction. Consequently, in the vicinity of 
percolation threshold o (x)/0o <P (x) (see 
Fig. 10). However, the spontaneous magnetiza- 
tion of a doped ferromagnetic in the neighbor- 
hood of percolation threshold is determined 
precisely by these dead-ends (see Chapter 3). 
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Universality of Critical Exponents 


We have encountered three critical exponents 
that describe the behavior of various quantities 
in the neighborhood of percolation threshold, 
i.e. v, t, and B. This behavior is also referred to 
as critical because the functions R (z), o (x) and 
P (x) have singularities al x Te. For instance, 
R (x) tends to infinity, while P (x) has a dis- 
continuous first derivative (al x). The derivative 
is zero on the left of Te and is infinitely large on 
its right. In the case of o (x), it is a second deriv- 
ative that is discontinuous. Many other quan- 
tities with critical behavior are known in per- 
colation theory, and correspondingly, many other 
critical exponents. 

For each of the above-discussed critical expo- 
nents we gave two values: one for 2D and another 
for 3D problems. However, the list of relevant 
2D and 3D problems would be extremely long. 
For instance, we know of three-dimensional sile-, 
bond-, and sphere problems, the problem of per- 
colation level in a random potential, and many 
others. What, then, were the exponents we dis- 
cussed so far? Now we come to what seems to be 
the most exciting feature of percolation theory. 
According to modern concepts, the critical expo- 
nents are identical for all problems in the space of 
a given dimensionality. (The only exceplion among 
the problems mentioned in this book is the prob- 
lem of directed percolation.) The statement of 
universality of critical exponents is rather a 
convention than a rigorously proved proposition. 
Nevertheless, numerous tests of this statement, 
carried out on computers, failed to refute it. 
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What are the physical reasons for the univer- 
sality of exponents? Presumably, the exponents 
are determined by the structure of clusters in the 
vicinily of percolation threshold. The properties 
that play the principal role here are the geometri- 
cal properties of clusters that are fell at large 
distances (of the order of correlation radius). 
In the neighborhood of the threshold these distan- 
ces hecome much greater than the lattice period (in 
lattice problems) or the radius of spheres (in sphere 
problems). Consequently, the geometry of clusters 
is independent of the type of lattice on which 
the problem is formulated. The structure of large 
clusters will remain unaltered even if a problem 
is given not on a lattice but on sites randomly 
distributed in space. But, of course, the dimen- 
sionality of space strongly affects the geometry 
of clusters because, for instance, it is much easier 
to ensure “bypassing” of curves in a three-dimen- 
sional space than in two dimensions. For these 
reasons, critical exponents do not depend on 
problem type, but do depend on space dimen- 
sionality. 

It is interesting to note 
change with increasing dime 
only up to dimensionality 6. Beginning with 
d = 6, exponents remain constant, and the expo- 
nent B = 1 as in the Bethe lattice. When d > 6, 
the critical exponent problem considerably sim- 
plifies and allows an exact solution. 

Critical exponents thus possess a certain uni- 
versality in contrast to percolation thresholds 
which strongly depend on the type of problem. 
This leads to a simple conclusion. If the results 
of a physical experiment are treated in terms 


that critical exponents 
nsionality of space 
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of percolation theory, and the microscopic struc- 
ture of the investigated system is not quile clear, 
the first characteristics to be compared with the 
theory must be critical exponents because they 
are almost invariant. This is the slralegy one 
uses in analyzing experimental data on electric 
conduction in heterogeneous materials (see Chap- 
ter 9). 

Percolation theory borrowed the idea of uni- 
versality of critical exponents from the theory of 
second-order phase transitions (e.g. among second- 
order transitions are the transition of a metal 
from the superconducting to the normal state and 
the transition of a ferromagnetic to the nonferro- 
magnetic state, both brought about by increasing 
the temperature of a sample). As in the vicinity 
of percolation threshold, large regions differing 
in properties from one another form close to the 
point of a second-order phase transition. The 
difference lies in that the boundaries of these 
regions are not “frozen” as they are in percolation 
theory, but change with time owing to thermal 
motion. The size of such regions is also called 
the correlation radius and also given by for- 
mula (4). 

The theory of phase transition supplied another 
important idea as well: the sealing hypothesis. 
We shall formulate it in terms of percolation 
theory. Assume that different parts of a net- 
work with some blocked sites have been photo- 
graphed at two values of z, namely, x, and zg 
Both these values are on one side of Zo and are 
close to it. Let, for instance, x, — Te > Ta — 
— zt, >0. At z = z, the correlation radius is 
R, = l (z, — ze)”, and at z = z it is Ry = 
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= l (£, — x,)~%. These conditions imply that 
R, < R,. The scaling states: if the photographs 
taken at x = z, are magnified by a factor R,/R,, 


they will not differ, “on the average”, from the 


photographs taken at 2 = Tə. The photographs 
ined that only 


are assumed to be so coarse-gral 
large blocks show up, and individual sites and 
bonds are not resolved. The scaling hypothesis 
thus states that the large-scale geometry of the 


system is transformed on approaching percolation 
ith all linear 


threshold in a self-similar manner, w1 
dimensions growing in proportion to the corre- 


lation radius. 
k Note that the Shklovskii-de Gennes model 
satisfies the scaling hypothesis, but that this 
hypothesis goes far beyond the model. It is neither 
restricted to the backbone of an infinite cluster 
nor does it assume that the network is subdivided 
into its backbone and dead-ends. 

The mathematical formulation of scaling hypo- 
thesis makes it possible to find a relation between 
the critical exponents v, B and a third exponent 


that we chose not to introduce here. Calculations 
show that this relation holds quite well. 
duced by 


Scaling concepts that were first intro 
the Soviet physicists . Z. Patashinsky and 
V. L. Pokrovsky and by the American physicist 
L. Kadanoff form the pasis of the modern 
theory of phase transitions and of percolation 
theory. 


Scaling concepts prought abo j 
new mathematical techniques for calculating 


critical exponents. A dramatic progress has been 
achieved in these methods in recent decade, so 
that certain perfection has been reached. At pre- 


ut a number of 


15-0240 


226 Part III. Quantities Near Percolation Threshold 


sent, these methods help to calculate almost all 
imporlant exponents. However, computations 
are so complicated that they could not be dis- 
cussed in this book. 


Chapler 13 
Hopping Electric Conduction 


The preceding chapters gave a detailed picture 
of how percolation theory is applied to calculate 
the electric conductivity of systems consisting 
of randomly coupled elements. Each element 
was either conducting or insulating, but the 
resistances of all conducting elements were as- 
sumed equal. Examples of such systems are nel- 
works with broken bonds or blocked sites, mix 
tures of metallic and dielectric spheres, and so 
forth. 

Let us turn now lo a different class of systems 
which are also composed of different elements, 
but the resistances of these elements can take 
on any values from extremely small to enorm- 
ously large. It was found that the resistance 
of systems composed of a very large number of 
such elements can also be calculated in terms 
of percolation theory. The theory of hopping 
electric conduction of semiconductors, based on 
percolation ideology, was constructed in 1974. 
In this chapter we make acquaintance with hop- 


ping electric conduction and present its mathe- 
matical description, 
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Mechanism of Hopping Conduction 


Consider an extrinsic semiconductor doped wilh, 
say, donor impurities. The dopant concentration 
will be considered small in comparison with the 
critical concentration Ne at which the Mott 
transition (metal-insulator transition, see Chap- 
ler 8) lakes place. In these conditions, outer 
electron shells of neighboring impurity atoms 
overlap very slightly. Each donor exists therefore 
as an isolated hydrogen-like atom whose outer 
electron is at a distance about aj from the nucleus 
and whose bonding energy is about E} (see Chap- 


ter 8). Assume that the temperature of the semi- 


conductor is so low that the thermal energy of 
vibrating atoms is not sufficient to ionize a donor. 
What could be the mechanism of electric con- 
duction in this situation? 

Imagine that some donors lost their outer elec- 
tron. Usually this results from impurity com- 
pensation (see Chapter 10). If the semiconductor 
contains both donor and acceptor impurities, 
each acceptor accommodates one electron from 
a donor. If the number of acceptors is smaller 
than that of donors, a fraction of donors will have 
retained an outer electron, while the other frac- 
tion will have lost their electrons and become 
positively charged. 

The mechanism of hopping electric conduction 
consists in the “hops” of an electron from one 
donor to another which before this hop had no 
outer electron. 

Now we are going 
outer donor electron 
depends only slightly on t 


15% 


to consider the case of an 
whose potential energy 
he spatial position 
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of the donor, that is, on the specificities of the 
configuration of charged impurities surrounding 
this donor. This case corresponds to a low con- 
centration of impurities. In conditions typical 
for hopping conduction, the spread in the ener- 
gies of outer electrons belonging to different 
donors is roughly 0.4 of the bonding energy £$- 

In this situation, the attraction of an electron 
to the donor to which it belongs at the initial 
moment is the main obstacle to a “jump” from 
one donor to another. From the standpoint of 
classical mechanics, an electron can be trans- 
ferred from the outer orbit of one donor lo the 
outer orbit of another if some work is done 
against attractive forces because it is necessary 
to remove the electron to one half of the distance 
between the donors. After this point the electron 
is attracted to the second donor. If donors are 
spread very thinly, this work is almost equal to 
the work required to remove the electron belong- 
ing to an isolated donor to infinity, that is, equal 
to the bonding energy Lx. 

What then is the advantage that hopping elec- 
tric conduction holds over the electric conduc- 
tion by free electrons? The point is that hopping 
electric conduction is a quantum phenomenon. 
Quantum mechanics allows an electron to pass 
from one donor to another without rising to a 
free state and without borrowing energy from the 
thermal motion of atoms. This is the tunnelling 
transition. Tunnelling transitions obey the law 
of energy conservation. The law imposes the con- 
straint that the electron energy be equal in the 
initial and final states, Consequently, if the 
energies of the electrons at the first and second 
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donors differ owing to the potentials of the sur- 
rounding impurities, a deficit of energy has cer- 
tainly to be borrowed, and a surplus of energy 
has to be dissipated. But this energy is ten 
times less than Æ$. As a result, hopping electric 
conduction defeats the electric conduction by free 
electrons at very low temperatures. 


Resistor Network 
Tt must not be overlooked, however, that a tun- 
nelling transition is a very low-probability event 
in conditions we discuss. As was said in Chapter 8, 
the probability of finding an outer electron of 
a donor at a distance r from the donor’s nuc- 
leus decreases with r as exp (—2r/aj). If two 
donors Z and 2 are at a distance Tyo" from 
each other, the probability of finding the outer 
electron of donor Z near the nucleus of donor 2 
equals exp (—2r,./a$). This exponential enters 
the probability of tunnelling transition. As fol- 
lows from Chapter 8, at a donor concentration 
much less than the critical concentration No 
of the metal-insulator transition, the mean dis- 
tance between donors is much greater than aj, 
so that the quantity — exp (—2ry./a8) = 
= 1/exp (27;2/ap) is, as a rule, very small. 
Nevertheless, tunnelling transitions do take 
place time and again between neighboring donors. 
If an electric field is applied to a semiconductor, 
these transitions will be more frequent in the 
direction against the field (along the applied 
force) than along the field. The result is an elec- 
tric current proportional to the electric field 
strength. This is the essence of the phenomenon 


we call hopping conduction. 
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The model used to calculate the resistivity of 
a semiconductor is the so-called “resistor net- 
work”. This model is formulated not in terms 
of atoms and tunnelling transitions but of con- 
ventional electric circuit resistors. Imagine that 
a resistor is connected between cach pair of 
donors. The donors themselves can be pictured 
as tiny metallic balls to which the wires from 
numerous resistors are soldered. The second end 
of each of these resistors is soldered to another 
ball. The result should resemble an irregular 
three-dimensional resistor network. Obviously, 
an intention to assemble an actual network model 
would require that the scale of the system be 
substantially magnified. Indeed, the mean dis- 
tance between donors is of the order of 10-5 cm. 

The resistors connecting two donors must be 
found by calculating the tunnel current that 
flows between these donors at a given electric 
field. This means solving a quantum-mechanical 
problem that will not be given here. Note only 
that, in view of the above description, the elec- 
tric current produced by tunnelling transitions 
between donors in a given field is the smaller, 
the greater the separation between these donors is. 

Accordingly, the resistance 9 connecting do- 
nors at a distance r from each other can be written 
in the form 
A (r) = Ry exp (2r/ag) (1) 
where #, is identical for 
of the order of 4 ohm). The problem is to find 
the resistivity of a system composed of a tre- 


mendously large number of donors (1019-101) 
with random patial distribution, 


all resistances (it can be 
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Properties of Resistor Network 


r : 

The main feature of the model we discuss is that 
the resistances described by formula (1) are 
spread over a fantastically large range. The 
mean distance Taon between donors is related to 
the donorconcentration Naon by (4/3) mrdonV aon = 
=, which means that on the average one donor 
is nested in a sphere of radius Taon: Typically, 
hopping electric conduction is observed in con- 
ditions in which Taon exceeds ağ by a factor of 
6 to 12. 

Assume, for instance, that Tygon = 10a}. Then 
the resistance connected between donors at a 
distance 1.5rgon is greater than that between 
donors at a distance qon by a factor 
exp (Taon/ah) =e = BDA OX, 

Donors at a distance Taon occur almost as fre- 
as those at a distance 1.5rgon: Conse- 
ge in separation b9- 
tremendous var ‘ation 


quently 
quently, a negligible chan, 
tween donors results in a 
of resistance connected between them. 

In principle, the resistor network model stip- 
ulates that each pair of donors is connected 
through a resistor. However, the resistances 
between remote donors are so large that they 
can undoubtedly be neglected. As a rule, the 
same two donors are connected through a chain 
of resistors connecting nearest-neighbor donors. 
Although the length of this chain is greater than 
the shortest distance between remote donors, the 
stance is much less than one resistance 


chain resis 

connecting these two donors. Such are the prop- 
erties of the exponential function: if 2 >1 
and 2, > 1, then exitXs = gg D> ex: p eis 
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Consequently, it is amply sufficient to retain 
in the resistor network only those resistors that 
connect each donor with four or five of its nearest 
neighbors. 


The Sphere Problem Revisited 


Our next step is to calculate the resistances, We 
Suggest the following line of reasoning. Let us 
disconnect all resistors Supposed to join the 
balls that stand for donors and start connecting 
them in a predetermined sequence. 

First we solder in the resistances that connect 
the donors Separated by a distance less than a 
certain length r’. This means that we include 
the lowest resistances smaller than #' = 
= Ro exp (2r'/a¥). If r is much smaller than 
the mean Separation Taon between donors, then 
at this stage we connect only the rare donors 
separated by anomalously short distance. As a 
rule, these resistances are not interconnected and 
Cannot sustain the flow of current through the 
System. 

Let us gradually increase r’, each time adding 
new groups of resistances, Beginning with a 
certain value 7’ = Te, the resistors form an in- 
finite cluster, At r’ = ra the system is electri- 
cally conductive and has a finite resistivity. 

The problem of finding r, is nothing less than 
the sphere problem (see Chapter 7). Indeed, let 
us construct a sphere of radius r’ around each 
donor. Then resistances will be connected only 
between those donors one of which lies inside 
the sphere Constructed around the other. Accord- 
ing to the results of Chapter 7, such donors 
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join into an infinite cluster when the condition 
By = (4/3) tNaonT2= 2.7 + 0.4 


is satisfied, whence 
Te = (0.865 =- 0.015) Naor? = 1.39r aon (2) 


Calculation of Resistivity 


As was mentioned in Chapter 12, the backbone 
of an infinite cluster can be pictured as a three- 
dimensional wire network (see Fig. 44) with knots 
spaced by distances of the order of the correla- 
tion radius. In the present problem each wire 
must be regarded as composed of numerous balls 
connected with resistors, and expression (1) of 
Chapter 12 must be rewritten in the form 
Ti ad 

R()=ra (qrr) (3) 

Expression (3) turns to infinity at percolation 
threshold, with the critical exponent v which, 
by virtue of universality, is the same as in other 
three-dimensional problems of percolation theory 
(v = 0.8-0.9). When r’ — re is of the order of 
we are far from the threshold, and the 
on radius becomes of the order of Fion: 

When r’ = re, the separation between network 
sites in an infinite cluster grows to infinity. At 
precisely this point, the density of the infinite 
cluster is still zero. However, if T S>7',,, the 
cluster forms channels going across the whole 
system and resulting in a finite resistivity. 

Let us continue with the procedure of adding 
resistances. If we switch on resistances between 
donors spaced by distances from re to Te + gap, 


Taon» j 
correlati 
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where g is a number smaller than unity (e.g. 
g = 0.2), the maximum one among the con- 
nected resistances will remain practically unal- 
tered because exp [(2r,/a%) + 2g] œ% exp (2r¢/ap)- 
On the other hand, an infinite cluster represents 
a network with site-to-site separation of about 


Re =Tqon (e i (4) 


gay, 


In the theory of hopping conduction this network 
is referred to as critical. It is through this network 
that the electric current flows. 

Indeed, resistivity will not be appreciably 
changed by any further increase of r’, that is, 
by adding the remaining resistances. As we find 
from formula (1), resistances for r — r, > ab 
are very large in comparison with resistances for 
r=ro. Infact, practically no current flows through 
these resistances since it “prefers” the critical 
network whose maximum resistances are 


Pmax = Ro exp (—2r,/ag) 


The last step to accomplish is to calculate the 
resistivity of a critical network. By analogy to 
what we did in Chapter 12, this network will be 
represented by a wire skeleton modelling a 
simple cubic lattice (see Fig. 12) with period Re. 

Each wire connecting two neighboring sites 
of the lattice consists of a large number of resist- 
ances given by formula (1) and connected in 
series. As individual resistances differ from one 
another enormously, the resistance of one wire 
must be set equal to the maximum one of the 
resistances it comprises, thal is, lo mas. The 
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cubic unit cell of a simple cubic lattice (see 
Fig. 27) consists of four wires connected in paral- 
lel, each wire belonging at the same time to 
four other cubic unit cells, so that there is one 
wire with resistance max per each unit cell. 

In order to find resistivity, it is necessary to 
multiply the total resistance may of a cubic 
unit cell by the area of ils face (R) and divide 
it by the edge length (Re). This gives 


P = FP maxr Ro = Ao: Re exp (2r,/ag) (5) 


We thus deal with a cubic unit cell as if it 
were filled by a homogeneous material with re- 
sistivity p. This is the meaning imparted to the 
concept of resistivity of a system that is in fact 
strongly inhomogencous. 


Discussion of the Result 


It must be borne in mind that different elements 
of formula (5) were obtained with different de- 
grees of accuracy. The factor Re, called the 
preexponential factor, was found to within a 
numerical coefficient. Indeed, we do not really 
know the quantity g in formula (4) for Re. Be- 
sides, the replacement of the network of an infinite 
cluster by a periodic lattice may well lead to an 
error in the numerical coefficient. It must be 
mentioned that no theory has been constructed 
yet that could reliably determine the numerical 
coefficient in the preexponential factor. 

As for the quantity in the exponent of the 
exponential expression (5), it is known to a good 
accuracy. In the true spirit of the derivation of 
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formula (5), re in the exponent could be replaced 
by re + apg. p 

Then the sum 2r,/ag + 2g would appear in- 
stead of 2r,/aġ. The second term of the sum ac- 
tually characterizes the uncertainty in our knowl- 
edge of the exponent. As has been mentioned, 
2r,/az > 1 at low donor concentrations, so that 
the relative error is small. In fact, it reflects the 
uncertainty in the numerical coefficient of the 
preexponential factor. 

The dependence of resistivity p on a donor 
concentration is determined solely by the expo- 
nent of the exponential. By virtue of (2) and (5), 


2-0.86 

Inp=In.%, Ret Tie (6) 

The first term on the right-hand side of (6) 
depends on Ngon much weaker than the second 
term, and thus can be considered constant. A com- 
parison between formula (6) and experimental 
data, carried out for numerous semiconductors, 
demonstrated that this formula gives a very 
good fit to the Ing vs Naon dependence; this 
agreement was an important achievement of the 
theory described above. 


The method discussed in this ch 
successfully applied to determining the lempera- 
ture dependence of hopping electric conduction, 
and to its dependence on the external magnetic 
field and on a number of other parameters, Fur- 
thermore, this method holds for any inhomogene- 


ous system whose resistance varies in a wide range 
as a function of coordinates, 


apler was 


Ch. 44. Final Remarks 237 
Chapter 14 
Final Remarks 


In this last chapter we shall briefly discuss some 
applications of percolation theory that have not 
been covered in the preceding chapters but never- 
theless appear to us as quite interesting, and we 
shall also make a summary, in an attempt to point 
at the features common for all the problems 
discussed in the book. 


Some Applications 


Flow of liquid through a maze. This problem is 
closely linked to that of gas penetration into 
carbon gas masks that served as the starting 
point for percolation theory. Imagine a porous 
body through which a liquid is forced under 
pressure. The liquid does not wet the material 
of the body so that capillary forces resist the 
penetration of the liquid into pores. 

A maze of capillaries permeates the body, but 
the capillary diameters are vastly different. They 
are wide in some places, and very narrow in 
others. 

First imagine that the body has a single cylin- 
drical capillary. At a certain pressure produced 
by a piston 2 a liquid is on the left of the body 7, 
and at atmospheric pressure air is on its right 
(Fig. 48). The body is fixed and cannot move. 
Because of surface tension, a nonwetting liquid 
penetrates the capillary only ifits pressure exceeds 
atmospheric pressure by the quantity 20/R, where 
R is the capillary radius, and o is the surface 
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tension. At lower pressures, the liquid forms a 
convex meniscus but cannot pass through the 
capillary. ; a 
Now imagine that the capillaries within the 
body Z are of different radii. Some of them (the 
widest) let the liquid flow through at a given 
pressure, but others (narrower ones) do nol. As 
pressure increases, the number of permeable 
capillaries grows. Al small pressures, when only 


Fig. 48. 


the widest capillaries are accessible t 
it cannot penetrate the bod 
surface layer. However, at a certain critical 
pressure, the permeable capillaries form a system 
that permeates the whole body. Beginning with 
this pressure, the liquid can be Squeezed through 
the body. 

Percolation theory calculates this pressure as 
well as some other characteristics of the process, 
which constitutes a very important practical 
application of the theory. 

Formation of polymer gels. A polymer consists 
of a very large number of elementary units (mo- 
nomers). Monomers in a solution may bind to 
one another and form a complex three-dimen- 
sional network that permeates the whole system. 


o the liquid, 
y deeper than its sub- 
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The result is a gel, thal is, a solid-like medium 
resembling jelly. 

A model is known that adequately describes 
the formation of a gel. It is essentially a problem 
of percolation theory in which while siles are 
the molecules of a monomer, and black sites are 
the molecules of a solvent. The bonds between 
white sites form with a temperature-dependent 
probability. A gel appears when an infinite cluster 
forms out of connected white sites. 

This problem of percolation theory is said to 
be mixed because both sites and bonds are ran- 
dom elements here. Let the probability for a site 
to be white be x, (it equals the concentration of 
monomer molecules), and the probability for 
a given bond to be intact be 9. We want to find 
the range of values of a, and z, in which an in- 
finite cluster of interconnected white sites is 
formed. 

By definition, 2, and £, vary within the inter- 
val from zero to unity. 

Ifa, = 4, that is, if no bonds are broken, an 
infinite cluster exists for all z} in the range 
zs <z <1, where z, is percolation threshold 
of the site problem. If sz; = 4, that is, if all 
sites are while, the condition under which an 
infinite cluster is formed is £p < Ta <1, where 
xp is percolation threshold of the bond problem. 

The square in Fig. 49 is the domain of the 
variables 2, and xy. The solid curve is the graph 
of the function amin (2) that describes the bound- 
ary of the domain in which an infinite cluster 
exists. The function Tmin (£2) gives the minimum 
value of 2, foreach value of x, within the interval 
£p <a, < 1, al which the infinite cluster exists. 
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It is readily understood that Tmin (1) = zs, and 
Zmin (ty) = 1. The domain of the infinite cluster 
in Fig. 49 is the hatched area. 


0 Xp Loy 


Fig. 49. 


If the probability T is known for each value 
of temperature, the function Zin (Cy) allows 
us to determine the range of temperatures and 
monomer concentrations in which a gel is formed. 


What Is Percolation Theory, After All? 


So far we h defining per- 
colation theory, This would be far from simple. 


age of the book, we 
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in awell-defined manner (e.g. by arandom-number 
generator with known properties). 

The common element of percolation theory 
problems is that all of them have identical geo- 
metry of bonded elements in the vicinity of 
percolation threshold. In order lo recognize this, 
it is necessary to disregard the small-scale 
structure depending on the type of bonding and 
the properties of elements, and monitor only the 
conneclivily of large blocks. The universal large- 
scale geomelry dictates the universal properties 
of physical quantities that depend on the struc- 
ture of large clusters. This is the feature that 
uniles percolation theory problems that look so 
dissimilar at first glance. 


16-0240 


Answers and Solutions 


Chapter 1 


1. The general rule for calculating the mean value 
states that each possible value of a random variable must 
be multiplied by the probability of this value, and the 
products be summed up. The probability for each of the 
faces of the cube to be on top equals 1/6. Therefore, 


a = (1/6) + (1/6)-2 + (1/6)-3 4. (1/6) -4 +. (1/6) +5 
-+ (1/6) -6 = 21/6 

2. The results of individual trials 
average value ze (4^) calcul 
mains the same because on the average the rightward mo- 
tion realizes with the same probability as the downward 
motion. Correspondingly, Te = lim ze (A) will not 

N= 

change cither. 

3. Let us denote the threshold values recorded in the 

$ , , 

ith trial by zi and z;, where xi corresponds to the new 
definition, and 2; to the old definition (rightward percola- 
lion). If a decrease in x first interrupted the rightward 
percolation and then the downward percolation, then 
vi = 2;. But if the sequence was reverse, then zi > ti 


t 
The average value of percolation threshold of the new 
definition is 


change, but the 
ated over many trials re- 


zitt +... +z? 
X(N) = ih . Hro 


and that of the old definition is 


(J) = Sette sta 


The total number of trials, Q, is considered very large in 
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both formulas. However, situations with zi >r; will 
inevitably realize if the number of trials is large. For this 


reason, z? (.4°) > te (4°). Nevertheless, 


lim ze( f )= lim Telf )= Te 


A ~o M+ 


The thing is that in an infinite system percolation 
threshold is not a random variable but a certain quantity, 
not varying from one trial to another. At the same time, 
differences between percolation thresholds for different 
directions represent a random phenomenon. From the 
standpoint of percolation, all directions are on the 
average equivalent. Hence, ze is independent of 
direction. 

4. The solution is quite similar to that of Problem 3. 

Let us denote the threshold values obtained in the 
trial by zi and z; with x} corresponding to the new 
definition of threshold, and z; to the old definition. It can 
be readily proved that zi < z; Arguments similar to 
those of the preceding problem give 


lM TEA, but lim lN) 
Mo 


5. Formula (8) yields 6 = 0.01. This signifies that 
“typical” deviations from the mean value are +0.01. 
Hence, the last numeral in the estimate 0.59 reported by 
Watson and Leath was very likely in error. The proba- 
bility of error in the first numeral after the decimal point 
is much smaller. Since only one trial was run, the 
authors of the paper could only evaluate the error with 
which they determined percolation threshold on the 
utilized sequence of blocked sites (it proved to be 0.005), 
However, they could not say anything about the possible 
change in the result if the experiment were rerun with 
a different random sequence. Later studies in which many 
trials were run at one valueof 4^, and studies with large 
values of jf’, led to formula (8). These studies also de- 
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monstrated that even the second decimal place is correct 
in the number 0.59; to a large extent, this must be re- 
garded as pure “luck”. 


Chapter 2 


1. The fraction of blocked sites is 1— r= M'IN 
and the fraction of nonblocked sitèsisz = (W =M IUA 
If Q sites are chosen at random, Qz of them will be non- 
blocked, and Q (1 — z) will be blocked (the greater is Q, 
the better the accuracy is with which this relation holds). 
Consequently, the probability for a randomly selected 
site to be blocked is Q (1 — z)/Q = 1 — z, while the 
probability for it to be nonblocked equals Qz/Q = zr. 
The sum of the probabilities equals unity because a 
ap can be either blocked or nonblocked: 1 — Fe r= 


2. The probability of any sequence of three fixed num- 
bers equals 1/6-1/6+4/6 = 1/216. The number of differ- 
ent sequences satisfying the formulated conditions for 
the numbers 4, 2, 3 equals 6 (123, 213, 321, 231, 132, 312), 
and for the numbers 1, 2, 2 it is 3 (122, 212, 221). The 
probability for one of the possible sequences to realize 
equals the sum of these probabilities. Therefore, in the 
first case the sought probability is 6-1/216 = 1/36, and 
in the second it is 3-1/216 = 1/72. 

3. (0.8)3-(0.9)! = 0.336. 

4. The solution is left to the reader. 

5. The distribution function of a random variable a 
is constant within the interval (—1, 1) and equals zero 


ii 
7 
a 7 
~1 -3/4 -1/4 0 1 K 
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everywhere else (Fig. 50). The total area of the rectangle 
bounded by the curve f (y) (in the present case by a hori- 
zontal straight line), the abscissa axis, and the perpen- 
diculars erected at points —1 and 4 must equal unity: 
consequently, f (y) = 1/2 for —1 <y <1. By virtue 
of the gencral rule, the sought probability equals the 
area of the rectangle bounded by the line f (y), the abs- 
cissa axis, and the perpendiculars erected at points y = 
= —3/4 and y = —1/4 (the rectangle is shaded in Fig. 50). 
The probability equals 
[—1/4 — (—3/4)] -1/2 = 1/4 

6. The variable y assumes all values from —co to 4-00. 
_ Consequently, we have to set A = —~, B = -+œ 
in formulas (3) and (4). According to formula (3), 


œ 


a= | uh pway 


-œ 


Substituting here formula (6), we find 


= ig 
<= a or j yexp (— 


The integrand contains an odd function. Substituting 
y = —t and comparing 


the result with the original for- 
mula, we note that a = —a, hence, a = 0. 
According to formula (4), the variance is 


o0 


ô= f wf gy (U) dy 


-%0 


Ti Sr j yžesp (— 


Substituting y = V2 ph we obtain 
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The integral over t equals 7/2, so that 62 = ô? 


Chapter 3 


1. By definition, P (1) =41. If z is nearly unity, 
the sites may not belong to an infinite cluster for two 
reasons: 

(i) They may contain nonmagnetic atoms. The frac- 
tion of such sites is 1 — z. S: 

(ii) Magnetic atoms may be isolated from an infinite 
cluster as, for example, atom B in Fig. 9. But if z is 
close to unity, so that nonmagnetic atoms are few, this 
factor is not so important because such isolation requires 
several nonmagnetic atoms to gather around one atom 
(four in the case of the plane lattice shown in Fig. 9). 
The probability of such an event is low if the number of 
nonmagnetic atoms is small. The second factor can thus 
be safely neglected and we can assume that the fraction of 
atoms belonging to the infinile cluster simply equals 
the fraction of magnetic atoms. We thus have P (z) = 2, 
provided 1 — z <1. 


„2. Each atom of a simple cubic lattice has six nearest 
neighbors located along the cube’s edges (see Fig. 12). 


The probability W, for all the nearest neighbors of an 
alom to be nonmagnetic equals the product of six proba- 
bilities: W, = (1 — x), The probability W of at least 
one of the atoms being a magnetic atom equals 
W(@)=1-W,=1-—(1 —2)8 

By virtue of formula (2) of Chapter 3, 

P, (z) = £W (z) = z [1 — (1 — z)°] 

so that if z< 1 
P, (x£) =% 62? 
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Obviously, 
P, (x) = x [1 — (1 — 2)7] 


for any lattice in which each atom has z nearest neighbors, 
so that for z <1 


Pa (x) = zz? 


3. Figure 51 shows 12 atoms in the neighborhood of 
atom 0. All of them can take part in forming a three- 
atom cluster. Such a cluster can be composed of, for 
instance, atoms J, 0, 2 if all these three atoms happen 
to be magnetic. The probability of this event equals the 
product of three probabilities: 2-a-+r = 23. The proba- 
bility for the cluster to be formed of atoms 0, 4, /2, or 
of any other trio of atoms, also equals 2°. 

First we have to answer the question about the number 
of such groups of three atoms. We begin with counting 


Fig. 51. 


how many of these groups include atoms 0 and 7. There 
are six of them: 015, 018, 019, 103, 102, 104. 

Now we consider three-atom groups including atoms 0 
and 3, but not atom Z. There are five of them: 036, 
03 11, 037, 203, 304. 

Likewise, there are four three-alom groups with atoms 
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5 2 10, 

9 and 2, but without atoms 1 and 3: 025, an (the 
204, and three-atom groups with atoms 0 and 4. u 
out atoms 7, 2, 3. These are 047, 048, 04 i P 

There are thus 6-5-4 4+ oe 18 pi or find 
groups, each having the probability 23. We want A 
the probability for at least one of them to ay Fone 
z < 1, the events consisting in the formation o aY ja, 
of the three-atom groups can be treated as TE: to 
Indeed, the probability for the groups 102 anc ets 
form simultaneously equals the Sagi a ppor aen 
0125, being magnetic, that is, the probabi ity a en 
tion of a four-atom cluster. The probability o eB lily 
equals a! = p.r < 23, Therefore, the Ane Ey Sey 
of these events is confirmed for z < 1. i ey ite b EA 
for at least one of the threc-atom Prous un orr 
equal to the sum of their probabilities, anc 


Py = 183 


i 0 
4. The probability for a randomly selected arom a 
belong to a cluster of not less than two aloms can Pn 
presented as the sum of the probabilities of incompati 
events: 


P, (2) = P; (x) + P, (x) (1) 


where P, (x) is the probability for an atom to belong to 
a cluster of two atoms. Therefore, 


Py (2) = P, (z) — Py (2) 2) 


The function P, (£) is given by formula (3) of Chapter 3, 
and hence, we only need to calculate Py (z). 
Atom 0 (sce Fi 


ig. 51) can form a two-atom cluster 
with atoms 7, 2, 3 or 4. 


The probability of the cluster 
being formed of atoms 0 and 7 equals the probability for 
both these atoms to be Magnetic times the probability 
for atoms 2, 3, 4, 8, 9, and 5 to be nonmagnetic, that is, it 
equals z? (1 — z). The events consis 
atom cluster being formed of atoms 0 
or 0 and 4 are essentially the same. 


ting in the two- 
and 2, or 0 and 3, 
All these events are 
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incompatible, and thus the probability Pa (x) equals the 
sum of four probabilities 
Pa (2) = 4a? (1 — x)" (3) 
Substituting formula (3) into formula (2), we obtain 
Py (2) = æ |1 — (L — 1)'] — 42? (1 — 2) (4) 
which is the solution of the problem. 
Making use of the binomial theorem, we can easily 
show that expression (4) contains no terms with power 


below 3. If «<4, then Pa (x) = 182, in accord with 
the result of the preceding exercise. 


Chapter 4 


1. 0.0085, 0.0072, 0.0051, 0.0026, 0.0006, 0.0000, 
0.0000...; 0.0067, 0.0044, 0.0019, 0.0003, 0.0000, 0.0000...; 
0.0032, 0.001, 0.0001, 0.0000, 0.0000... 

2. The number b is in fact an n-digit number. Con- 
sequently, b < 10”. In order to generate the next number 
b’, we have to find b°, divide it by 10", and take its inte- 
gral part. Therefore, b’ < b?/10”. But b2/10” = b-(b/10"). 
Furthermore, (b2/10") < b because (b/10") < 1. We thus 
conclude that b’ < b, which was to be proved. 

3. 5, 15, 5, 15, 5, 15... 

4. 5, 16, 9, 8, 5, 16, 9, 8, 5.... 

eB, 475 13, 4, Dy 1% 13, 4, Waes 

6. Condition (c) is nowhere sa tisfied. In Exercise 3, 
c = 0, in Exercise 5, condition (a) is not satisfied either, 
and so forth. 

7. Let Xo = 0. This gives us 0, 3, 4, 4, 2, 0, 3, 4, 4, 
2. and X, coincides with one of the numbers of 
this sequence whatever Xo is. 

_8. The distribution function of random numbers ob- 
tained from the generator is shown in Fig. 4. These num- 
bers form the array V. The fraction of nonblocked sites 
in the array V equals the fraction of random numbers in 
the array V, which satisly the inequality V < t. There- 
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fore, the mean fraction of nonblocked sites equals the 
probability for a random number to be less than ż. By de- 
finition of distribution function, this probability equals 
the area bounded by the curve f (y), the abscissa axis, 
and two perpendiculars erected at points 0 and ¢. In the 
present case this is the area of a rectangle, equal 


to t. Hence, the mean fraction of nonblocked sites x 
equals ż. 


Chapter 5 


1. When z is close to unity, 
an infinite cluster. Only those 


ster whose bonds connecting them with the surrounding 
system are all broken. The probability that one specific 
bond is broken equals 1 — z. If the lattice is square, it 
is necessary to interrupt four bonds originating from 
a site in order to isolate this site from the system (Fig. 52a). 
The probability of this event equals the product of the 
appropriate probabilities, that is, it equals (1 — a)‘. 
In order to isolate two sites from the system, six bonds 
must be interrupted (Fig. 52b). The robability of this 
event is ({ — z). If (1 —z)« 1, this probability is 
much smaller than that of one isolated site. Therefore, 
we can assume in the limiting case which is of interest 
now that all isolated sites are single, and the probability 


almost all sites belong to 
sites stay out of this clu- 


l 1 1 
| 


| l 
tw) ‘ ; (h) 


Fig. 52. (a) One isolated site; (b) two isolated sites. The 
intact bond is shown by the solid line, and the broken 
bonds by the dashed lines, 


for a randomly selected site to be isolated is (1 — z)*. 
The probabili 


1 ty for a randomly selected site to be non- 
isolated equals 4 — (1 — z)4, that is, in a square lattice 
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it is 

P? (z) = 1 — (1 —2)! 

Similar arguments show that in a triangular lattice 
P» (z) = 1 — (1 — 2)8 

and in a honeycomb lattice 


P? (z) = 1 — (1 — T}? 


These results hold if (1 — z) <1. 

2. The problem is solved by analogy to the preceding 
one. If (1 — z) < 1, almost all sites belong to an infinite 
cluster. A site chosen at random is isolated from the in- 
finite cluster if all its nearest neighbors are nonmagnetic 
atoms (for definiteness, we use the terminology of the 
ferromagnetic problem). As in the preceding problem, 
the probability for a site to have for one of its neighbors 
a magnetic atom isolated from the infinite cluster is low. 
Therefore, it is sufficient to calculate the probability 
for all neighbors of the given site to be nonmagnetic, The 
probability that a site is occupied by a nonmagnetic 
atom is 1 — x. The number of nearest neighbors equals 
the number of bonds that originate from a given site. 
Consequently, the results are the same as in Problem 1: 


in a square lattice: PS (z) = 1 — (1 — z)! 
in a triangular lattice: PS (z) = 1 — (1 — z)8 
in a honeycomb lattice: PS (x) = 1 — (1 — z) 


Hence, PS (z) = P? (z) for (1 — z) <1, in agreement 
with formula (2). 

3. Consider a honeycomb lattice with a fraction x of 
white bonds, and a triangular lattice with a fraction y 
of white bonds. (We remind the reader that the term 
“white bond” is synonymous to “unbroken bond”, and 
“black bond” is synonymous to “broken bond”. The phrase 
“a site is connected with another site” must be interpret- 
ed, unless expressly stated otherwise, as equivalent to 
“a site connected through unbroken, i.e. white, bonds”.) 
Let us superpose these lattices as shown in Fig. 53. Such 
sites as A, B, C are common for the two lattices, while 
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sites such as D belong only to the honeycomb lattice. 
The idea of the arguments that follow lies in that the 
percolation problem on a honeycomb lattice is reducible 
to the corresponding problem on a triangular lattice. 


Vig. 53. Star-della transformation. The dashed lines show 


the honeycomb lattice, and the solid lines show the trian- 
gular lattice. 


The probabilities for siles A, B, C to be connected 
with one another must be expressed in terms of the frac- 
tion z of white bonds in the E lattice. We need 
to make use of the geometry and statistical properties 
of the bonds that originate from D-type sites. After this 
we can discuss only the triangular lattice shown in Fig. 53 
by the solid lines and forget that cach triangle has the 
dashed lines and D-type sites. 

This operation, known as the star-delta transforma- 
tion, is frequently used in designing electric circuits. 

Tn fact, we shall need the following quantities: 

(1) W__ is the probability for site A not to be connect- 
ed with cither site B or site C. This probability equals 
the sum of the probabilities of two incompatible events. 
The first event consists in bond AD being black, with 
arbitrary colors of bonds BD and DC. Its probability 
equals 1 — x. The second event consists in bond AD 
teing white, and both bonds BD and DC being black. 
The probability of this event equals the product of three 


Answers and Solutions 253 


probabilities, < ad —a) (i —-2). Asa result, 
W..@=1—2- a a a) 


(2) W4- is the probability for site A to be connected 
with B but not with C. It equals the probability of bonds 
AD and DB being white, and bond DC being black, and 
is found as the product of the probabilities of all three 


events: 
Wi- = — 2% (2) 


(3) W 4 is the probability for site A to be connected 
with C but not with B. It is readily found that Wy = 


= Wy 

(4) Waa is the probability for site A to be connected 
with both B and C. 1t equals the probability for all the 
three bonds AD, DB, DC to be white and is calculated as 
the product of their probabilities: 


Waa (e) = (3) 


With these four probabilities known, we can leave 
aside the honeycomb lattice and solve our percolation 
problem on the triangu ir lattice. If this problem is solved, 
we shall have the critical value xy (H) for the honeycomb 
lattice. 

Of course, this transformation does not in the least 
facilitate the solution. However, the same probabilities 
W can be expressed in terms of the fraction y (it repre- 
sents the fraction of white bonds on the triangular lattice). 
At percolation threshold these probabilities assume quite 
definite values, yel unknown; however, equating the 
probabilities expressed as functions of x to the probabil- 
ities as functions of y, we can derive a relation between 
the thresholds of the honeycomb (ap (H)) and triangular 
(x (T)) lattices. 

Our next problem is thus to derive all four probabili- 
lies as functions of y- 

1) Wo. For site ‘A to be unconnected with both B 
and C, the color of bond BC is immaterial, but bonds AB 
and AC must be black. The probability of this event equals 
the product of probabilities: 


Wey) = (@— (4) 
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(2) W4._. Site A is connected with B but unconnected 
with C only if bond AB is white and both bonds BC and 
AC are black: 


We =o =g (5) 
(If e.g. bond BC were white, site A would connect with C 
via path ABC.) 

(3) W+ = W}, as in the preceding case. 

(4) W44. The probability for site A to be connected 
with both B and C equals the sum of the probabilities of 
four incompatible events. The first event consists in all 
the three bonds AB, BC, AC being white. Its probability 
is y®. The other three events consist in only one among 
the three bonds being black. For instance, if the black 
bond is AB, then site A is connected with C via a white 
hond AC, and with B via path ACB. The probability of 
each of the three events equals y? (1 - y). Finally, we 
obtain 
Was (y) = y+ 3y2(1—y) (6) 

At percolation threshold all W (x) must be equal to 
the respective W (y). This yields a system of equations 


W__(2)=W__(y); 1a +2 (1— z} = (1—y)? (7) 
Ws-(z)=We_(y); 2 (1—2)=y(1—y) (8) 
Waa (t)=Was(y)s =y 4 3y2(1—y) (9) 
These equations must he satisfied by z= xp (H) and 
y = xy (T). Furthermore, percolation thresholds satisfy 


relation (23) of Chapter 5, which yields that zp (H) = 
= 1 — zp (T). 


The substitution of 
y = ay (T), z = 1 — zp (T) 
into Eqs. (7), (8), and (9) transforms (8) into an identity, 
(L — ay (T))? zp (T) = ay (T) (1 = zp (T)}2 
and (7) and (9) into the identical cubic equation: 
xi (T)— 3x (T)-+1=0 
This equation has 


Ti quati a single root within the interval 0< 
S rp (T)< 1; 
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zp (T) = 2 sin (1/18) œ% 0.347 296 
Correspondingly, 
zp (H) = 1 — ap (T) = 1 — 2 sin (n/18) œ 0.652 704 


4. The first step is to find the area per site for each 
lattice shown in Fig. 16, under the condition that the 
spacing between nearest neighbors is a. 

Square lattice. Each unit cell includes four sites, but 
each of these sites belongs to four adjacent unit cells. 
Hence, there is one site per each square unit cell; in other 
woes the area per site is that of one square, that is, 
S (8) =a. 

Triangular lattice. Each unit cell includes three sites, 
but cach of these sites belongs to six adjacent unit cells. 
Hence, we find half a site per each triangular unit cell, 
and the area of two triangles per cach site. The area of 
an equilateral triangle with side a is y3a?/4. Therefore, 
the area per site is S (T) = V3a?/2. 

Honeycomb lattice. Each hexagon includes six sites, 
but each of these sites belongs to three adjacent unit cells. 
Hence, there is one hexagon per two sites. The area of 
a hexagon equals six times the area of an equilateral tri- 
angle with side a, that is, equals 3 V3a?/2. Hence, the 
area per each site is S (H) = 3 V3a?l4. 

The period a is determined in each lattice via a pe 
determined function a (z), with percolation threshold zp 
of the appropriate lattice used for x. This gives 


S(T) = (V 3/2) [a (0.35)? 
S (S) = [a (0.5)? 
S (11) = (3 374) [a (0.65)? 


Obviously, the function a (x) decreases with increasing 
r. (If the trees infect each other easier, the separation 
within a pair correspondingly decreases.) It is readily 
noticed, however, that this proposition is not sufficient 
in itself to derive even one of the above-written inequa- 
lities for areas. Indeed, in the honeycomb lattice the period 
a is the shortest, but the numerical coefficient in the 
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expression S (H) for area is the greatest, while in the 
triangular lattice the situation is reverse. The lattice 
with minimum area thus cannot be found unless we know 
the function a (z) in greater detail. 


Chapter 6 


1. The filling factors for plane lattices are easily 
calculated from the results of Exercise 4 of Chapter 5. 
We have illustrated this in the text for the honeycomb 
lattice. We shall, therefore, restrict the analysis to 3D 
lattices, or rather to two of them, leaving the remaining 
checks to the reader. x 

Simple cubic lattice. Each cubic unit cell (see Fig. 12) 
contains eight sites, but each of the sites is shared by 
eight adjacent unit cells. Hence, there is one site per 
one unit cell, and the volume per each site equals the 
volume of the cubic unit cell, i.e. a3. The spheres drawn 
around the sites have a radius a/2 and a volume 4na?/24. 
The fraction of volume occupied by the spheres equals 
the ratio of the volume of one sphere to the volume per 
site. Hence, in a simple cubic lattice 


f (SC) = 4n/24 ~ 0.52 


Body-centered cubic lattice. Each cubic unit cell 
(see Fig. 28b) contains nine sites, eight of which are at the 
vertices of the cube and one is at the center point. The 
site at the center belongs to one cubic unit cell exclusive- 
ly, while each of the vertex sites belongs to eight adjacent 
unit cells. Hence, there are two sites per one unit cell, 
and the volume per each site is half the volume of the 
cubic unit cell, i.e. «3/2. The nearest neighbor of each 
site is at a distance of half the body diagonal of the unit 
cell, i.e. V3a/2. The radius of a sphere drawn around each 
site is 34/4 and the volume of the sphere is (47 -39/2/3. 
+64) a = (n 73/16) a3. The filling factor equals the ra- 
tio at the volume of the sphere to the volume per site, 
so that 


{ (BCC) = a 3/8 = 0.68 
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Chapter 7 


1. The first coordination group of site 0 consists of 12 
type-I sites (Fig. 54), the second consists of six type-2 


Fig. 54. Neighbors of site 0 in the fee lattice. 


sites, and the third of 24 type-3 sites (uot all of them are 
shown in the figure); 12 +- 6 -+ 24 = 42. 


Chapter 8 


1. From formulas (3) and (5) we find af = 636 -10-8 cm 
= 636 Å, Ne= 7.810% cm. This is a very low 
critical concentration. The transition to metallic 
conduction takes place when there is one impurity 
atom per 108 atoms of the semiconductor matrix! A very 
sophisticated technique of purification is required to ma- 
nufaclure semiconductors with such a low impurity con- 
centration. 


Chapter 11 


1. Let us introduce a function @ (x), defining it as 
the probability for an information that has appeared al a 
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randomly chosen site to reach only a finite number of 
other sites. As before, P (z) = 1—Q. The number of 
independent channels through which the information 
emerges from each site is q. Let us determine the proba- 
bility for one of these channels to be interrupted at some 
step. The interruption can result from one of two incom- 
patible events, a or b: a consists in the first bond of the 
channel being broken, and b consists in the first bond being 
intact, but the site to which this bond leads being capable 
of passing the information to only a finite number of other 
sites. The probability of event a is 1 — z, and that of 
event b is xQ (x). The probabilities of incompatible events 
are added. Consequently, the probability for one channel 
to be interrupted equals 1 — x -+|- zQ (x). All the channels 
being independent, the probability for all of them to be 
interrupted is [1 — z -}- zQ (x)]7. 
This gives an equation for Q (x): 


Q (a) = [1 — z 4- 20 Ia 


A This equation reduces to Eq. (4) by the substitution 


Q' =1— az- rQ. Finally, we obtain 
p = E210) 2% 
q—1 


Note that percolation thresholds in the site and bond 
problems on the Bethe lattice are identical (zg = 1/q). 
In fact, this could be predicted at the start. Indeed, let us 
assume that the site to which the broken 


bond leads is 
blocked. This allows us to assume that all bonds are in- 
tact, so that 


$ our bond problem has been transformed to 
the site problem. Hence, these problems necessarily have 
identical percolation thresholds, 


Chapter 12 


1. We need to calculate the resistance of a cube with 


unit-long edge. The number of wires connected in parallel 
is, as before, 1/R?, but the length of each wire does not 
any more equal one unit of length. It is longer by a factor 
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equal to the ratio CEIR) = (z — re) 7675. Correspond- 
ingly, the resistance of one wire is not pọ but Po (£ — 
— ze) 767. The necessary result is therefore obtained 
by replacing the quantity pọ in formula (6) for os by 
Po (© — ze) 77V. The substitution gives 


o = 03 (2 — mT 


where 63 = pọ™l-2. Hence, t= b + v- 
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G. 1. Kopyloy, D.Sc. (Phys.-Math.) 


Elementary Kinematics of Elementary Particles 


This book tells a fascinating story of one of the basic 
goals of physics today: the discovery of the primary 
building blocks of matter. This field of science is called 
particle, or high-energy, physics, and is one of the 
frontiers of present-day physical research. 

llow is it possible to detect particles a hundred thou- 
sand times smaller than the atom, which itself is as 
many times smaller than an apple as the apple is 
smaller than the carth? How can we follow the motion 
of particles that have a velocity almost that of light? 
llow can we measure the lifetime of these particles 
when it is of the order of 0.00000000000000000000001 
second? What kind of a clock can we use? How can 
we investigate the properties of these astonishing and 
clusive bits of matter? 
AIL these and many other questions are comprehensive- 
ly answered in this book written for the layman by 
tħe late Dr. Gertsen Kopylov, who was a prominent 
scientist, well known in the world of particle physics. 
This book was written to be understood even by read- 
ers having only a secondary school education and it 
requires a knowledge of only elementary algebra and 
geometry. Nevertheless, in the author’s treatment, the 
material is in no way oversimplified or distorted. 


Ya. A. Smorodinsky 


Temperature 


This book starts with a historical background on the 
notion of temperature and the development of the tem- 
perature scale. Then Ya. A. Smorodinsky covers the 
fundamentals of thermodynamics and statistical phys- 
ics, only using concepts that will be familiar to high- 
school students, Having built a solid foundation, he 
exposes the reader to a number of phenomena that are 
essentially quantum-mechanical, but for which the con- 
cept of temperature “works”, and works very well. 
These include the spins in crystal lattices, inverse pop- 
ulation of energy levels, microwave background radia- 
tion, black holes, and cooling antiproton beams, 

Although it has been written for high-school students, 
the book contains a minimum amount of mathematics. 
Nevertheless, Ya. A. Smorodinsky compensates for this 
Severe restriction by the lucid manner in which he 
discusses very complicated effects, 


L. V. Tarasov 
A. N. Tarasova 


Discussions on Refraction of Light 


The rainbow and the Galilean telescope, the spectre of 
the Brocken and illuminated fountains—what can all 
these have in common? The answer is that in all of 
them the refraction of light has an important role to 
play. That is why you will read about them in this 
book written by the well-known popular science authors 
Lev and Aldina Tarasov. Professor Tarasov is a prom- 
inent figure in the field of quantum electronics and 
optics; Aldina Tarasova is a teacher of physics and 
studies methods for teaching the basics of science. 

The book treats refraction of light from two points of 
view. Firstly, it describes the many, and sometimes 
quite surprising, forms that the refraction of light can 
assume. Secondly, the book traces the long and diffi- 
cult path mankind covered before it came to under- 
stand some of the mysteries of Nature. Recent inven- 
tions such as the laser and the transmission of images 
over optic fibres are also dealt with. The book was 
first published in Russian in 1982, and has come 
out in English. We are certain it will be well received 
by those who are interested in the history of science 
and the present development of theoretical optics and 


its related technology as well as by all lovers of Na- 
ture. 


Th. Wolkenstein 


Electrons and Crystals 


The increasingly important field of solid-state physics 
concerns the behaviour of electrons in various crystals. 
Problems of solid-state physics, which include specific 
differences between metals and dielectrics and the re- 
markable properties of semiconductors, are particularly 
topical in today’s ‘electronic’ society. 

Electrons and Crystals by Dr. Theodore Wolkenstein 
covers the fundamentals of solid-state physics in an 
engaging way. Written in am easy, readable style, the 
book is intended as a supplement to textbooks in sec- 
ondary-school physics courses, and the approach to 
certain topics in the volume is, therefore, unique. The 
material is presented in terms of models and requires 
no special additional knowledge. Suitable for the gen- 
eral reader with a good command of elementary phys- 
ics and mathematics, this book can also serve as a 
useful study guide for high-school students, 


SCIENCE 
FORE'VERYONE 


This book is about percolation theory and its various 
applications, which occur mostly in physics and chemistry. 
The book is self-sufficient in that it contains chapters on 
elementary probability theory and Monte Carlo simulation. 
Most attention is paid to the relationship between the geo- 
metrical and physical properties of systems in the vicinity 
of their percolation thresholds, The theory is applied to 
examples of impurity semiconductors and doped ferromag- 
netics, which demonstrate its universality. Although writ- 
ten Tor students at high schools, the book is very good 
reading for college students and will satisfy the curiosity 


‘of a physicist for whom this will be a first encounter with 


percolation theory. 


